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EXECUTIVE  SUMMARY 


The  Workshop  on  the  Measurement  of  Transient  Pressure  and 
Temperature  was  held  at  the  National  Institute  of  Standards  and 
Technology  (NIST) , Gaithersburg,  MD,  on  23-24  April  1991.  There 
were  77  attendees,  55%  from  government  and  45%  from  the  private 
sector.  Aerospace,  automotive,  and  instrumentation  manufacturers 
were  present  along  with  engineers  from  private  sector  and 
government  (Army,  Navy,  Air  Force,  and  NASA)  research 
organizations.  The  Army,  Navy,  and  NASA  primary  standards 
laboratories  also  were  represented.  The  workshop  drew 
international  attention  with  participants  from  Taiwan,  Sweden, 
Austria,  and  France.  Only  France  currently  offers  transient 
pressure  calibrations  as  part  of  their  national  measurement  system. 

The  purpose  of  the  workshop  was  to  identify  current  and  emerging 
measurement  problems  and  needs  for  transient  pressure  and 
temperature  measurements,  to  explore  scientific  and  technical 
barriers  and  opportunities  for  achieving  measurement  accuracy 
goals,  and  to  assist  in  the  formulation  of  an  appropriate  role  for 
NIST. 

The  scope  of  this  workshop  was  the  measurement  of  transient 
pressure  and  transient  temperature;  however,  the  participants  were 
much  more  concerned  with  pressure  than  with  temperature. 
Temperature  was  discussed  briefly  by  two  speakers:  Damion  described 
the  use  of  a shock  tube  to  provide  a transient  temperature  stimulus 
to  study  a special  platinum  resistance  thermometer  and  Rosasco 
described  how  spectroscopic  techniques  can  be  used  to  measure  both 
transient  temperature  and  pressure. 

I.  Measurement  Problems  and  Needs 

The  range  of  industrial  and  scientific  interests  represented  at  the 
workshop  is  a reflection  of  the  importance  of  transient 
measurements  in  a wide  range  of  applications,  such  as  aeronautics, 
blast-wave  analysis,  internal  combustion  and  jet  engine  testing, 
ballistics,  health,  and  safety  (see  for  example  Winiarz,  Ball, 
Estes  and  Smith,  Hartenbaum  and  Sturtevant)  . The  accuracy  of  these 
transient  measurements  is  crucial,  and  a common  thread  in  many  of 
the  presentations  was  the  need  for  traceability  to  national 
standards  at  NIST.  In  many  cases,  agency  policy  and  government 
regulations  require  such  traceability  (Ball,  Estes  and  Smith). 
There  currently  are  no  national  primary  standards  for  the 
measurement  of  either  transient  pressure  or  transient  temperature 
in  the  metrology  system  of  the  United  States,  and  hence  there  is  no 
traceability  to  NIST. 

As  one  might  expect  from  the  range  of  applications,  the  amplitude 
and  frequency  of  transient  pressure  calibration  needs  are  widely 
varied  and  range,  for  example,  from  13  Pa  at  50  kHz  for  the 
national  aerospace  plane  to  1 GPa  at  100  MHz  for  studies  on 
explosives.  Estes  and  Smith  documented  the  requirements  of  NASA. 
We  present  these  in  Fig.  1 as  an  outline  in  pressure  and  frequency 
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space,  with  the  numbers  inside  the  outline  being  the  desired 
accuracies  in  that  pressure-frequency  region,  expressed  as 
percents.  At  the  conclusion  of  the  workshop,  the  participants  were 
invited  to  indicate  the  pressure  and  frequency  ranges  of  the 
calibration  services  they  needed  with  the  corresponding  accuracies. 
Figure  2 is  a collage  of  the  replies  from  organizations  other  than 
NASA.  Again,  the  numbers  inside  the  outline  represent  the  needed 
accuracies  expressed  as  percents. 

II.  Current  Practice  in  Transient  Pressure  Calibrations 

Because  there  is  no  traceability  to  NIST  for  transient  pressure  or 
temperature  measurements,  laboratories  that  need  to  make  such 
measurements  are  left  to  their  own  resources  to  develop  calibration 
and  testing  programs.  Calibration  generally  is  understood  to  mean 
"...  a test  during  which  known  values  of  measurand  are  applied  to 
a transducer  and  corresponding  output  readings  recorded"  (see  Lally 
who  quotes  A Guide  for  the  Dynamic  Calibration  of  Pressure 
Transducers,  ANSI  MC88. 1-1972,  published  by  The  American  Society  of 
Mechanical  Engineers,  New  York,  NY)  . Because  of  the  technical 
problems  associated  with  providing  a known  measurand  for  transient 
pressure,  an  approach  often  taken  is  to  assume  a calibration  with 
static  pressure  standards  and  methods  to  be  valid  for  pressure 
transducers  used  under  transient  pressure  conditions.  As  widely 
documented  by  the  speakers,  this  assumption  is  incorrect  because 
the  pressure  transducer  is  a mechanical  structure  which  will  behave 
differently  at  different  frequencies.  The  static  calibration 
yields  no  frequency  behavior  information. 

The  ranges  of  the  needed  calibration  services  indicated  in  Figs.  1 
and  2 span  more  than  eight  orders  of  magnitude  in  frequency  and  six 
orders  of  magnitude  in  pressure.  It  is  not  surprising  that 
multiple  techniques  are  required  to  meet  the  needs.  The  current 
state-of-the-art  in  transient  pressure  calibrations  was  presented 
by  a number  of  speakers.  From  the  discussions,  we  observe  that  as 
the  frequency  and  pressure  increase,  the  difficulty  of  generating 
a known  measurand  for  the  calibration  of  transducers  also 
increases . 

A variety  of  transient  pressure  generators  were  discussed  (Lally, 
Newhall,  Riegebauer,  Damion)  including  the  pistonphone,  deadweight 
tester  with  a quick-opening  valve,  pneumatic  and  hydraulic  quick- 
opening valve  devices,  a column  of  liquid  mounted  on  a shaker  head, 
a modified  air  compressor,  a weight  falling  on  a piston  which 
compresses  a confined  liquid  (drop-hammer) , and  the  shock  tube.  For 
the  very  high  frequency  regime,  i.e.  for  nanosecond  rise-time 
devices,  Hartenbaum  and  Sturtevant  suggested  the  use  of  an 
electrically  driven  flyer  plate  to  produce  the  requisite  (i.e.  in 
frequency  and  amplitude)  shock  waves. 

These  generators  may  be  classified  as  periodic  or  aperiodic. 
Periodic  generators  yield  sinusoidally  varying  pressures,  are 
limited  to  lower  frequencies,  and  suffer  decreased  amplitude  as  the 
frequency  increases.  The  pressure  amplitudes  are  generally  not 
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known  by  absolute  means  and  thus  periodic  generators  are  used  as 
comparators  rather  than  as  primary  standards  (Damion) . The  quick- 
opening valves,  the  drop  hammer,  and  the  shock  tube  are  aperiodic 
generators  and  were  considered  in  greater  detail. 

Aperiodic  generators  yield  a pressure  step  or  a pressure  impulse. 
These  generators  subject  the  transducer  to  a range  of  frequencies, 
the  high  frequency  limit  of  the  generator  usually  depends  upon  the 
rise  time  of  the  pressure;  the  low  frequency  limit  depends  upon  the 
duration  or  width  of  the  pressure  stimulus.  The  description  of  the 
complete  frequency  content  of  the  stimulus  depends  on  the  details 
of  the  temporal  evolution  of  the  step/impulse  which,  in  general, 
are  not  known.  We  observe  from  many  of  the  presentations  that  it 
is  difficult  to  assess  the  sources  of  the  high  frequency  content  in 
the  response  of  a transducer  under  calibration.  Some  oscillations 
observed  with  a quick-opening  valve  step-generator  are  thought  to 
be  a common  characteristic  of  such  devices  (Newhall) . Defective 
transducers  or  defects  in  the  electrical  cables  and  connections 
also  can  give  rise  to  strong  oscillations  in  the  step  or  impulse 
response  (Lally,  Newhall) . Shock  waves  also  have  high  frequency 
content  which  is  not  well  understood.  Under  these  conditions,  the 
frequency  information  derived  from  the  response  of  a transducer 
used  with  such  generators  is  essentially  qualitative. 

The  steady-state  amplitudes  for  the  step  generator,  or  peak 
amplitudes  for  the  impulse  generator,  are  used  to  determine  the 
transducer  sensitivity  (the  ratio  of  the  transducer  output  voltage 
to  unit  pressure) . These  amplitudes  are  derived  either  from  direct 
measurement  of  steady-state  pressures  in  quick-opening  valve 
devices  (Newhall,  Lally,  Damion)  or  from  calculations,  for  example, 
ideal  shock-wave  theory  (Lally,  Damion)  or  a quasi-empirical  model 
for  the  drop-hammer  generator  (Riegebauer) . In  the  latter  paper, 
two  different  approaches  are  used  to  arrive  at  the  peak  pressure 
and  good  agreement,  within  4%,  is  demonstrated. 

France  is  the  only  nation  currently  offering  transient  pressure 
calibration  services  as  part  of  their  national  measurement  system. 
The  laboratory  is  located  at  Ecole  Nationale  Superieure  d'Arts  et 
Metiers  (ENSAM)  in  Paris  and  is  under  the  direction  of  J.  P. 
Damion.  The  service  covers  the  frequency  range  from  a few  Hz  to  a 
few  tens  of  kHz  at  pressures  up  to  2 0 MPa.  Quick-opening  valve 
devices  are  used  for  the  lower  frequencies  and  shock  tubes  for  the 
higher  frequencies  (Damion).  Damion's  discussion  gives  good 
account  of  the  practical  limits  and  assumptions  implicit  in  the  use 
of  the  various  calibration  devices. 

III.  Pressure  Transducers 

Many  of  the  presentations  and  the  discussions  following  them 
indicated  that  modern  transducer  technology,  based  on  a wide 
variety  of  physical  phenomena,  is  highly  developed.  Transducer 
resonance  frequencies  greater  than  100  kHz  are  routinely  achieved, 
and  thus  measurements  for  frequencies  up  to  10 's  of  kilohertz  can 
be  made.  Lally  presented  sensitivity  data  for  a transducer 
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obtained  by  5 different  techniques:  the  dead-weight  tester,  a 
hydraulic  drop-hammer,  two  varieties  of  quick-opening  valve 
devices,  and  a shock-tube.  Excellent  consistency,  at  the  2%  level, 
was  demonstrated  for  the  sensitivity  of  this  transducer,  which  has 
a resonance  frequency  in  the  range  of  500  kHz.  Newhall  compared  a 
piezoelectric  gage  and  a strain-type  gage  subjected  to  a 700  MPa 
pressure  step  in  a quick-opening  valve  device  and  reported 
virtually  identical  responses.  Riegebauer, demonstrated  transducer 
repeatability  at  the  0.5%  level  for  a drop-hammer  device  with  a 
rise  time  in  the  1-2  ms  regime.  Perhaps  the  most  widely  tested 
transducers  are  the  Navy's  tourmaline-based  piezoelectric  devices. 

The  Navy's  underwater  explosion  shock  measurements  program  has  used 
gages  made  of  single-crystal  tourmaline  for  over  40  years. 
Tourmaline  is  bulk  or  hydrostatically  sensitive.  It  does  not 
require  any  mechanical  shielding  to  permit  only  a certain  direction 
of  strain  as  does  quartz,  thus  avoiding  the  complexities  of 
diaphragms,  mechanical  resonances,  and  reflections.  The  Navy 
calibrates  the  gage  statically  by  means  of  a dead  weight  tester  and 
assesses  the  time  response  with  a quick-opening  valve.  The  gage 
calibration  constant  is  stable  and  linear  over  the  span  of  0.1  to 
350  MPa.  Many  gages  used  and  recalibrated  over  a 15  year  period 
show  no  aging  effects.  The  gage  shows  no  signs  of  hysteresis.  The 
tourmaline  gage  remains  the  Navy's  gage  of  choice  (Tussing) . 

IV.  Signal  Transfer  and  Processing 

Modern  high  speed  digital  data  recording  systems  are  routinely  used 
for  transient  measurements,  and  indications  are  that  they  are  quite 
adequate,  in  terms  of  speed,  accuracy,  and  data  processing 
capability,  for  the  needs  of  the  community.  In  any  application 
there  are  important  issues  with  respect  to  spurious  signal 
generation  and  pick-up  and  with  the  noise  associated  with  the 
measurement  environment;  these  issues  were  addressed  by  many 
authors  (Winiarz;  Ball;  Hartenbaum,  Morrison  and  Gord;  Stein). 
Such  signal/noise  problems  are  ubiquitous.  Some  are  associated 
with  the  transducer,  for  example,  the  effects  of  temperature 
transients  are  familiar.  Stein  also  cautioned  us  to  be  wary  of 
strain-induced  voltages  in  impedance-based  measurement  devices. 
Other  problems,  perhaps  the  most  pervasive,  are  associated  with  the 
transfer  of  the  electrical  signal  to  the  recording  device. 
Hartenbaum  et  al.  described  advances  in  noise  reduction  by  means  of 
a new  signal  conditioner  based  on  hybrid  optical-electronic 
circuitry.  These  types  of  problems  must  be  considered  in  any 
experimental  design,  including  that  of  a calibration  facility.  In 
general,  care  must  be  exercised  in  analyzing  and  interpreting  the 
frequency  content  of  one's  signal,  especially  the  high  frequency 
components. 

V.  Technical  Opportunities  and  Barriers 

In  assessing  the  current  state  of  the  national  calibration  system 
for  transient  pressure,  we  note  first  that  there  are  dozens  of 
laboratories  active  in  this  area.  Most  of  these  have  transient 
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sources  which  are  used  mainly  for  qualitative  assessment  of 
transducer  response  functions,  e.g.,  no  unexpected  resonance 
behavior,  lag,  overshoot,  or  decay.  Some  laboratories  have 
transient  sources  where  the  amplitude  of  the  step  stimulus  or  of 
the  peak  impulse  stimulus  is  known. 

We  also  observe  that  transient  pressure  transducer  technology  is 
quite  reliable.  This  implies  that,  once  calibrated,  transducers 
can  be  used  to  transfer  national  standards  to  calibration 
laboratories  and  to  maintain  the  accuracy  of  a national  calibration 
system.  These  observations  provide  a technical  basis  for 
establishing  a national  calibration  system  in  order  to  provide  the 
needed  traceability  for  transient  pressure  measurements. 

At  the  present  time,  no  transient  generator  is  sufficiently  well 
understood  to  allow  it  to  serve  as  a primary  standard.  Assurance 
of  the  accuracy  of  our  knowledge  of  the  pressure-amplitude  versus 
frequency  for  a transient  source  becomes  increasingly  difficult  for 
frequencies  above  1 kHz.  The  pressure  wave  in  a shock  tube  (the 
highest  frequency  source  for  gas  phase  calibrations)  can  only  be 
assumed  to  follow  ideal  shock  wave  theory.  Detailed  information  on 
the  rise  time,  shock  front  curvature,  shock  attenuation,  or  other 
non-ideal  behavior,  although  of  concern,  has  not  yet  been  obtained 
for  calibration  sources  (Damion) . At  the  present  time  we  cannot 
resolve  this  issue  by  recourse  to  a first  principles  theory  for  the 
transient  phenomenon.  Furthermore,  because  of  the  high  frequency 
content,  conventional  transducers  cannot  be  used  to  answer  the 
questions  experimentally.  This  situation  limits  the  ability  to 
characterize,  even  qualitatively,  the  transducer  response  function 
above  a few  10 's  of  kHz. 

These  observations  underlie  the  approach  NIST  is  following  to 
produce  a primary  standard  for  the  measurement  of  transient 
pressure  and  temperature  (Rosasco) . In  the  NIST  approach,  the 
primary  standard  is  the  pressure  and  temperature  sensor,  a sample 
of  diatomic  gas  which  is  part  of  the  pressurizing  gas.  Diatomic 
gas  molecules  have  a fundamental  vibrational  motion  whose  frequency 
is  a function  of  pressure.  These  molecules  also  have  well  defined 
rotational  energy  levels  whose  populations  provide  a reliable 
measure  of  temperature.  The  vibrational  frequencies  and  the 
populations  of  the  rotational  energy  levels  can  be  determined  by 
laser  spectroscopy.  The  pressure  standard  can  ultimately  be 
regarded  as  primary  in  that  the  shift  and  widths  of  the  spectral 
features  used  to  determine  the  pressure  can  be  obtained  from  ab 
initio  theory.  The  temperature  measurement  approach  is  also 
primary  in  that  it  is  based  on  the  Boltzmann  equation  for  the 
thermodynamic  system  in  local  equilibrium. 

The  gas  molecule  transducer  offers  significant  advantages  over 
conventional  transducers.  It  does  not  contaminate  the  pressure 
information  through  resonance,  phase  lag,  ringing,  damping,  limited 
frequency  response,  and  limited  rise  time  as  is  the  case  with 
conventional  transducers.  The  signals  carrying  the  pressure  and 
temperature  information  are  transferred  to  the  recording  device  by 
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means  of  laser  beams;  the  traditional  sources  of  noise  and  pick-up 
are  thereby  eliminated.  Significantly,  this  transducer  provides 
simultaneous  pressure  and  temperature  information  which  is 
necessary  to  completely  specify  a transient  source. 

In  order  to  be  useful  for  calibration  of  conventional  transducers, 
the  proposed  primary  standard  needs  to  be  realized  in  a transient 
source.  The  source  to  be  used  is  a shock  tube  capable  of 
generating  reflected  shock  pressures  up  to  20  MPa.  The  shock  tube 
will  be  calibrated  by  means  of  the  molecular  transducer  to  avoid 
the  necessity  of  assuming  ideal  shock  wave  behavior  for 
calculating  the  pressure  and  temperature  from  shock  tube  theory. 
The  calibrated  shock  tube  will  serve  as  a reference  standard  for 
the  calibration  of  conventional  pressure  transducers  and  will 
provide  a method  for  decoupling  pressure  and  temperature  effects  in 
the  transducer  response. 

VI.  Suggestions  for  the  Role  of  NIST 

There  was  considerable  discussion  on  the  role  of  NIST  in  solving 
these  measurement  problems.  The  consensus  view  which  emerged  from 
the  discussions  oan  be  summarized  as  follows: 

1.  NIST  should  develop  standards  to  provide  traceability  for  the 
calibration  of  transient  pressure  transducers. 

2.  The  participants  were  generally  willing  to  collaborate  on 
these  measurement  problems  through  a working  group  organized  to 
provide  a forum  for  discussion,  and  through  measurement  round- 
robins  for  transient  pressure.  More  than  15  laboratories  indicated 
a willingness  to  be  involved  in  such  round-robins. 

3.  The  shock  tube  is  the  instrument  of  choice  for  obtaining 
pressure  transducer  transfer  functions  because  of  the  fast  rise 
time  of  the  pressure  step. 

4.  There  was  genuine  interest  in  the  proposal  by  NIST  to  develop 
a national  standard  based  on  the  transfer  of  a molecular-level 
primary  standard  to  a reference  shock  tube. 


V.  E.  Bean 
G.  J.  Rosasco 
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Figure  1.  The  outlined  areas  in  pressure- frequency  space  represent 
the  transient  pressure  calibration  needs  of  NASA.  The  numbers  in 
the  areas  represent  the  needed  accuracies  expressed  in  percent. 


7 


cd 

CU 

03 

(72 

CO 

O 

ia£) 

O 


12 

11 

10 

9 

8 

7 

6 

5 

4. 


Dynamic  Pressure  Calibration  Needs 


-9  -8  -7  -6  -5  -4  -3  -2  -1 

log  rise  time  (s) 


■ ; 

0^ 

/ 

) 3 

25 

r 

1 

1 

I 

5 

0.. 

2 

o 

0.5 

7 

6 

5 

4 

3 

2 

1 

0 

-1 


P(psi) 


147,000 


Cd 


a; 

1-^ 

(72 

(72 

02 

ti£) 

O 


0 1 


1,470 


14.7 


Figure  2.  The  outlined  areas  in  pressure-frequency  space  represent 
the  transient  pressure  calibration  needs  of  organizations  with 
representatives  attending  the  workshop  other  than  NASA.  The 
numbers  in  the  areas  represent  the  needed  accuracies  expressed  in 
percent . 
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PREFACE 


Twelve  technical  talks  were  presented  during  the  workshop. 
Speakers  were  encouraged,  but  not  required,  to  publish  in  these 
proceedings.  All  twelve  talks  are  here  represented  in  one  of  three 
formats;  traditional  conference  papers,  abstracts  with  view  graphs, 
or  reprints  from  archival  sources.  The  last  names  in  parentheses 
in  this  summary  identify  the  speaker's  published  material  herein 
where  more  information  on  the  topic  can  be  found. 
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ABSTRACT 


To  be  presented  at  the  Workshop  on  the  Measurement  of 
Transient  Pressure  and  Temperature  at  the  National 
Institute  of  Standards  and  Technology  (NIST),  on  the 
23-24  Apri 1,  1991. 


Dynam i c High  Pressure  Measurement 
I ssues  and  I nstrumentat i on 


Marek  L.  Winiarz,  P.  E. 


The  world  of  dynamic,  transient  pressure 
measurements  at  high  regimes  is  full  of  assumptions 
and  uncertainties.  First  section  of  the  paper  will 
address  the  key  issues  of  metrology  including 
validation,  calibration,  and  sensor  inclusion. 
Attempts  at  solutions  such  as  calculat ional  modeling, 
calibration  facilities  will  be  illustrated. 

The  second  part  will  review  instrumentation 
(sensors)  used  to  measure  high  pressure  and  stress 
dynamically.  Gage  examples  will  include  strain-based 
devices,  piezo-effect  devices,  and  others.  Mounting 
issues,  inclusion  problems  and  data  acquisition  will 
be  addressed  for  each  gage. 


The  concluding  section  will  provide  an  overview 
of  the  state-of-the-art  techniques  used  to  measure 
pressure  in  turbine  engine  design  and  development  at 
General  Electric  Aircraft  Engines.  It  will  illustrate 
how  engineers  approach  these  critical  measurements  in 
the  presence  of  temperature  and  air-flow  fluctuations. 
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GE  Aircraft  Engines 

Quality  Technology  Center,  Metrology  Subsection 
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DATA  PRESENTATION 
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Example  of  Zero  Shift  Correction 
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STATIC  AND  STAGNATION  MEASUREMENT 
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Typical  Installation  at  a Test  Site 
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36 


37 


PIEZOELECTRIC  MEASUREMENTS 
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TOURMALINE  GAUGE  FOR  NAVY'S  UNDERWATER 


EXPLOSION  SHOCK  MEASUREMENTS 
(UNCLASSIFIED) 

RONALD  B.  TUSSING 

NAVAL  SURFACE  WARFARE  CENTER 
WHITE  OAK  LABORATORY 
(301)  394-1187 

1.  INTRODUCTION 

There  are  many  considerations  to  be  made  in  selecting  a transducer  or 
piezoelectric  gauge  for  underwater  shock  wave  phenomena  measurements.  The 
accuracy  of  the  total  measurements  system  is  the  primary  factor,  of  which  the  gauge  is 
only  a part.  However,  the  total  system  can  be  no  better  than  the  gauge  or  transducer 
selected. 

This  paper  addresses  the  characteristics  of  tourmaline,  the  Naval  Surface 
Warfare  Center  or  NAVSWC  (Navy)  tourmaline  gauge  and  its  response,  and  the 
methods  used  to  determine  the  response.  Tourmaline  has  been  the  transducer  of 
choice  for  over  40  years,  and  still  remains  such.  The  NAVSWC  gauge  has  been  used 
to  measure  the  data  to  characterize  all  of  the  Navy's  explosives  for  over  20  years, 
which  are  archived  at  the  NAVSWC  White  Oak  Laboratory.  The  NAVSWC  gauge  has 
been  made  available  for  data  comparisons  internationally  to  over  a dozen  countries, 
as  well  as  to  national  laboratories  and  companies.  The  NAVSWC  gauge  has  been 
used  on  innumerable  tests:  charge  weights  ranging  from  a few  grams  to  20  tons;  test 
sites  including  pressure  vessels,  ponds,  quarries,  rivers,  and  the  deep  ocean;  and 
environments  ranging  from  the  tropics  to  the  arctic. 

2.  NAVSWC  TOURMALINE  GAUGE 
2.1  Characteristics  of  Tourmaline 

Tourmaline  is  a hard  natural  crystal  that  is  bulk  or  hydrostatically  sensitive.  Its 
piezoelectric  constant  is  high  enough  to  be  useful  and  is  constant  over  an  adequate 
temperature  range. 
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A survey  of  piezoelectric  activity  of  minerals  by  W.  L Bond  (1948)  revealed  that 
only  quartz  and  tourmaline  were  practical  for  extensive  piezoelectric  work.  Since  then, 
lithium  niobate  and  PVDF  or  PVF2  have  shown  some  promise,  though  PVDF  is 
presently  characterized  well  in  only  one  dimension  and  its  response  Is  nonlinear. 
Some  of  the  characteristics  for  tourmaline  are: 

Piezoelectric  Constant:  K = 1 1 Picocoulombs/psi 
Linear  Pressure  Range:  1 psi  - 50,000  psi 
Dielectric  Constant:  7.44 

Density:  190  - 200  Ib/ft^  (3.0  - 3.2  g/cm^) 

Hardness:  7.3  (Mohs'  Scale) 

Sound  Speed:  19,685  ft/  sec  (6000  m/sec) 

Piezoelectric  Moduli:  d31,d33 

Tourmaline  must  be  cut  perpendicular  to  the  z axis  or  optic  axis.  Charge  is  generated 
on  the  faces  perpendicular  to  this  axis  by  pressure  applied  to  these  faces  as  well  as 
the  parallel  faces  or  edges:  K = d33  + 2d31 . The  properties  of  quartz  are  similar  to 
tourmaline  except  that  quartz  is  not  hydrostatically  sensitive,  and  its  edges  must  be 
constrained  or  protected  from  the  shock  wave.  Tourmaline  has  the  hydrostatic 
sensitivity  advantage  and  does  not  require  any  mechanical  shielding  to  permit  a 
certain  direction  of  strain.  This  avoids  the  accompanying  problems  of  diaphragms, 
mechanical  resonances,  reflections  and  complexity. 

2.2  History  of  the  NAVSWC  Gauge 

The  NAVSWC  (Navy)  gauge  is  fabricated  and  calibrated  currently  at  the  R15 
Dahlgren  site.  The  gauge  evolved  from  the  Wood's  Hole  gauge  used  by  R.  H.  Cole.  A 
contingent  of  scientists  from  Wood's  Hole  involved  with  underwater  explosions 
research  moved  to  the  Naval  Gun  Factory  in  Washington,  D.  C.  in  1949,  and  then  to 
Indian  Head,  Maryland  in  about  1950.  Some  of  this  group  moved  to  the  Naval 
Ordnance  Laboratory  (presently  NAVSWC)  in  White  Oak  Maryland  when  it  opened; 
however,  all  testing  was  done  in  the  Potomac  River  at  Indian  Head  site.  The  best 
coating  used  at  this  time  was  bees-wax  and  it  became  the  standard  internationally. 
Many  gauges  were  purchased  from  a small  company  Crystal  Research,  although 
gauges  were  always  fabricated  in-house,  and  all  gauges  were  calibrated  in-house 
using  a dead  weight  tester.  When  Crystal  Research  went  out  of  business  in  the 
1960's,  the  Naval  Ordnance  Laboratory  bought  all  their  equipment  and  tourmaline. 
Some  of  the  tourmaline  went  to  UERD(Underwater  Explosions  Research  Division)  of 
DTRC(David  Taylor  Research  Center).  When  their  gauge  person  retired,  NAVSWC 
was  the  only  one  in  the  U.  S.  fabricating  an  underwater  tourmaline  gauge.  The 
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NAVSWC  gauge  was  then  made  available  to  other  agencies  and  internationally  to 
other  countries  as  a test  standard  for  comparing  data,  and  it  still  is  made  available. 

In  the  mid  1 960's  R.  S.  Price  and  J.  B.  Dempsey  experimented  with  gauge 
coatings  to  improve  the  variations  in  testing,  using  a statistical  approach  to  minimize 
the  standard  deviation  or  sigma.  The  myriad  of  coatings  tried,  included  various 
plastics,  waxes,  epoxies  and  various  rubber  compounds.  In  1967,  a gauge  immersed 
in  silicone  oil  in  a plastic  boot  proved  to  diminish  the  variations  noted  by  all  of  the  other 
coatings.  Final  and  conclusive  results  were  reported  in  1972.  The  problem  with  other 
coatings  seemed  to  be  that  each  generated  a charge  in  an  inconsistent  way  that 
added  to  the  charge  generated  by  the  tourmaline  crystals.  These  materials  essentially 
changed  the  values  of  the  gauge  constants,  inconsistently  causing  variations  with  use 
and  time.  Other  materials  rounded  the  shock  peaks  drastically.  In  the  early  1970’s,  R. 
B.  Tussing  demonstrated  a transient  method  of  evaluating  a gauge's  transfer  function 
which  was  more  definitive  and  required  fewer  tests,  than  the  statistical  approach.  The 
technique  was  not  reported  upon  until  1982. 

2.3  NAVSWC  Gauge  Construction 

The  tourmaline  must  be  sliced  perpendicular  to  the  z axis  and  then  the  pieces  are 
lapped  (see  Fig.  1 ).  The  discs  are  cut  with  a diamond  hole  cutter  and  assembled 
between  metal  tabs.  Disc  sizes  range  from  1/8"  to  2",  with  one  to  eight  discs 
paralleled.  The  four  disc  gauge  is  used  more  in  sizes  of  1/4",  3/8",  1/2",  3/4",  and  1". 
Single  disc  gauges  are  used  for  close-in  measurements  at  high  pressures  where 
gauge  damage  is  expected  and  high  output  is  not  desired.  The  assembled  crystal  disc 
elements  are  enameled  with  a high  dielectric  paint  and  mounted  to  an  oil  barrier  of 
nylon,  after  the  element  has  been  calibrated  using  the  dead  weight  tester  which 
employs  a quick-opening  valve.  The  gauge  element  and  oil  barrier  is  then  immersed 
in  Dow  Corning  DC200  high  dielectric  silicone  fluid  (100  centistokes)  in  a Tygon  boot 
heated  and  pinched  closed  at  one  end  (See  Fig.  2).  The  flattened  end  was  found 
experimentally  to  cause  less  flow  perturbations  than  a hemispherically  molded  end.  A 
lighter  viscosity  oil  (1000  m/sec)  is  used  than  that  which  would  match  the  sound 
velocity  in  water.  Though  the  sound  velocity  would  match  better  with  a heavier  fluid, 
the  bubbles  are  nearly  impossible  to  remove,  particularly  in  the  field.  Air  bubbles  in 
the  fluid  cause  reflections  and  result  in  erroneous  data. 

2.4  Characteristics  of  the  Oil-booted  Gauge 

The  gauge  calibration  constant  KA  (pico-coulombs/psi)  is  a constant:  it  is  linear 
and  stable  from  about  1 psi  to  50,000  psi.  This  has  been  confirmed  by  calibration  tests 
at  NAVSWC  and  reported  in  the  literature.  It  is  reasonable  to  assume  that  tourmaline 
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is  linear  over  a much  wider  range.  This  seems  to  be  confirmed  by  the  excellent 
agreement  in  experimental  results  obtained  in  controlled  explosion  tests  at  NAVSWC. 

Many  gauges  have  been  used  and  re-calibrated  for  a span  of  15  years  showing 
no  aging  effects  with  a standard  deviation  of  less  than  1 percent.  Pyro-electric  effects 
have  been  reported  to  be  negligible  for  underwater  tests  for  the  signal  durations 
encountered,  even  for  large  charges. 

The  booted  tourmaline  gauge  shows  no  signs  of  hysteresis.  The  upper  frequency 
response  is  dependent  upon  the  gauge  size  or  transit  time  of  the  shock  wave  across 
the  gauge,  while  the  low  frequency  is  controlled  by  the  input  time  constant  of  the 
system  to  which  the  gauge  is  connected.  With  input  impedances  of  100  to  300 
Megohms,  and  cable  lengths  normally  required  by  the  safe  distance  from  the  charges, 
the  input  time  constant  is  100  or  more  times  the  duration  of  the  shock,  usually  yielding 
a lower  response  of  less  than  0.1  HZ.  The  gauge  rise  or  transit  time  is  considered  in 
detail  In  the  next  section. 

Adequate  sensitivity  with  good  signal-to-noise  (S/N)  ratio  and  an  adequately 
small  transit  time  can  generally  be  achieved.  Sometimes  it  is  necessary  to  place  a 
preamplifier  near  the  gauge  to  meet  these  conditions. 

The  gauge’s  ruggedness  has  been  proven  in  innumerable  tests,  in  some  tests  to 
pressures  of  50,000  psi,  others  where  the  gauge  was  attached  to  a model's  steel  hull 
that  was  damaged  and  torn,  and  on  ship  shock  tests  attached  to  the  sail  of  an 
operational  submarine. 

Minimum  distortion  includes  a great  deal.  Any  measuring  device  of  finite  size  will 
disturb  the  shock  medium.  Some  initial  diffraction  of  the  shock  front  occurs  as  it 
passes  over  the  gauge.  Some  shock  reflection  occurs  as  the  shock  front  encounters 
the  denser  medium  of  the  tourmaline.  This  also  sets  up  oscillations  or  multiple 
reflections  within  the  gauge  discs.  The  Bernoulli  flow  around  the  gauge  causes  the 
initial  pressure  to  decrease.  However,  this  has  been  shown  to  be  a negligible  effect 
for  shock  waves  in  water  up  to  30  Kpsi.  Diffraction  effects  and  internal  gauge 
reflections  are  considered  in  the  next  section.  Gauges  are  used  with  the  disc  edge-on 
to  minimize  and  smooth  distortions.  Face  on  gauges  disturb  the  flow,  result  in  internal 
multiple  modes  of  vibration,  and  consequently  may  not  actually  shorten  the  rise  time. 
The  rise  is  more  ragged  and  may  actually  be  longer  for  gauges  of  four  discs  face-on. 
Testing  experience  has  confirmed  these  findings. 
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3.  GAUGE  RESPONSE 


3.1  The  Shock  Wave 

The  free  field  underwater  shock  wave  is  generally  described  by  the  simple 
exponential  equation: 


p = p^e-^®  t>0  (p  = 0,t<0) 


where,  P is  the  shock  wave  peak  pressure  at  any  time,  Pm  is  the  peak  shock 
pressure:  t is  time  from  the  arrival  of  the  shock  front: 

and  © is  the  shock  wave  time  constant  defines  as  the  time  required  for  the  shock 
to  decay  to  P^^/  e,  where  e = 2.71 83. 

This  equation  is  a fair  approximation  out  to  one  theta,  but  for  longer  times  the  real 
shock  wave  decay  does  not  continue  to  decrease  as  rapidly  as  described  by  this 
equation.  This  is  not  important  for  the  upper  frequency  response  of  the  recordings 
system  or  to  the  gauge  size  selection.  It  does  have  to  be  considered  in  selecting  the 
lower  frequency  limit  of  a recording  system,  setting  the  system’s  input  time  constant 
adequately  high,  and  in  determining  the  error  introduced  by  the  electronics  system  to 
measurements  of  energy  and  impulse. 

3.2  Gauge  Selection 

The  selection  of  a gauge  for  underwater  pressure  measurements  generally 
becomes  a compromise  in  choosing  one  that  is  physically  small,  and  thus  has  a short 
transit  time  relative  to  the  shock  wave  duration  or  time  constant  0,  and  yet  large 
enough  to  have  an  adequate  gauge  constant  (KA)  to  provide  a recordable  signal. 
Other  factors  also  become  important  in  this  selection:  the  recording  system's 
sensitivity,  noise  levels,  and  frequency  response;  gauge  cable  length;  availability  of 
shock  resistant  preamplifiers  near  the  gauge;  the  physical  geometry,  setup  and 
rigging,  charge  size,  range,  cable  signal,  desired  recording  time,  boundary  and  rig 
reflections,  i.e.,  the  total  experiment.  However,  there  is  a criterion  or  gauge  selection 
that  can  be  divorced  from  the  above  factors.  This  criterion  is  based  upon  the 
acceptable  accuracy  in  the  measurement  made  by  the  finite  gauge  immersed  in  the 
underwater  shock  wave  field. 
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3.3  Selection  Criterion 


The  finite  gauge  size  affects  the  accuracy  of  measurement  of  the  shock  wave 
peak  pressure  to  the  greatest  degree.  The  peak  pressure  affects  the  measurement  of 
© directly  (Pm/e)-  Impulse  is  the  integrated  area  under  the  curve  to  an  agreed  upon 

point  (typically  5 or  10  ©),  and  energy  flux  density  is  the  integrated  area  under  the  p^ 
curve  to  the  same  points  (with  the  inverse  of  the  acoustic  impedance  as  a multiplier  out 
in  front  of  the  integral  or  1/pc). 

Both  impulse  and  energy  are  rather  insensitive  to  small  errors  in  the  shock  wave  peak 
pressure.  Therefore,  since  peak  pressure  is  the  most  sensitive  and  is  an  important 
parameter  for  explosives  comparison  and  effects,  it  is  selected  as  the  measurement 
criterion  in  gauge  size  selection.  The  peak  pressure  recorded  is  the  apparent  peak 
Papp.  The  actual  or  real  peak  pressure  is  P^.  The  response  ration,  Rp,  is  defined  as 

the  ration  of  the  Papp  Pm*  This  ratio  is  dependent  upon  the  duration  0 of  the 
exponential  shock  wave  and  the  transit  time  of  the  gauge,  t^: 

Rp  = f(0/tQ)  = Papp/^m  ^ 


3.4  Geometric  Response 

The  thin  circular  discs  that  make  up  the  gauge  have  a geometric  transit  time,  the 
time  for  the  shock  wave  to  cross  the  bare  disc  when  facing  edge-on  to  the  shock  front 
disregarding  ail  flow,  interference  or  distortion  effects.  A computer  program  was 
written  to  get  the  exponential  plane  wave  response  of  a circular  disc  edge-to-edge 
across  the  disc.  Separate  runs  were  done  with  0/tQ  ranging  from  1/16  to  2400  in 
small  Increments.  The  normalized  or  response  ratio  Rp  values  were  tabulated  and 
selected  values  were  plotted  in  Fig.  3.  The  rise  with  the  value  of  2400  is  essentially 
the  step  response. 

3.5  Actual  Measured  Response 

The  geometric  consideration  is  certainly  a simplification  of  the  real  gauge.  In  fact, 
the  gauge  covering/coating,  the  oil  boot,  must  be  considered  as  part  of  the  actual 
gauge.  No  covering/coating  has  been  found  that  Is  transparent  to  the  shock  wave.  It 
would  have  to  have  the  same  acoustic  impedance  and  properties  as  the  water,  and  yet 
afford  protection  and  insulation.  The  diffraction  of  the  shock  front  around  the  real 
gauge  will  also  tend  to  make  it  appear  larger  or  its  transit  time  longer.  The  strain  and 
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signals  induced  in  the  tabs  for  connection  to  the  cable  will  also  make  the  gauge  look 
longer  and  lengthen  the  transit  time,  these  strains  may  also  cause  transverse  modes 
of  vibration  or  oscillation. 

If  the  gauge  coating  has  little  deleterious  affect,  then  the  shape  of  the  rise  of  the 
real  gauge  should  be  controlled,  or  effectively  be  the  rise  of  the  bare  crystals 
themselves.  Consider  the  real  gauge  to  be  a black  box  with  a transfer  function 
equivalent  to  the  bare  gauge.  If  the  actual  signal  recorded  is  the  same  shape  as  the 
signal  predicted  by  computation  except  for  the  duration,  then  the  coating  has  done 
nothing  but  make  the  real  gauge  look  like  a larger  diameter  bare  gauge. 

High  frequency  response  recordings  were  made  with  an  analog  oscilloscope 
system  using  fast  sweeps  for  gauges  of  size  1 /4-inch.  The  shock  peaks  on  the  film 
records  for  the  1/4  inch  gauges  are  about  4 cm  in  height.  These  were  digitized  and 
plotted  to  an  expanded  scale  five  times  larger  (20  cm).  The  1 /2-inch  gauge  records 
were  recorded  on  a a Nicolet  digital  scope  and  plotted  to  the  same  scale.  Some 
records  have  been  obtained  for  1 /8-inch  gauges  from  large  shaped  charges  and  tiny 
cylinders.  The  few  records  appear  to  have  the  same  shape  as  those  for  the  1 /4-inch 
and  1 /2-inch  gauges,  as  they  should.  However,  these  have  not  been  plotted  to  this 
expanded  scale. 

The  results  of  these  plots  will  be  discussed  in  detail  in  the  next  section.  It  is 
sufficient  to  say  that  the  plotted  real  gauge  rise  shapes  were  excellent  matches  to  the 
computed  geometric.  It  should  also  be  noted  that  the  gauge  constant  of  the  oil-booted 
gauge  was  the  same  as  for  the  bare  gauge,  and  that  the  standard  deviation  of  a set  of 
measurements  was  reduced  with  use  of  the  oil-booted  gauge.  Separate  pressure  pot 
calibrations  of  the  calibration  constant  KA  for  the  booted  gauges  and  then  their  bare 
crystal  elements  were  identical.  The  rise  responses  of  an  oil-booted  gauge  and  a 
wax-covered  gauge  subjected  to  an  underwater  shock  are  compared  in  Fig.  4.  Notice 
that  the  wax  not  only  slows  the  rise  time  but  changes  the  wave-shape  by  rounding  and 
smoothing. 

3.6  Response  Comparisons 

The  measured  transit  time  of  the  oil-booted  gauge  is  about  half  again  as  long  as 
the  computed  geometric  transit  time  of  the  bare  gauge  element  alone,  i.e.  the  time  for 
the  shock  wave  to  traverse  the  diameter  of  the  bare  gauge  element:  tQ  (measured)  = 
1.5 10  (geometric).  Therefore,  in  order  to  compare  recorded  gauge  response  with  the 
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computed  geometric  response,  the  recorded  responses  were  reduced  by  their 
respective  transit  times,  i.e.  an  actual  gauge  rise  was  reduced  by  its  measured  transit 
time. 


The  measured  © was  divided  by  the  measured  tp  to  get  the  values  of  this  ratio. 
The  computer  runs  of  geometric  or  computer  responses  were  identified  by  this  same 
ratio. 


The  comparisons  of  reduced  actual  response  and  geometric  response  were 

plotted  for  matching  values  of  0/Tp.  The  vertical  scale  for  pressure  response  ratio  was 

also  normalized,  thus,  the  ordinate  Rp 

^norm 

Many  fast  rise  recordings  have  been  made  of  various  size  explosions  and 
compositions.  Typical  charge  weights  were:  1 oz,  1 lb,  5 lb,  10  lb,  15  lb,  50  lb,  and 
1300  lb.  Compositions  were:  C-4,  Pentolite,  TNT,  HBX-1,  PBXN-103,  and  Minol.  The 
ranges  of  0/tp  for  these  were  5-90.  A limited  number  of  fast  rise  recordings  are 
shown  in  Fig.'s  4 - 6,  for  1 /4-inch  gauges.  The  expanded  and  reduced  plots  for  a 1/4- 
inch,  and  a 1/2-inch  gauge  digitally  recorded,  are  shown  in  Fig.'s  7 - 8.  Note  how 
closely  the  actual  rises  for  the 
real  gauges  match  their  appropriate  geometric  curves. 

This  comparison  of  the  actual  rise  shape  to  that  predicted  for  the  bare  crystals  by 
the  geometric  response  proves  to  be  an  excellent  technique  for  evaluating  gauge 
coatings  and  various  other  configurations.  This  technique  along  with  the  comparison 
of  calibration  constants  (KA)  for  the  bare  crystals  and  for  the  covered  gauge  defines 
the  characteristics  of  the  completed  gauge. 

The  small  perturbations  or  oscillations  on  the  plotted  rises  of  Fig.'s  7 - 8,  are  seen 
in  the  actual  oscilloscope  records.  Fig's  4-6.  These  are  the  internal  reflections  of  the 
shock  front  within  the  tourmaline  crystals.  The  shock  speed  within  tourmaline  is  about 
four  times  faster  than  in  water.  Note  that  there  are  inflection  points  at  about  the  1/4 
points.  The  reflections  in  the  crystal  would  be  expected  to  be  damped  oscillations  with 
the  same  S-shaped  rise  and  fall. 

Note  also  the  damped  oscillations  after  the  peak.  These  are  at  lower  frequencies 
that  correspond  well  with  the  gauge  size  and  the  shock  speed  in  the  water.  They  are 
probably  mechanically  induced  stresses  or  reflections  off  gauge  tabs  or  connections. 
Further  investigations  with  modified  gauge  constructions  would  shed  more  light  on 
these  conjectures.  Whatever  the  source,  the  oscillations  are  part  of  the  real  response 
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of  the  gauge;  the  frequency  of  oscillations  is  directly  related  to  the  crystal  size  and  the 
shock  speed  in  water. 


4.  THE  NAVSWC  GAUGE  IN  USE 

The  NAVSWC  oil-booted  tourmaline  gauge  has  been  in  use  for  over  20  years  on 
hundreds  of  underwater  explosion  tests  of  various  types  and  weights  of  charges  and 
weapons.  The  gauges  have  been  in  close  proximity  to,  mounted  on  and  within 
models,  ships,  and  submarines.  Tests  have  been  conducted  in  fresh  and  salt  water 
with  a myriad  of  thermoclines  and  sound  velocity  profiles.  Tests  have  been  conducted 
in  small  chambers,  ponds,  quarries,  rivers,  bays  and  in  the  ocean  at  a multitude  of 
locations  and  environmental  conditions.  The  gauges  have  been  used  by  the  NATO 
nations  for  a means  of  data  comparisons  and  by  over  a dozen  other  countries,  as  well 
as  by  many  national  laboratories  and  companies. 

The  gauge  characteristics  have  been  well  defined  and  confidence  in  the  results 
have  been  demonstrated  repeatedly.  In  fact,  most  of  the  current  underwater  explosion 
theory  and  the  Navy's  archives  have  been  based  upon  the  analyses  of  data  measured 
with  these  gauges.  Before  one  can  believe  and  utilize  data  recorded  under  highly 
refracted  environmental  conditions  or  highly  reflected  conditions,  a gauge’s  response 
must  be  well  documented  for  the  more  standard  free  field  tests;  these  would  be  tests 
where  the  exponential  nature  of  the  shock  wave  and  bubble  pulse  are  not  interfered 
with  by  reflections  from  the  surface  and  the  bottom,  or  from  ships  or  targets  near-by,  or 
from  other  boundaries. 


4.1  Explosive  Comparison  Fresh  Water  Tests 

Free  field  (unbounded)  measurements  are  the  more  conventional  measurements 
used  to  characterize  different  types  of  explosives  and  to  compare  new  ones  with 
standards  and  others.  Generally  four  gauges  at  each  recording  station  spaced 
logarithmically  at  four  ranges  are  required  to  get  the  statistics  and  accuracies 
necessary  for  such  comparison  tests.  Some  focus  has  been  placed  on  increasing  the 
energy  in  the  bubble  pulses  in  the  past  few  years,  requiring  integration  under  the 
bubble  pulses.  Thus,  high  gain  with  good  low  frequency  response  (.01-  0.1  HZ)  must 
be  achieved  along  with  confidence  in  the  gauge's  performance,  to  faithfully  record 
some  of  the  "odd"  and  unexpected  wave-shapes.  Fig.  9 shows  the  total  record  with  the 
shock  wave  and  bubbles  compressed  in  time,  and  then  the  shock  and  each  bubble 
pulse  expanded  in  the  following  views. 
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4.2  Model  Tests 


Model  tests  typically  require  that  the  gauges  be  mounted  near  or  on  the  model 
and  sometimes  inside  flooded  compartments  that  are  air-backed.  The  reflections  and 
resulting  pressure  waves  cannot  be  predicted  and  can  be  almost  anything  (Fig.  10). 
The  point  is,  that  without  confidence  and  experience  in  a gauge's  performance  under 
more  controlled  testing,  one  cannot  differentiate  gauge  characteristics  or  problems 
from  the  "real"  phenomena.  This  could  result  In  erroneous  interpretation  and  damage 
theories.  A recent  example  emphasized  the  importance  of  this,  when  several  finite 
element  codes  predicted  "after-the-fact"  radically  different  results  with  their  damage 
explanations  without  regard  to  the  data.  The  measurements  were  then  the  only 
common  link  that  at  least  had  to  be  matched  by  the  code  before  further  extrapolations 
of  "what  happened"  could  be  made  or  believed  by  a code. 
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5.  SUMMARY  AND  CONCLUSIONS 


The  NAVSWC  oii-booted  gauge: 

• has  been  well  characterized  and  has  demonstrated  repeatable  and  consistent 
response  under  all  kinds  of  conditions  and  types  of  tests; 

• has  been  the  standard  gauge  for  recording  all  of  the  Navy's  archival  data  and 
similitude  parameters: 

• has  been  made  available  at  cost  as  a standard  for  comparison  both 
internationally  and  nationally; 

• is  rugged,  reliable,  easily  fabricated,  calibrated  and  refurbished; 

• is  linear,  bulk  or  hydrostatically  sensitive  and  is  omnidirectional; 

• effectively  does  not  change  the  amplitude  response  of  the  bare 
crystals; 

• appears  to  be  half  again  as  large  as  the  diameter  of  the  crystal,  extending  the 
response  time  to  about  1 .5  times  the  geometric  rise; 

• Is  easily  impedance  matched  to  instrumentation  systems,  the  high  frequency 
response  determined  by  the  transit  time  of  the  shock  front  crossing,  and  the  low 
frequency  or  decay  response  determined  by  the  input  or  over-all  time  constant  of 
the  system; 

• can  simultaneously  cover  extremely  wide  dynamic  pressure  ranges  simply  by 
paralleling  several  channels  of  instrumentation  at  various  gains  across  it; 

• can  be  used  with  an  in-line  preamplifiers  close-in  to  increase  the  S/N  ratio, 
with  the  resulting  limitation  In  dynamic  pressure  range  caused  by  the 
preamplifier; 

• and  remains  as  the  transducer  of  choice  for  underwater  explosion  pressure 
measurements. 
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FIG.  1 NAVSWC  TOURMALINE  GAUGE  ASSEMBLY 
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FIG.  2 NAVSWC  TOURMALINE  GAUGE  DIAGRAM 
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FIG.  3 GEOMETRIC  GAUGE  RISE  RESPONSES 
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5 ^lsec/cm 

FIG.  4 GAUGE  RISE  COMPARISON 
OIL  BOOTED  VS.  WAX 


5 fisec/cm 


FIG.  5 GAUGE  RISES  - OIL  BOOTED 
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lOjLt  sec/cm 


FIG.  6 A GAUGE  RISES 


200  /isec/cm 

FIG.  6B  SHOCKWAVE  RECORD  FOR  GAUGE  RISES 
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^pNORM 


FIG.  7 ACTUAL  vs.  GEOMETRIC  RISE  1/4-INCH  GAUGE 
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RpNORM 


FIG.  8 ACTUAL  vs.  GEOMETRIC  RISE  1 /2-INCH  GAUGE 
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PRESSURE  vs  TIME 


PRESSURE  vs  TIME 


PRESSURE  vs  TIME 


FIG.  10  MODEL  TEST- INTERNAL  GAUGES 
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ARMY  DYNAMIC  PRESSURE  MEASUREMENTS 
By  John  M.  Ball 
U.S.  Army  IMDE  Activity 
Redstone  Arsenal,  AL 

ABSTRACT:  This  paper  briefly  describes  dynamic  pressure  measurements  in  the  Army, 
tests  performed,  ranges  of  interest,  and  transducers.  Army  calibration  program  and 
traceability  requirements  are  outlined  along  with  measurement  concerns.  Several 
actual  calibration  systems  are  described.  The  need  for  a national  calibration 
service  and  measurement  program  in  this  parameter  is  emphasized. 


I.  INTRODUCTION. 

A national  survey  of  Army  dynamic 
pressure  measurement  activities  was 
conducted  in  1986  [1].  This  census, 
while  not  all  inclusive,  identified 
quantities  of  transducers,  ranges  of 
interest,  and  locations  at  which 
dynamic  pressure  measurements  were 
conducted.  Test  ranges,  proving 
grounds,  and  development  centers 
constituted  the  principal  users  of 
dynamic  pressure  transducers.  While 
ranges  of  interest  varied  with  the 
application,  measurements  from 
thousands  of  atmospheres  to  tens  of 
thousands  of  atmospheres  were  then 
common  and  are  currently  being 
performed  on  a routine  basis  in  a 
number  of  locations.  Higher  range 
measurements  usually  involve  combustion 
generated  pressures:  chamber  pressures, 
blast  pressures,  and  so  on.  Lower 
range  measurements  are  applied  to  such 
things  as  engine  pressures.  Products 
from  most  major  manufacturers  of 
dynamic  pressure  transducers  and 
unique.  Army  designed  and  manufactured 
transducers  are  represented.  Hundreds 
of  these  transducers  are  formally 
supported  through  the  Army  calibration 
process.  The  typical  transducer  had  a 
range  of  about  6800  Atm  (100,000  PSI), 
with  an  expected  error  of  1%  to  2%. 

A study  conducted  in  1990 
demonstrated  that  certain  of  these 
dynamic  pressure  machines  generated 
pressures  which  agreed  with  precision 
static  measurements  to  about  one 
percent. [5]  It  was  not  possible, 
however,  to  directly  validate  the 
performance  of  these  devices  in  their 
dynamic  mode. 


II.  ARMY  CALIBRATION  PROGRAM 

By  regulation  [2],  the  U.  S.  Army 
operates  a calibration  program 
designed  to  ensure  the  accuracy  of 
measurements  in  all  parameters  and  to 
guarantee  the  correct  performance  of 
Army  test,  measurement,  and  diagnostic 
equipment.  The  test  ranges,  proving 
grounds,  RD&E  Centers,  and 
laboratories  which  perform  dynamic 
pressure  measurements  are  supported  by 
calibration  and  repair  services  from 
mobile  calibration  vans  and 
sophisticated,  fully  equipped 
measurement  laboratories.  These 
activities  are  controlled,  operated, 
and  supported  through  a hierarchy  of 
calibration  laboratories,  including 
the  Army  Primary  Standards  Laboratory, 
which  functions  as  the  "bureau  of 
standards"  for  the  Army. 

Program  management,  engineering, 
and  logistics  for  these  calibration 
functions  and  for  general  purpose  Army 
test  instrumentation  and  automated 
test  equipment  are  centralized  in  the 
U.  S.  Army  TMDE  Activity  headquartered 
at  Redstone  Arsenal,  AL. 

III.  TRACEABILITY 

Calibration  generally  requires  a 
measurement  system  with  known  accuracy 
and  a test  unit  which  is  certified 
against  this  system.  Regulations 
govern  accuracy  relationships  between 
measurement  systems  and  items  to  be 
calibrated  [3].  In  many  cases,  a 
transducer  which  has  been  calibrated 
is  used  as  an  "artifact  standard"  to 
transfer  accuracy  to  other  items. 
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Periodic  calibration  detects  changes  in 
performance  and  permits  correction  or 
adjustment.  The  accuracy  of  a 
calibration  is  dependent  upon  a 
reference  which  may  be  another  piece  of 
hardware,  a physical  phenomena,  or  a 
process.  The  National  Institute  of 
Standards  and  Technology  (NIST),  is  the 
organization  in  this  country  to  which 
all  measurements  except  time  are 
legally  referenced.  "Traceability"  is 
the  process  by  which  the  relationship 
of  field  measurements  to  NIST  is 
documented. 

Dynamic  pressure  measurements, 
however,  are  not  directly  traceable  to 
NIST.  This  is  the  most  significant 
issue  in  this  parameter  today.  Neither 
does  the  NIST  provide  dynamic  pressure 
measurement  services.  Without 
traceability,  the  accuracy  of  such 
measurements  cannot  be  directly 
validated  at  any  level  in  the 
calibration  hierarchy,  calibration 
support  is  reduced  to  the  certification 
of  repeatability  at  some  inexactly 
known  value,  and  the  Army's  testing 
program  is  seriously  compromised. 

Army  users  and  calibration 
operations  have  applied  several 
different  calibration  approaches  in  an 
attempt  to  provide  measurement  services 
and  some  degree  of  traceability. 
Several  will  be  described  briefly 
later  in  this  paper. 

IV.  THE  PHENOMENA 

Dynamic  pressure  measurements  are 
almost  always  measurements  of  rapid 
pulses.  In  many  applications,  the 
pressure  pulse  is  accompanied  by  a 
dynamic  temperature  pulse,  the  fluid 
being  measured  is  a combusting  gas,  and 
powerful  mechanical  shocks  and 
vibrations  exist  in  the  test  item. 
Consider  as  an  illustration  a 
piezoelectric  transducer  mounted  to 
measure  chamber  pressure  in  a cannon, 
the  nature  of  the  fluid,  the 
simultaneous  mechanical  shock  pulses. 


vibrations,  environmental  effects, 
mounting  conditions,  and  the 
construction  of  transducers  commonly 
applied.  Transducers  in  such  an 
installation  experience  mechanical 
stresses,  electromagnetic  pulses,  and 
temperature  shocks  in  addition  to  the 
parameter  of  interest.  To  physically 
protect  the  transducer  from  the 
process,  grease,  modeling  clay,  or 
some  other  material  is  often  packed 
into  the  mounting  hole. 

The  user  is  generally  interested 
not  only  in  the  amplitude  of  the 
pressure  pulse  but  in  capturing  the 
pulse  for  analysis  (Figure  1). 

V.  TRANSDUCER 

Because  the  process  being  measured 
is  such  a strain  upon  the  transducer, 
it  has  become  a common  practice  to 
"calibrate"  these  devices  prior  to  and 
just  after  a test  to  detect  drastic 
performance  changes  which  may  have 
occurred  during  the  measurement 
process. 

Dynamic  pressure  transducers 
exhibit  hysteresis,  nonlinearity, 
nonrepeatability,  zero  shift, 
resonances,  and  sudden,  catastrophic 
failure.  Their  output  is  affected  by 
strain  on  the  cables,  thermal 
transients  and  shocks,  mechanical 
vibrations,  electrical  effects  (ground 
loops,  lead  resistances,  etc.), 
ambient  conditions,  photo  flash,  over 
pressure,  and  warm  up  time. 

Transducers  are  often  specified  in 
terms  of  the  root  mean  square  (RMS)  of 
linearity  to  a "best  straight  line 
fit"  of  pressure  data  to  transducer 
response  based  upon  static  pressure 
values,  repeatability,  and  hysteresis. 
Frequency  response  is  generally 
indicated  by  specifying  the  resonant 
frequency  of  the  device.  Other 
sources  of  potential  error  (see 
Table),  are  well  documented  by 
manufacturers  and  researchers  [4]. 
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Figures  2 and  3 illustrate  typical 
temperature  related  effects  for  one 
type  transducer.  Even  though  these 
effects  are  relatively  large  compared 
to  the  accuracy  of  the  measurement 
desired,  the  effect  of  transient 
temperatures  can  be  much  greater  and  is 
not  so  predictable.  Figure  4 is 
based  upon  illustrations  from  a study 
of  the  performance  of  miniature 
pressure  transducers.  When  used 
without  a specially  constructed  thermal 
damper,  40%  performance  shifts  were 
observed.  The  same  study  found  the 
effects  of  mechanical  vibrations  to  be 
as  much  as  25%  of  full  scale  transducer 
output  [6]. 

The  grease  or  modeling  clay  which 
protect  the  sensor,  the  orifice 
snubber,  protective  shield,  diaphragm, 
liquid  cell,  and  mechanical 
construction  of  the  device  itself  all 
act  as  low  pass  filtering  elements, 
guaranteeing  that  the  signal  from  the 
transducer  is  not  identical  to  the 
actual  process  pulse  of  interest.  An 
accurate  characterization  of  the 
relationship  between  the  process  and 
the  senor  response  is  required  to 
reproduce  process  information  from 
transducer  signal.  Not  only  pressure 
information  but  also  frequency  response 
and  repeatability  data  are  essential  to 
the  utilization  of  dynamic  test 
results.  The  validity  of  such  a 
characterization  is  a direct  function 
of  the  calibration  process. 

VI . MEASUREMENT  SUPPORT 

Whether  NIST  support  exists  or  not. 
Army  dynamic  pressure  transducers  must 
be  certified  and  their  performance 
verified.  A variety  of  Army  systems 
and  commercial  approaches  to  providing 
this  measurement  support  are  in  common 
use.  Because  direct  NIST  traceability 
does  not  exist,  each  approach  depends 
upon  indirect  means  and  calculations 
for  estimating  system  accuracy.  Of 
course,  stated  accuracies  can  often  be 
neither  independently  verified  nor 
disproved. 


Direct  calibration  traceability  to 
NIST  through  an  artifact  is  not 
necessary  for  a measurement  to  be 
traceable  to  NIST.  If  a particular 
measurement  depends  in  a well 
understood  manner  upon  other 
parameters  which  are  directly  NIST 
traceable,  it  is  theoretically 
possible  to  determine  the  accuracy  of 
the  measurement  of  interest  through  an 
error  analysis  calculation.  Such  is 
the  justification  for  the  systems 
below  being  identified  as  calibrators 
rather  than  testers. 

VII.  SAMPLE  SYSTEMS 

Several  actual,  Army  calibration 
systems  are  described  below.  None  of 
these  devices  claims  to  adequately 
simulate  actual  usage  conditions. 

Shock  tubes  have  been  in  use  for 
years  to  produce  known  pressure  pulses 
in  gas  (Figure  5).  The  chamber  on  the 
left  is  pressurized  behind  the 
diaphragm,  which  is  burst  when  a 
predetermined  pressure  level  has  been 
reached.  A pressure  wave  travels  down 
the  tube,  stimulating  the  transducer. 
Calculations  from  gas  theory  are  a 
possible  basis;  however,  pressure 
data  is  generally  obtained  through 
comparison  of  the  transducer  under 
test  to  another  transducer,  a 
calibration  standard.  Shock  tubes 
yield  good  frequency  data.  Weaknesses 
in  the  method  include  difficulty  in 
reproducing  the  same  pulse  at  the 
transducer,  very  slow  repetition  rate, 
and  a generally  low  maximum  pressure 
(hundreds  rather  than  thousands  of 
atmospheres).  Also,  at  high  pressures, 
such  devices  are  somewhat  dangerous. 

In  contrast  to  the  situation  in 
dynamic  pressure,  static  pressure 
measurement  services,  calibration,  and 
NIST  traceability  are  readily 
available.  The  accuracy  of  such 
measurements  is  parts  per  million  up 
to  about  7,000  Atm  for  liquid 
pressure.  The  device  illustrated  in 
Figure  6 takes  advantage  of  the 
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extreme  accuracy  and  repeatability  of 
the  dead  weight  piston  gage  to  produce 
calibrations.  A precise  pressure  head 
is  generated  upon  the  transducer  under 
test  and  then  the  valve  is  thrown, 
rapidly  releasing  the  pressure. 
Transducer  response  is  measured  for  a 
series  of  pressure  levels.  This  device 
is  ideally  suited  for  performing 
repeated  tests  at  very  well  defined 
pressure  levels,  for  example,  before 
and  after  test  range  shots.  Concerns 
include  the  use  of  constant  temperature 
oil  to  simulate  the  exploding  gas 
actually  present  under  use  conditions, 
and  the  assumption  that  the  transducer 
under  test  responds  in  the  same  manner 
to  a static  pressure  suddenly  bled  off 
as  to  a positive  step  pressure.  Such 
a system  has  excellent  repeatability. 

Another  design,  also  using  a piston 
in  cylinder,  generates  an  impulse  in 
the  fluid  by  dropping  an  anvil  on  the 
piston  (Figure  7).  Force  applied  is 
calculated  from  anvil  weight  and  height 
dropped,  pressure  from  force  divided  by 
cross  sectional  area.  While  the 
pressure  generated  is  less  well  known 
than  in  the  previous  system,  a true 
pulse  is  produced. 

A third  liquid  based  design  is 
built  around  a very  fast  acting  valve 
(Figure  8.)  The  chamber  is  pressurized 
and  the  static  pressure  measured  by  a 
pressure  sensor  which  need  not  be  fast, 
but  should  be  accurate.  The  valve 
fires,  exposing  the  transducer  to  a 
sudden,  positive  pressure  step.  The 
pressure  gage  records  the  new,  slightly 
reduced  pressure  to  which  the 
transducer  was  exposed.  It  should  be 
noted  that  the  pressure  gage  is  never 
exposed  to  significant  stresses.  This 
design  shares  the  weaknesses  of  the 
other  liquid  calibrators,  also  assumes 
symmetry  of  response,  and  has  the 
additional  problem  of  the  fast  acting 
valve  which  experiences  terrific 
stresses  and  must  be  rebuilt  after 
every  few  dozen  uses.  It  does, 
however,  repeatably  generate  very  high 
pressure  calibration  steps  and 


correlates  well  with  other  measurement 
methods. 

A final  example  is  the  pneumatic 
step  function  generator  illustrated  in 
Figure  9.  This  device  uses  a uniquely 
actuated,  very  fast  poppet  valve  to 
generate  step  pressure  pulse  on  a 
static  pressure.  It  has  proven  useful 
for  generating  accurately  known, 
repeatable  performance  data  with 
either  positive  or  negative  steps. 
This  calibrator,  however,  is  limited 
to  about  2,000  PSI. 

VII.  SUMMARY  AND  CONCLUSIONS 

Requirements  currently  exist  in 
the  Army  for  dynamic  pressure 
measurements  and  calibration  support. 
Future  requirements  will  probably  far 
exceed  those  of  today  [6]. 
Calibrators  which  are  in  use  do  not 
very  well  simulate  actual  use 
conditions  and  accuracies  stated  are 
based  upon  estimations,  comparisons, 
repeatability,  and  experience,  rather 
than  direct  traceability  to  national 
standards.  On  the  other  hand, 
commercial  and  custom  made  calibrators 
probably  perform  sufficiently  well  to 
be  practical  in  most  applications,  and 
could  be  made  much  more  accurate  if 
transfer  standards  existed  which  could 
certify  and  characterize  their 
performance. 

A national  measurement  service 
with  a solid,  scientific  basis  is 
sorely  needed  and  an  associated 
calibration  service  for  transducers. 
The  national  measurement  system  should 
measure  dynamic  pressure  and  dynamic 
temperature  in  gas.  The  system  should 
be  capable  of  generating  a variety  of 
pulse  and  step  widths.  There  are 
numeous  legal  and  practical, 
operational  requirements  for  this 
service.  Transducers  which  currently 
exist  might  well  prove  adequate  as 
transfer  artifact  standards  which 
could  be  used  to  calibrate  many  of  the 
calibration  systems  now  in  use, 
providing  real  NIST  traceability. 
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Measurement  accuracy  could  well  be 
improved  an  order  of  magnitude  or  more 
without  significant  changes  to  the 
current  hardware  base  in  the  Army  with 
true  traceability. 

Development  of  a service  for 
calibrating  dynamic  pressure 
transducers  will  not  only  improve 
Armysupport  in  this  parameter,  but 
provide  a significant  national  resource 
for  military  and  civilian  industries. 
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Figure  5 
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RESPONSE  OF  MINIATURE  DYNAMIC 
PRESSURE  TRANSDUCERS 


Figure  4 
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Input  Resistance 

Output  Resistance 

Zero  Error  Due  to  Temperature 
-18  ° to  65  ®C 
-54  ° to  74  *C 
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-18  ° to  65  ®C 
-54  ° to  74  °C 

Diaphragm  Resonant  Frequency 
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50  psia 
100  psia 
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Flat  Frequency  Response  of  Trans- 
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Warmup  Time  (1%  Accuracy) 
Acceleration  Sensitivity 
Zero  Shift  with  Mounting  Torque 
Nonlinearity  at  2X  Range 
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ABSTRACT 

The  National  Aeronautics  and  Space  Administration  (NASA) 
performs  a variety  of  research  requiring  the  measurement  of 
transient  pressures.  To  verify  the  accuracy  of  these 
measurements,  NASA  presently  must  use  calibration  techniques 
that  were  developed  by  pressure  transducer  manufacturers,  but 
these  techniques  are  unsatisfactory  and  fail  to  meet  all  NASA 
requirements. 

NASA  and  the  National  Institute  of  Standards  and  Technology 
(NIST)  metrologists  have  collaborated  to  develop  an  alternative 
calibration  approach  that  will  resolve  problems  with  current 
calibration  techniques  and  will  meet  critical  research  requirements. 
This  paper  will  discuss  current  calibration  limitations,  NASA 
research  requirements,  and  the  proposal  for  the  alternative 
calibration  approach. 
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INTRODUCTION 


NASA  research  is  dynamic  and  challenging,  both  in  the  type  of  research 
conducted  and  in  the  type  of  measurements  made.  One  area  that  has  become 
especially  challenging  is  the  measurement  of  transient  pressures.  The  challenge 
is  to  provide  accurate  measurements  under  actual  use  conditions. 

To  ensure  accuracy,  NASA  requires  that  all  instruments  used  to  gather  data  be 
calibrated  and  that  all  calibrations  be  traceable  to  national  standards.  NIST  does 
not  provide  a transient  pressure  standard,  so  NASA  presently  must  use 
calibration  techniques  that  were  developed  by  manufacturers. 

One  manufacturer's  technique  requires  two  distinct  devices:  a vacuum-driven 
shock  tube  and  a hydraulic  drop  tube.  The  shock  tube  measures  the  time- 
related  responses,  such  as  rise  time.  The  drop  tube  establishes  sensitivity. 
Between  them,  the  pair  of  devices  can  generate  pressures  with  amplitudes  in 
excess  of  100  MPa  {15,000  psi)  and  rise  times  shorter  than  one  microsecond. 

Another  calibration  technique,  based  on  the  Aronson  calibrator,  measures  both 
time-related  responses  and  sensitivity.  This  calibrator  generates  pressure  pulses 
up  to  7 MPa  (1000  psi)  with  rise  times  of  about  35  microseconds. 

Neither  calibration  technique  is  satisfying.  In  the  first  technique,  the  shock 
tube's  vacuum-to-atmospheric  pressure  step  is  quite  different  from  the 
atmospheric-to-pressurized  step  usually  encountered  in  actual  testing. 

Calibration  with  the  drop  tube  is  tedious  and  time-consuming,  and  its  resolution 
is  insufficient  for  the  lower  pressure  ranges.  Also,  the  properties  of  the 
hydraulic  fluid  in  the  drop  tube  are  quite  different  from  those  of  the  gas 
encountered  during  actual  measuring.  The  technique  based  on  the  Aronson 
calibrator  fails  to  produce  high  enough  pressures  or  fast  enough  rise  times  to 
meet  all  NASA  requirements.  Finally,  neither  of  these  techniques  can  measure 
other  important  factors,  such  as  acceleration  sensitivity  or  flash  sensitivity,  nor 
can  they  be  used  as  research  tools  in  studies  of  phenomena  such  as  supersonic 
shock  waves. 

An  ideal  calibration  system  would  combine  the  short  rise  times  of  a shock  tube 
and  the  variable  pressures  of  a drop  tube,  without  the  drawbacks  of  the  latter. 


REQUIREMENTS 

NASA  research  requiring  transient  pressure  measurement  falls  into  three  main 
categories:  explosion,  aerodynamic,  and  engine.  The  specific  requirements  of 
each  category  follow. 
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Explosion  Research 

This  research  investigates  blasts  and  associated  shock  phenomena,  such  as 
would  be  found  in  hydrogen  and  oxygen  mixing  explosions,  ballistic  propellant 
measurements,  or  liquid  fuel  explosions  in  air.  Explosion  research  uses  only 
piezoelectric  transducers  for  the  measurements.  Typically,  the  transducers  must 
measure  a single,  high-pressure  pulse  with  a short  rise  time.  A special  need  is 
separation  of  the  thermal  response  from  the  pressure  response  in  flash 
environments. 

Table  1 illustrates  some  general  requirements  of  explosion  research. 


Table  1 

Explosion  Research  Requirements 


Parameter 

Range 

Accuracy 

Sensitivity 

500  MPa  (70,000  psi) 

1% 

Resolution 

0.1%  Full  Scale 

Not  Applicable 

Rise  Time 

1 /vs 

5% 

Resonant  Frequency 

500  kHz 

10% 

Linearity 

Not  Applicable 

1% 

Repeatability 

Not  Applicable 

1% 

Temperature 

Coefficient 

-160  to  150  °C 

0.05%  per  °C 

Aerodynamic  Research 

This  research  investigates  the  interaction  of  air  frames  with  gas  streams,  such 
as  would  be  found  in  studies  of  shock  fronts,  boundary  layers,  and  turbulent 
zones.  Aerodynamic  research  generally  uses  strain  gage  transducers,  but 
piezoelectric  and  piezoresistive  transducers  are  becoming  more  common. 
Typically,  the  transducers  must  measure  relatively  small  pressures  that  fluctuate 
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rapidly.  At  times,  they  will  need  to  be  able  to  measure  both  transient  and  static 
responses  during  the  same  test  run.  Transducers  are  usually  exposed  to 
moderate  temperatures,  although  tests  involving  vehicle  re-entry  into  the 
atmosphere  have  very  high  temperatures. 

Table  2 shows  some  general  requirements  of  aerodynamic  research. 


Table  2 

Aerodynamic  Research  Requirements 


Parameter 

Range 

Accuracy 

Sensitivity 

700  kPa  (100  psi) 

0.1% 

Resolution 

1 % of  Reading 

Not  Applicable 

Rise  Time 

1 //s 

5% 

Resonant  Frequency 

1 MHz 

0.1% 

Linearity 

Not  Applicable 

0.1% 

Repeatability 

Not  Applicable 

0.1% 

Temperature 

Coefficient 

-170  to  1,200  °C 

0.05%  per  °C 

Engine  Research 

This  research  investigates  the  performance  of  rocket,  jet,  and  heat  engines. 
These  tests  use  strain  gage,  piezoelectric,  and  piezoresistive  transducers. 
Typically,  the  transducers  must  measure  sustained,  high  pressures  that  rapidly 
fluctuate.  They  are  usually  exposed  to  severe  environments.  Often  the 
structures  required  to  protect  the  transducers  complicate  their  response,  so 
actual  use  conditions  require  transducer  cooling  and  temperature  compensation. 
Another  special  need  is  determination  of  the  phase  lag  introduced  by  the 
transducer  in  heat  engine  research. 
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Table  3 shows  some  general  requirements  of  engine  research. 


Table  3 

Engine  Research  Requirements 


Parameter 

Range 

Accuracy 

Sensitivity 

100  MPa  (15,000  psi) 

1% 

Resolution 

0.1%  Full  Scale 

Not  Applicable 

Rise  Time 

5 fjs 

5% 

Resonant  Frequency 

225  kHz 

1% 

Linearity 

Not  Applicable 

0.5% 

Repeatability 

Not  Applicable 

0.25% 

Temperature 

Coefficient 

-50  to  300  °C 

0.02%  per  °C 

Additional  Requirements 

In  addition  to  the  requirements  listed  in  the  tables,  ringing  and  acceleration 
sensitivity  are  also  important  factors  when  measuring  the  responses  of  a 
transducer. 


PROPOSAL 

NASA  and  NIST  metrologists  collaborated  to  develop  a joint  proposal  that 
addresses  the  critical  NASA  requirements.  The  objective,  approach,  and 
benefits  of  the  proposal  follow. 
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Objective 


The  objective  is  to  develop  a calibration  system  for  characterizing  transient 
pressure  instruments. 

The  calibration  system  will  include  a gas-driven  shock  tube  capable  of  producing 
shock  pressures  up  to  100  MPa  (15,000  psi)  and  spectroscopic  instrumentation 
to  determine  the  pressure  using  inherent  properties  of  the  gas  in  the  tube.  It  will 
be  designed  so  that  it  can  be  maintained  by  calibration  laboratories  and 
compared  to  the  calibration  system  maintained  by  NIST. 


Approach 

The  collaborative  development  will  be  conducted  in  three  phases. 

The  first  phase  will  be  the  design  and  fabrication  of  the  calibration  systems. 
Both  NASA  and  NIST  will  construct  calibration  systems  for  future  cross- 
correlation. 

The  second  phase  will  be  the  testing  and  evaluation  of  the  performance  of  the 
calibration  systems  over  their  full  pressure  ranges.  Accuracy  limits  will  be 
determined. 

The  third  phase  will  be  the  cross-correlation  of  the  calibration  systems  using 
calibrated  transducers.  The  cross-correlations  will  verify  the  traceability  of  the 
calibration  process. 


Benefits 

The  most  important  benefit  is  that  NIST  will  have  established  a national  standard 
for  transient  pressure.  NIST  will  provide  traceable  calibrations  to  NASA,  which 
will  maintain  a reference  standard  for  calibration  of  NASA  instruments. 

Other  agencies  and  organizations  will  benefit  as  well.  Once  this  research  is 
complete,  other  shock  tubes  can  be  fabricated  and  calibrated  against  the  NIST 
national  standard. 

An  additional  benefit  is  that  the  calibration  system  has  the  potential  of  being 
developed  into  an  intrinsic  standard,  which  would  reduce  the  dependence  on 
NIST  for  calibration  of  pressure  transducers. 

Finally,  this  calibration  system  will  contribute  to  a better  understanding  of  the 
extremely  dynamic  phenomena  observed  in  wind  tunnels,  rocket  engines, 
explosions,  and  supersonic  shock  waves. 
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FAST  PRESSURE  MEASUREMENTS: 

THE  NEED  FOR  PVDF  CALIBRATION  STANDARDS 


Bruce  Hartenbaum 
H-Tech  Laboratories,  Inc. 
California  Institute  of  Technology 

Bradford  Sturtevant 
California  Institute  of  Technology 


Fast,  nanosecond  street  measurements  in  solids  can  be  made  by  using  a variety  of 
transducers  including  manganin,  ytterbium,  quartz,  carbon,  and  pol3rvinlyidene 
fluoride  (PVDF).  In  water,  PVDF  is  the  transducer  of  choice  because  of  its  close 
impedance  match  to  water. 

There  are  several  questions,  however,  concerning  the  accuracy  of  PVDF  pressure 
measurements  in  water  at  less  than  1 kb,  which  is  our  primary  area  of  interest: 

(1)  What  are  the  effects  of  the  strength  of  PVDF  on  charge  production 
when  the  transducer  is  placed  in  a fluid  subjected  to  uniaxial  strain?  In 
the  simplest  case  the  d33  axis  is  aligned  with  the  direction  of  strain  in 
the  fluid,  but  the  PVDF  experiences  triaxial  stress  and  the  d^^  respond 
accordingly. 

(2)  What  is  the  response  of  PVDF  to  an  oblique  plane  wave?  This  case  is 
similar  to,  but  more  complex  than  case  (1). 

(3)  What  are  the  causes  of  observed  offsets?  Is  the  offset  caused  by  the 
PVDF  itself  or  by  the  data  recording  methods? 

To  answer  these  questions  we  are  investigating  two  calibration  and  several  data 
reduction  techniques  that  should  be  of  interest  to  the  NIST: 

(1)  The  use  of  well-established  electrically  driven  flyer  plate  and  shock 
Hugoniot  technology  to  produce  a planar  shock  of  known  pressure  in  a water 
cell  for  calibration  of  PVDF  transducers. 

(2)  The  use  of  differential  interferometry  in  a water  cell  to  infer  the 
pressure  field  behind  a nonuniform  shock,  to  provide  an  alternate  measure 
of  pressure  for  evaluating  PVDF. 

(3)  Three  data  reduction  methods:  (i)  numerical  analysis  of  digital 
data,  (ii)  analog  integration,  and  (iii)  direct  recording  as  used  in  the 
acoustics  community. 

It  is  believed  that  the  NIST  could  provide  a needed  role  by  developing  standards 
for  calibration  of  fast  nanosecond  pressure  transducers  and  by  ensuring  the 
accuracy  and  the  equivalence  of  shock  Hugoniot  and  acoustic  reciprocity  methods. 
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• POLYVINYLIDENE  FLUORIDE  (PVDF)  THIN  FILM  TRANSDUCERS  FOR  FAST 
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OTHER  APPLICATIONS 


Pressure  (MPa)  Pressure  (MPa) 


LITHOTRIPTER  SHOCK  WAVE  PRESSURE  WAVEFORMS 
AT  SHOCK  WAVE  GENERATOR  EXIT  PLANE  & 

AT  SHOCK  FOCUS  IN  WATER  BATH 
(ELECTROHYDRAULIC  GENERATOR) 


104  106  108  110  112  114  116 

Time  {fjs) 

CL  Exit  Plane 


88 


WHAT'S  THE  PROBLEM? 


89 


PVDF  SHOCK  CALIBRATION  DATA  (SANDIA) 
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SHOCK  STRENGTH  (PASCAL) 
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CHARGE  EFFICIENCY  = POLARIZATION  CHARGE  / REMANENT  POLARIZATION 


CEA  UNLOADING  {&  LOADING!)  DATA 
PVDF  COMPARED  WITH  CARBON  GAGE 
DETONATOR  SOURCE  — WHAT  IS  THE  INPUT  WAVEFORM? 


CD 

L_ 

C/5 

CO  I— 
CD  CD 
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time,  M s 


H-TECH  UNLOADING  (8  LOADING!)  DATA 
SHOCK  PRESSURE  AT  LITHOTRIPTER  GENERATOR  EXIT  PLANE 
VARIATION  IN  AMPLITUDE  OF  NEGATIVE  PHASE 
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TO  COMPLICATE  MATTERS... 


NONUNIFORM  PRESSURE  FIELD  IN  LITHOTRIPTER  — 
NONUNIFORM  LENGTH  SCALE  COMPARABLE  TO  SIZE  OF  TRANSDUCER 


CALCULATED  PRESSURE  FIELD.  STRONG  INPUT  SHOCK  WITH  CONSTANT  STRENGTH. 


SHOCK  rocus.^ 


CALCULATED  PRESSURE  FIELD.  STRONG  INPUT  SHOCK  WITH  VARIABLE  STRENGTH. 
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TO  FURTHER  COMPLICATE  MATTERS... 

ACCURACY  OF  PEAK  PRESSURE  MEASUREMENT  DEPENDS  ON  SIGNAL  CONDITIONING 


STEP  RESPONSE  OF  MODEL  OF  CIRCUIT  COMMONLY  REPORTED  IN  LITERATURE 
(65  CM  CABLE,  1 MEGOHM  INPUT  IMPEDANCE)  3 50  NS/DIV 


STEP  RESPONSE  OF  H-TECH  COMPENSATED  CIRCUIT  9 50  NS/DIV 
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SHOCK  HUGONIOT  CALIBRATION  OF  PVDF  IN  WATER 


SHOCK  STATE  IN  WATER: 

PV^  = PwCy^Upy^ 

SHOCK  STATE  IN  FLYER  PLATE  AFTER  IMPACTING  WATER: 

Pf  = PfCf(Uf  - Upw) 

UNIFORM  PRESSURE  CONDITION:  P^  = Pf 


Pressure 

Kbar 


water 
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SCHEMATIC  OF  ELECTRICALLY  DRIVEN  FLYER  PLATE  SHOCK  GENERATOR 

FOR  CALIBRATING  PVDF  IN  WATER 
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water  cell 


SHOCK  HUGONIOT  CALIBRATION 
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CALIBRATION  AND  ULTRASONIC,  FREQUENCY  DOMAIN  CAL 


DYNAMIC 

STEP-PRESSURE 

CALIBRATION 


Presented  at:  NIST  (NBS)  Workshop,  Gaithersburg,  Md. 
MEASUREMENT  OF  TRANSIENT 

PRESSURE  AND  TEMPERATURE 

23-24  April  1991 


By:  J.F.  Lally,  President 

PCB  Piezotronics,  Inc. 

3425  Walden  Avenue 
Depew,  NY  14043 


ptff 

PIEZOTRONICS 


New  "ULTIMA" 

DYNAMIC  PRESSURE  CALIBRATOR 


Zi 

n 

ri 

ri 

c: 

1.1 

ILI 

JLI 

Environmental 

Chamber 


PROVIDES  A CONTINUOUSLY  ADJUSTABLE  STEP,  PULSE  OR  SINE  WAVE 
PRESSURE  OF  KNOWN  AMPLITUDE  FROM: 


0.00001  psi  toSOOOOOpsi 
DC  to  10  MHz 
Absolute  0 to  +2  000  °F 


• COMPUTER  CONTROLLED 

• AUTOMATIC,  SELF-COMPENSATING 

• COMPLETE  WITH  ALL  BELLS  & WHISTLES 

• ECONOMICALLY  PRICED  AT  $99.95 

• DELIVERY  FROM  STOCK 

• NIST  TRACEABLE  Manufactured  by 

"VISIONARY"  Calibrator  Co. 
UTOPIA,  USA 
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Obviously,  it  would  be  nice  if  the  "ULTIMA"  dynamic  calibration  device  existed. 

Even  if  it  did  exist,  I'm  not  sure  how  it  could  be  "calibrated".  If  the  "ULTIMA"  calibrator 
were  dismantled,  you  would  find  it  is  made  up  of  several  individual  special  purpose  cali- 
bration devices  - each  designed  to  dynamically  calibrate  pressure  transducers  over  a 
limited  pressure,  frequency  or  temperature  range.  Each  calibrator  would  pretty  much 
represent  the  current  state-of-the-art  of  dynamic  calibration  over  the  specific  operating 
range  of  the  device. 

In  1972,  considerable  effort  was  expended  by  the  American  Society  of  Mechanical 
Engineers  to  establish  "A  Guide  for  the  Dynamic  Calibration  of  Pressure  Transducers". 
The  guide  was  originally  published  as  American  National  Standard  ANSI  B88.1-1972; 
however,  the  current  number  is  ANSI  MC88.1-1972  (see  page  3). 

The  ANSI  calibration  guide  defines  calibration  as  follows: 

"Calibration  means  a test  during  which  known  values  of 
measurand  are  applied  to  a transducer  and  corresponding 
output  readings  are  recorded." 

The  ANSI  document  then  proceeds  to  qualify  expectations  for  lower  accuracy  and 
the  somewhat  more  difficult  nature  of  dynamic  calibration; 

"The  degree  of  accuracy  associated  with  these  dynamic  tests 
is  generally  lower,  and  the  manner  in  which  the  results  are  used 
is  generally  less  rigorous  than  in  the  conventional  and  more  easily 
controllable  field  of  static  pressure  calibration." 

Much  of  the  dynamic  pressure  calibration  equipment  used  today  has  evolved  out 
of  the  transducer  manufacturer's  necessity  to  test  and  calibrate  his  products,  and  by  the 
customers'  need  for  in-house  capability  to  test  and  calibrate  transducers  more  closely  to 
the  customers'  actual  application  requirements.  A wide  variety  of  calibration  devices  exist, 
many  of  which  are  offered  as  commercial  products.  Several  are  described  in  the  ANSI 
dynamic  calibration  guide: 

• PISTONPHONE 

• PULSE  CALIBRATORS  USING  FAST-ACTING  VALVES 

• SINUSOIDAL  PRESSURE  GENERATORS 

• IMPULSE  DROP  CALIBRATORS 

• SHOCK  TUBE 

Most  of  the  above  are  traceable  to  the  NIST  through  an  accurate  DC  reference 
gage,  comparison  transfer  standard  transducer  or  calculations  of  physical  quantities  of 
mass,  gravity,  velocity  and  temperature. 
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AMERICAN  NATIONAL  STANDARD 


A Guide  for  the 
Dynamic  Calibration 
of  Pressure  Transducers 

ANSI  B88.1  - 1972 

[Note:  The  current  number  for  this 
standard  is  ANSI  MC  88.1-1972] 


SECRETARIAT 

THE  AMERICAN  SOCIETY  OF  MECHANICAL  ENGINEERS 


PUBLISHED  BY 

the  AMERICAN  SOCIETY  OF  MECHANICAL  ENGINEERS 
United  Engineering  Center  345  East  47  th  Street  NewYork,N.Y.  10017 
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SENSOR  BEING  CALIBRATED 


PISTONPHONE 

The  pistonphone  is  designed  to  provide  a 
precision  sound  reference  source  at  a fixed 
frequency  and  amplitude  for  calibration  of  low 
pressure  acoustic  sensors.  It  is  usually  pack 
aged  in  a compact,  portable  battery-powered 
configuration.  Calibration  accuracy  is  claimed  to 
be  ± 0.2  dB  with  certain  specified  microphones. 


Pistonphone  Output  124  dB,  250  Hz 

DEAD  WEIGHT  TESTER  USED  AS  A 

PRESSURE  RELEASE  CALIBRATOR 


The  dead  weight  tester  is  designed  and  most 
commonly  used  for  static  calibration  of  DC  type 
pressure  transducers  such  as  strain  gages, 
capacitive,  piezoresistive,  or  quartz  piezotypes 
exhibiting  long  discharge  time  constants. 

The  dead  weight  tester  can  also  be  used 
dynamically  in  the  pressure  release  mode  for 
transducers  whose  output  characteristics  in 
response  to  pressure  rise  and  pressure  drop 
have  been  determined  to  be  the  same. 


DIGITAL 
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PNEUMATIC  PULSE  CALIBRATOR  USING 

FAST-ACTING  VALVES 


0 

1 

1 

ri  ri 

JLI  JLI 

rzi 

^3* 

Pneumatic  Pulse  Calibrator 
PCB  Model  903A 


A number  of  hydraulic  and  pneumatic  pulse 
calibration  devices  have  been  built  using 
manual  or  solenoid  valves.  The  device  at  left 
uses  a digital  pressure  gage  to  set  an  accurate 
static  reference  pressure.  A fast-acting  valve 
quickly  switches  the  pressure  sensor  to  the 
reference  pressure  which  produces  a positive- 
going pressure  pulse  of  known  amplitude  with 
Sms  rise  time.  It  is  usually  more  difficult  to 
produce  a positive-going  step  pulse  than  a 
negative  pressure  drop.  Major  concern  with 
pulse  calibrators  involves  reducing  rarefaction 
waves  in  the  pressure  volume  associated  with 
the  valves. 


Step  Pressure  1 00  psi 
Rise  Time  5 ms 

NIST  TRACEABILITY  - DC  Reference  Gage 


VIBRATING  COLUMN  OF  LIQUID 

This  calibration  method  adapts  vibration  trans- 
fer standard  techniques  to  pressure  calibration. 
An  electrodynamic  shaker  vibrates  a column  of 
fluid  which  generates  a sinusoidal  pressure 
wave.  Output  from  sensor  being  calibrated  is 
compared  with  transfer  standard. 


NIST  TRACEABILITY  - Secondary  Pressure 
Transfer  Standard 
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TEST 

PRESSURE^ 
TRANSDUCER 
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SINUSOIDAL  PRESSURE  GENERATORS 

A small  air  compressor  has  been  modified  to 
generate  a repeatable  sine  wave  output  for 
comparison  calibration  with  a transfer  standard. 
Pressure  amplitude  can  be  adjusted  by  opening 
or  restricting  the  discharge  valve.  Frequency  is 
fixed  at  motor  speed.  This  method  of  calibration 
conveniently  adapts  to  temperature  coefficient 
calibration  by  plumbing  test  transducer  inside 
oven  and  transfer  standard  outside. 

Pressure  Adjustable  3 to  30  psi 
Frequency  Fixed  29  Hz 

NIST  TRACEABILITY  - Dynamically  Calibrated 
Secondary  Transfer  Standard. 


2.0  IN.  DIA.  X 6.0  FT. 
DROP  TUBE 


Pressure  969.3  psi 
Rise  Time  3 ms 
Duration  6 ms 


003A10L0W  NOiSe 
CABLE  ion. LONG 


DROP  PULSE  CALIBRATOR 

The  drop  pulse  calibrator  uses  a tourmaline 
transfer  standard  to  measure  the  hydraulic 
pressure  pulse  from  a mass  impacting  a piston 
and  cylinder  manifold.  The  drop  calibrator  has 
capability  to  generate  a wide  range  of  dynamic 
pressures  from  about  100  to  20,000  psi  with 
relative  ease.  Rise  time  is  about  3ms  and 
pulse  duration  6ms.  Another  manufacturer  offers 
a commercial  high  pressure  version  of  this  type 
of  calibrator  which  uses  calculations  of  mass 
based  on  gravity,  velocity,  and  piston  area  to 
determine  reference  pressure. 


Impulse  Calibrator 
PCB  Model  91 3A 


I 1 

I I 

I I 


012A03  OUTPUT  I OUPI- trace 
CABIT  BNC/BNC  , storage  SCOPE 
/ I (NOT  SUPPIID) 


NIST  TRACEABILITY  - Through 
Dynamically  Calibrated  Model  136A 
Tourmaline  Transfer  Standard. 


lT I 


46ZA36  CHARGE 
AMPIIEIER 
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Instrumented  Shock  Tube,  Model  901 A02 
Atmospheric  Driven  (3  ft.  x 1 1/2  in.  diameter) 


SHOCK  TUBES 

Simple  atmospheric-driven  plexiglas  shock 
tube  generates  low  pressure  shock  waves  for 
resonant  frequency  determination  of  pressure 
transducers  and  gas  passages  associated  wilth 
recessed  diaphragms  or  recessed  mounting. 
NIST  traceability  for  determination  of  transducer 
or  passage  resonance  through  electrically  cali- 
brated measurement  instruments. 


Shock  Response  from 

Accel-compensated 

Transducer 


CALIBRATION  SHOCK  TUBE 


Shock  Response  from  jh©  helium  or  air-driven  shock  tube  is  capable 
Non-compensated 

Transducer  of  producing  shock  waves  with  nanosecond 


rise  times  at  pressures  from  3 psi  using  foil 
diaphragms  to  >1 ,000  psi  with  sheet  aluminum 


diaphragms.  Shock  pressure  calibration 


Incident  Pressure 
134A  Pressure  Bar 
Pressure  123  psi 
Rise  Time  3 )is 


Incident  Pressure 
113A21  ICP  Transducer 
Pressure  5 psi 
Rise  Time  8 


agrees  within  3%  of  other  calibration  methods. 

A cone  shaped  adaptor  installed  at  the  end  of 
the  shock  tube  test  section  amplifies  reflected 
pressures  to  >10,000  psi.  The  cone  adaptor, 
while  not  suited  for  calibration  purposes,  is  use- 
ful for  developing  high  pressure,  high  tempera- 
ture shock  waves  fortesting  thermal  and 
mechanical  shock  capability  of  pressure  sensors. 


HELIUM 


DRIVER 


ALUMINUM 

DIAPHRAGM 


SENSOR 
TIME  OF 
ARRIVAL 


SENSOR 


Calibration  Shock  Tube,  Model  901 A10 
(6  ft.  X 3 in.  diameter) 


Amplitude  of  shock  pressure  step  is  calculated  from: 
Velocity 

Initial  absolute  pressure 
Temperature 
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ARONSON  SHOCKLESS  PRESSURE  STEP 

GENERATOR 

The  "Aronson"  Shockless  Pressure  Step  Gene- 
rator is  a step  toward  combining  capabilities  of 
several  of  the  special  purpose  calibrators  Into 
one  device.  Capabilities  include  the  use  of  high 
accuracy  digital  reference  pressure  standards 
and  fast-acting  poppet  type  valves  to  produce 
positive  or  negative  step  pressures  with  rise 
times  approaching  those  of  shock  tubes. 


The  Aronson  Step  Pressure  Generator  was  developed  by  Phil  Aronson  and 
Robert  Waser  at  the  U.S.  Naval  Ordinance  Laboratory.  Their  primary  objective  was  to 
develop  a device  capable  of  performing  dynamic  step  pressure  calibration  with  greater 
accuracy,  speed  and  ease  than  was  possible  with  the  shock  tube.  Phil  Aronson  dedicated 
much  of  his  professional  career  to  transient  pressure  measurements  and  dynamic  calibra- 
tion. After  acquiring  a license  agreement  with  the  Navy,  PCB  named  the  calibration 
device  after  Phil  who  passed  away  In  the  early  1980's. 

The  Aronson  Step  Pressure  Generator  consists  of: 

Main  Pressure  Reservoir  "A" 

Transducer  Adaptor  Plug  (three  sizes) 

Spring-loaded  poppet  valve  (three  sizes) 

A secondary  pressure  reservoir  "B"  at  poppet  valve/diaphragm  interface 
Impact  weight  for  poppet  valve  release 
Trigger  output 

Precision  Valves  for  adjusting  the  gas  pressure  in  both  reservoirs 

A durable  housing  and  support  system 

Two  each:  digital  reference  pressure  gages  (optional) 

Helium  gas  pressure  source  is  provided  by  the  customer 
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The  concept  and  operation  of  the  Aronson  Step  Pressure  Generator  is  quite  funda- 
mental. It  simply  involves  applying  to  the  transducer  an  accurately  known  static  pressure 
very  quickly.  This  is  accomplished  by  pressurizing  the  main  reservoir  "A"  with  an 
accurately  measured  static  pressure  and  then  exposing  the  transducer  to  the  reference 
pressure  by  releasing  the  quick-opening  poppet  valve.  The  pressure  drop  in  the  main 
reservoir  due  to  the  added  volume  between  the  transducer  diaphragm  and  poppet  valve 
is  negligible  with  flush  diaphragm  pressure  sensors.  Any  pressure  drop  would  be  indica- 
ted by  the  digital  pressure  gage  which  monitors  pressure  in  the  reservoir. 

The  "Aronson"  Step  Pressure  Generator  has  two  separate  pressure  reservoirs 
(A  & B)  in  which  the  pressure  is  controlled  by  accurate  static  digital  pressure  gages.  In 
addition  to  the  main  pressure  reservoir  "A",  pressure  between  the  poppet  valve  and  the 
diaphragm  of  the  transducer  (reservoir  "B")  can  be  pressurized  with  a known  static 
pressure.  Having  control  of  a known  pressure  in  each  reservoir  offers  unique  capability 
to  compare  static  vs.  dynamic  response,  and  provide  small  incremental  pressure  steps  at 
higher  static  levels. 

Rise  time  of  a step  pressure  depends  on: 

1 . type  of  gas  used  (helium  recommended  for  fastest  rise  time) 

2.  diameter  of  poppet  valve 

3.  initial  pressure  difference  across  the  poppet  valve 

4.  design  of  transducer  diaphragm  (flush  or  recess) 


Traceability  to  NIST  is  through  calibration  of  accurate  DC  reference  gages,  used 
to  set  the  known  static  pressure  level  which  the  transducer  under  calibration  is  being 
rapidly  switched  to. 
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SIGNAL 

CONDITIONER 


Digital  Storage  Scope 


TRANSDUCER  BEING 
CALIBRATED 


REGULATED 

PRESSURE 

SOURCE 


Heise  Digital  Reference  Gages 
One  Each  for  Controlling 
Pressure  In; 


'Main  Reservoir  8< 

Poppet  Volume  Reservoir 


mmM 


SHOCKLESS  PRESSURE 
STEP  GENERATOR 


CONSOLE 

(OPTIONAL) 


TYPICAL  CALIBRATION  SYSTEM  USING 
SHOCKLESS  PRESSURE-STEP  GENERATOR 
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ARONSON  SHOCKLESS 

PRESSURE-STEP  GENERATOR 
Model  907A 


A.  MAIN  PRESSURE  RESEVOIR 

B.  POPPET  VALVE  VOLUME 

C.  POPPET  VALVE 

D.  TRANSDUCER  MTG.  ADAPTOR 

E.  DROP  WEIGHT 

F.  IMPACT  PLATE 

G.  PIEZO  TRIGGER 

H.  POPPET  LOCK  TAB 

I.  HOUSING 


3 


MECHANICAL  SYSTEM  DIAGRAM 
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Calibration  and  test  capabilities  of  the  Aronson  Step  Pressure  Generator  include: 


1 .  Dynamic  calibration  of  pressure  trans- 
ducers. Oscillations  of  71  k Hz  at  peak  of 
the  pressure  is  due  to  the  rarefaction 
waves  at  the  interface  of  the  poppet  valve 
and  transducer  diaphragm. 

Dynamic  calibration  of  1 13A24  ICP  Pressure  Trans- 
ducer using  Aronson  Step  Pressure  Generator: 
Pressure  1000  psi 

Transducer  Sensitivity  5.3  mV/psi 


2.  Determination  of  the  RC  time  constant 
of  an  entire  pressure  measurement 
system  including  the  transducer,  signal 
conditioner  and  readout. 

Discharge  Time  Constant  (DTC) 

Time  required  for  a transducer  or  measurement 
system  to  discharge  its  signal  to  37%  of  the  original 
value  from  a step  change  of  pressure. 


3.  Determination  of  response  time  and 
resonance  of  recessed  cavities  and  pas- 
sages. Model  1 13A24,  ICP®  Pressure 
Transducer  recess  mounted  0.050  inch 
with  150  psi  step  change  in  pressure. 

Resonant  frequency  of  0.050  recessed  passage 
Frequency  24.8k  Hz 
Rise  Tiem  71  ps 
Pressure  Step  150  psi 
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4.  Calibration  of  incremental  pressure  steps  above  or  below  preset  static  levels.  It 
is  possible  to  simulate  the  static  pressure  environment  on  the  sensor  diaphragm  that 
might  be  encountered  in  an  underwater  environment  while  applying  a step  pressure 
of  known  amplitude. 


5.  Comparison  of  static  and  dynamic  calibra- 
tion of  a pressure  transducer  using  the 
same  pressure  amplitude  in  the  same  de- 
vice without  moving  the  transducer. 


Slowly  applied  pressure,  150  psi 


Step  pressure  of  same  150  psi 
Rise  time  < 50  |is 


6.  Accurate  dynamic  calibration  of  pressure 
transfer  standards  used  in  other  types 
of  dynamic  calibrators,  e.g.,  sinusoidal  and 
drop-pulse  calibrators. 


Dynamic  calibration  of  1 13A24  ICP  Pressure 
Transducer  using  Aronson  Step  Pressure  Generator 
1000  psi,  5.3  mV/psi 
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7.  Vacuum  step  pressure  calibration. 


Vacuum  step  pressure  calibration,  Model  1 13A24 
ICP  Pressure  Transducer,  -10  psi,  5.34  mV/psi 


8.  Checkout  and  troubleshooting  of  dy- 
namic pressure  transducers,  expecially 
those  used  for  shock  and  blast  wave  mea- 
surements. 


step  pressure  calibration  detects  noise  created  by 
defective  transducer  diaphragm. 


9.  Pulse  calibration  detects  signal  breakup 
due  to  erratic  connection  and  resultant  zero 
shift. 
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During  the  process  of  accumulating  calibration  data  for  this  report,  sensitivity  of 
Model  113A24  ICP®  Pressure  Transducer  was  recorded  and  plotted  using  five  different 
methods  of  calibration.  Sensitivity  deviations  for  the  methods  was  less  than  ±11/2  %. 


Legend: 

A Dead  Weight  Test 
□ Hydraulic  Drop  Test 
X Shock  Tube  Calibrator 
O Aronson  Shockless  Calibrator 
# Pneumatic  Pulse  Calibrator 


Mean  (yj=5. 20942 

Standard  Deviation  (y)=0. 06667 


Calibrations  of  a 113A24  Quartz  Transducer  using  five  different  calibrators. 


While  the  Aronson  Calibrator  has  not  equaled  the  performance  capability  of  the 
"Ultima"  calibrator,  it  is  a step  in  that  direction.  The  digital  storage  oscilloscope  has  signi- 
ficantly improved  the  accuracy  of  dynamic  calibration.  Dynamically  calibrated  pressure 
transfer  standards  have  also  contributed  to  improving  the  accuracy  of  sinusoidal  and 
drop-test  impulse  calibrations. 


0000 
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Improved  Accuracy  in  Piezoelectric  Measurements 

Using  a 

Low  Noise,  Ultra-High  Isolation  Signal  Conditioner 

by 

Bruce  Hartenbaum 
H-Tecb  Laboratories,  Inc. 

California  Institute  of  Technology 

James  Morrison 

Lawrence  Livermore  Laboratory 

John  Cord 
NavDesign 


Commercial  piezoelectric  pressure  transducers  and  accelerometers  are  used  mainly 
in  either  the  voltage  or  charge  mode.  Two  shortcomings  of  existing  commercial 
signal  conditioners  are:  (1)  Their  single-ended  circuitry  has  zero  common  mode 
noise  rejection  capability,  and  (2)  if  the  transducer  test  point  is  at  a 
different  ground  potential  from  data  records,  several  undesirable  conditions  can 
occur,  including  noise  feedthrough,  voltage  offsets  in  the  data  waveforms, 
burnout  of  the  transducer  and/or  the  data  records,  and  shock  hazard. 

As  a result,  noise  contamination  of  piezoelectric  data  often  degrades  the 
accuracy  of  the  data.  In  many  cases  spectral  analysis  is  incapable  of  sorting 
the  true  signal  from  the  noise.  In  other  cases  the  data  are  filtered  with 
attendant  loss  of  signal  components.  The  resultant  data  are  misleading. 

We  have  developed  a new  piezoelectric  transducer  signal  conditioner  that  under 
actual  field  conditions  reduces  noise  contamination  by  20  to  60  db  below  that  of 
conventional  systems.  The  new  signal  conditioner  employs  hybrid  optical- 
electronic  circuitry,  while  maintaining  full  compatibility  with  commercial 
transducers  and  most  data  transmission  systems.  The  new  signal  conditioning  has 
enabled  unusually  clean  data  to  be  recorded  at  full  bandwidth  in  noisy 
environments.  Optical  decoupling  provides  120  db  of  isolation  at  100  Khz  and  up 
to  2500  VRMS  potential  difference  between  measurement  point  and  data  records. 
Differential  conversion  provides  60  db  of  common  mode  noise  rejection  and  allows 
conventional,  twisted  pair,  long  signal  lines  to  be  driven  at  their  bandwidth 
limit,  i.e. , 60  kHz  for  low  cost  cable  and  500  Khz  for  low  loss  cable.  In  noisy 
environments,  the  signal  to  noise  ratio  is  typically  63  db.  An  additional 
improvement  in  measurement  accuracy  is  provided  by  an  end-to-end  precision  square 
wave  voltage  calibrator,  either  under  manual  or  computer  control. 

Used  in  shock  and  vibration  testing,  ordnance  testing,  shock  tube  research,  and 
biomedical  research,  the  new  signal  conditioner  has  been  shown  to  be  quieter  than 
either  conventional  piezoelectric  or  piezoresistive  systems.  An  important 
benefit  to  NIST  is  measurement  accuracy.  An  additional  benefit  to  commercial 
users  is  decreased  set-up  time  and  lower  test  costs. 
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SKILLED  TEST  ENGINEER  CAN 
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TROUBLESHOOTING 
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A NEW  SIGNAL  CONDITIONER 

END  RESULT  OF  DEVELOPMENT 
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SOURCES  OF  NOISE 
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SOURCES  OF  GROUND  VOLTAGE  VARIATIONS  AND 
NOISE  CURRENTS 
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FLOATING  TRANSDUCER 
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NOISE  COUPLING  IN  FLOATING  SYSTEM  VIA  PARASITIC  CAPACITANCES 


MAINTENANCE  AND  TEST  RESPONSE 
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INSULATOR  FAILURE  BNC  INSULATOR  SLEEVE  FAILURE  CABLE  FAILURE/LEAKAGE 
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LOW  NOISE,  ULTRA-HIGH  ISOLATI 
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PROTOTYPE  SIGNAL  CONDITIONER  WITH 
TWO  EQUIPOTENTIAL  GROUND  PLANES 
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HYBRID  ELECTRO-OPTIC/ELECTRONIC 
PIEZOELECTRIC  SIGNAL  CONDITIONER 
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KIEASUREMENT  GROUND  PLANE  | RECORDER  GROUND  PLANE 
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MEASUREMENT  GROUND  PLANE  I ^ RECORDER  GROUND  PLANE 
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MEASUREMENT  GROUND  PLANE  | RECORDER  GROUND  PLANE 
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MEASUREMENT  GROUND  PLANE  | RECORDER  GROUND  PLANE 
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ACTUAL  TEST  WAVEFORM 
NOISE  COMPARISONS 
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ACTUAL  TEST  SHOCK  SPECTRA 
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ACTUAL  TEST  SHOCK  SPECTRA 
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SHOCK  SPECTRA  USING  CONVENTIONAL  SHOCK  SPECTRA  USING  HYBRID 

SIGNAL  CONDITIONING  IS  OUT-OF-SPEC  ELECTRO-OPTIC/ELECTRONIC  SIGNAL 

CONDITIONING  IS  IN-SPEC 
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THREE  CASE  STUDIES  OF  STRAIN/TEST-INDUCED  VOLTAGES 

WHICH  OBSCURED  SIGNALS 

Peter  K.  Stein 

Stein  Engineering  Services,  Inc. 

5602  East  Monte  Rosa 
Phoenix,  Arizona  85018-4646 

Among  the  self-validation  checks  to  which  all  data 
taken  with  impedance-based  transducers  must  be  sub- 
jected, is  a documentation  of  the  self-generating 
(voltage)  responses  from  these  impedance-based  trans- 
ducers. Metallic  or  semi-conductor  strain  gages  and 
strain-gage-based  transducers  such  as  load  cells, 
pressure  transducers  and  accelerometers  are  known  to 
generate  voltages  due  to  strain  (or  due  to  other 
test  conditions).  A simple  way  to  document  the  pre- 
sence or  absence  of  these  noise  levels  is  to  record 
a check  channel  which  is  not  powered  from  an  exter-. 
nal  source  such  as  bridge  supply  voltage  or  current. 
It  is  also  possible  (but  not  in  transient  testing) 
to  turn  off  such  a supply  during  the  test,  on  a 
single  channel.  The  effect  is  DYNAMIC  only  not  static. 

Three  case  studies  have  recently  come  to  the  author's 
attention  and  are  described  below. 

Case  Study  1:  Charles  P.  Wright,  Manager,  Measure- 
ments Engineering  Group,  TRW  S & T Div. 

"A  shock  separation  test  was  run  at  McUonnell- 
Douglas  in  Huntington  Beach,  California  in  Feb- 
ruary 1990.  We  had  1/8“  constantan  quarter 
bridges  mounted  on  the  structure  to  check  the 
survival  of  the  installation  technique  under 
high  shock  loadings..  The  shock  acceleration 
levels  present  at  the  gage  locations  was  about 
10,000  Gpk  in  the  direction  normal  to  the  gage 
pattern.  Levels  in  the  other  two  directions 
were  not  measured. 

“In  one  area  we  had  duplicate  gages  mounted 
right  next  to  each  other  - same  gage,  same  di- 
rection, same  installation,  same  technician, 
same  materials,  same  everything.  One  was  powered 
at  5 volts  DC  in  a Dynamics  7600A  Signal  Condi- 
tioners and  recorded  on  analog  tape  DC-75KHz 
bandwidth  (-3dB).  That  is  Fig.  10.  The  first 
120  milliseconds  of  the  record  are  shown  in  the 
top  part. ..The  lower  left  is  the  first  20  milli- 
seconds to  an  expanded  scale... 

"Fig.  11  is  the  same  data  for  the  gage  next 
door  run  with  the  excitation  OFF.  Beautiful, 


isn't  it!"  (WRIGHT  1990)  TRW  Space  & Technology 
Note  that  the  strain-induced  voltage  noise  is  more 
than  10%  of  signal.  Such  noises  cannot  "be  swept 
under  the  rug".  As  in  the  Pratt  & Whitney  Case  Study, 
they  are  the  size  of  the  rug! 

Case  Study  2 : In  1987,  Richard  BILLIA  was  involved 
in  a test  on  a i-scale  model  test  vehicle  in  Sandia 
National  Laboratory's  Albuquerque  19  ft  diameter 
blast  tube  (see  photo).  Billia  is  with  Lawrence 
Livermore  National  Labs.  He  reported  (BILLIA  1987) 
the  use  of  constantan  CEA-062-LIT-350  Measurements 
Group  strain  gages  on  a cylinder  with  the  use  of 
parallel  gages  again  - one  powered  and  one  un-powered. 
Figure  12  shows  the  two  records.  Test  - induced  volt- 
ages are  again  some  10%  of  signal. 

Case  Study  3 : Rodney  MAY  of  Sandia  National  Labora- 
tory,  Albuquerque  reports  in  MAY  (1991): 

"Along  with  the  Midgetman  Program  on  which  Sandia 
was  working,  the  19-ft  diameter  blast  tube  was 
being  expanded.  Before  the  full  test,  one  with  a 
i-charge  was  run.  Strain  gages  were  put  on  the 
blast  tube  to  see  if  the  stress  calculations  had 
been  correct. 

"They  were  getting  near  75%  of  yield  stress 
and  were  worried  that  the  tube  would  fail  at 
full  charge.  I asked  them  if  they  had  run  any 
noise-check  gages.  Their  reply  was  a classic: 
'No,  no  ! We  needed  all  the  channels  available 
to  get  our  data! ' 

"A  couple  of  days  later  we  put  2 check  chan- 
nels on:  one  un-powered  gage  mounted  to  the  tube, 
and  one  un-powered  and  un-strained,  on  a 'dummy' 
block.  That  way  we  could  observe  the  strain- 
induced  voltages  and  any  antenna  effects.  After 
the  shot  they  showed  me  the  results,  but  only 
briefly,  and  documentation  was  ever  issued. 
It  turned  out  that  80%  of  the  signal  they  had 
obtained  were  due  to  antenna  effects,  which  the 
un-powered,  un-strained  gage  documented  very 
nicely. 

"Noise-check  channels  do  give  you  information, 
but  you  have  to  do  battle  for  them!" 


PHOTO;  SANDIA  NATIONAL  LABORATORY'S  ALBUQUERQUE  19  FT 

DIAMETER  BLAST/SHOCK  TUBE.  An  explosive  charge  deton- 
ated at  the  far  end  ^ the  shock  tube  produces  a 
shaped  blast  wave  used  in  testing  responses  of  a 
structure  under  such  a wave. 


REFERENCES: 


All  3 case  studies  are  reported  in: 
Proc.  Western  Regional  Strain  Gage 
Committee  Meeting,  February  5-6,  1991, 
Tempe,  AZ.  From:  42  N.  Benson  Ave, 
Upland,  CA  91786. 

WRIGHT  (1990);  Letter  of  March  21  to 
Peter  Stein  - also  reporting  Case 
Study  1. 

The  cause  for  these  strain-induced 
voltages  is  not  as  yet  understood  but 
has  been  found  in  all  metals,  especial- 
ly strain  gages.  Nickel  (as  in  leads 
and  resistance  thermometers,  and  in 
Alumel,  as  in  thermocouples.  The  author 
has  collected  many  references  to  this 
phenomenon. 
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FIGURES  10  AND  11:  NOISE-DOCUMENTATION  DURING  A SHOCK-SEPARATION  TEST 
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FIGURE  11;  SELF-GENERATING  NOISE  LEVELS.  Same  signal  conditioning  but 
SUPPLY  VOLTAGE  IS  ZERO.  At  least  3 sources  are  suspected:  self-gen- 
erating ability  of  strain  gage;  rapid  gage  grid  movement  in  ambient 
magnetic  fields;  charge-response  of  gage  backing.  Test  planner:  TRW 
& T Group,  Measurements  Engineering  Dept.  Charles  Wright,  Manager. 
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From  WRIGHT  (1990) 
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FIGURE  12:  NOISE  DOCUMENTATION  IN  SANDIA  LABORATORY'S  ALBUQUERQUE  19  FT  DIAMETER  BLAST  TUBE 


From  BILLIA  (1987) 
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Reprinted  with  permission  from 
Measurements  and  Control,  (146)  April  1991 

CALIBRATION  OF 
BALLISTIC  PRESSURE 
TRANSDUCERS 


ARP  AD  A.  JUHASZ,  CHARLES  D.  BULLOCK  AND  DONALD  H.  NEWHALL 


Techniques  for  calibrating  ballistic  pressure  transducers  include  the  ballistic 
pulse  method,  shock-tube  methods  and  pressure-step  techniques.  A positive- 
going, high-pressure-step  generator  is  useful  in  calibration,  servicing  and 
development  of  dynamic  high-pressure  transducers.  A positive- going  step  cali- 
bration device  allows  accurate,  safe  and  simple  dynamic  calibration  and  evalua- 
tion of  ballistic  pressure  gages.  The  gage  response  obtained  can  be  related  to 
its  static  ( deadweight)  behavior.  Several  gages  may  be  evaluated  simultaneous- 
ly, relative  to  a common  dynamic  event.  The  calibrator  may  also  be  used  as  a 
diagnostic  tool  in  analyzing  and  developing  experimental  pressure  gages. 


.Ajthovigh  ballistic  pressxire  transducers  are 
used  to  meastire  events  which  occur  in  milli- 
seconds, the  determination  of  their  response 
characteristics  has  been  limited  to  static  cali- 
bration against  a deadweight  pressure  stand- 
ard. Although  the  DWT  is  a primary  standard, 
the  assumption  is  made  that  the  static  and 
dynamic  responses  of  the  gage  in  question  are 
equivalent.  Obviously,  any  differences  in  gage 
response  to  static  and  dynamic  events  can  lead 
to  serious  measurement  errors.  There  has 
been  a general  agreement  that  dynamic  techni- 
ques are  needed  to  supplement  current  static 
calibration  methods,  ^veral  techniques  deve- 
loped to  address  this  problem  are  discxissed^’^ 
— (1)  the  ballistic  pulse  method,  (2)  shock-tube 
methods,  (3)  negative-going  pressure  step  and 
(4)  positive-going  pressure  step, 

BALLISTIC-PULSE  METHOD 

In  one  version  of  the  ballistic-pulse  techni- 
que, the  gage  is  mounted  at  the  end  of  a tube  in 
contact  with  a hydraulic  fluid  confined  by  a 
movable  piston.  The  tube  guides  a projectile 
which  impacts  the  piston  to  create  a pressure 
pulse  in  the  fluid.  Different  pressures  may  be 
achieved  by  varying  the  compressibility  of  the 
fluid,  the  mass  of  the  piston,  and  the  mass  and 
velocity  of  the  projectile.  FVom  these  variables 
the  reference  pressure  may  be  calculated. 

The  pulses  rise  within  milliseconds,  and 
mimic  the  characteristic  rising  and  falling  of  a 
ballistic  pressure  pvilse. 


One  such  device,  capable  of  operating  to  a 
pressure  of  100,000  psi,  is  operational  at  the 
Combat  Systems  Testing  Activity  (CSTA), 
Aberdeen  Proving  Ground,  Maryland.^  This 
method  is  useful  only  for  dynamic  comparison 
of  different  pressure  gages  because  variations 
in  projectile  velocity,  frictional  effects  on  the 
moving  piston,  and  other  energy  losses  make  it 
difficult  to  accurately  compute  the  actual  deli- 
vered pressures.  Because  a projectile  is  fired 
during  the  calibration  process,  this  method  re- 
quires extensive  safety  provisions,  not  readily 
available  in  most  laboratories. 

SHOCK-TUBE  METHODS 

There  are  two  general  approaches  to  shock- 
tube  calibration.  In  the  first,  the  test  gage  is 
movmted  in  the  end  wall  of  a tube  and  sub- 
jected to  a reflecting  shock  wave.  The  gage 
output  is  monitored  as  the  shock  front  arrives 
at  and  reflects  from  the  end  wall. 

In  the  second  approach,  the  gage  is  moimted 
in  the  side  wall  of  the  tube  and  its  output  is 
monitored  as  the  shock  front  passes.  Both 
methods  generate  rapidly  rising  pressure  pul- 
ses that  are  readily  calcvilated  by  thermody- 
namic principles  vising  velocity  measurements 
and  known  gas  properties.  Shock-tube  me- 
thods are  useful  in  establishing  the  dynamic 
resp>onse  characteristics  of  pressure  gages. 
However,  calibration  is  generally  limited  to 
pressures  below  1000  psi,  whereas  ballistic  ap- 
plications require  much  higher  pressures. 
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FIGURE  2.  Pressure-Time  History  of  a 75,000  psi 
dynamic  pressure  step;  10-second  window. 


NEGATIVE-GOING  PRESSURE-STEP 
METHOD 

In  this  technique,  the  gage  is  exposed  to  a 
high  pressure  under  static  conditions,  using  a 
hydraxilic  fluid.  The  pressure  on  the  gage  is 
then  relieved  using  a fast-acting  dump  valve, 
bringing  the  system  to  atmospheric  pressure. 
The  response  of  the  gage  output  obtained  dur- 
ing the  depressurization  is  assumed  to  be  the 


inverse  of  the  corresponding  positive  pressure 
step. 

This  method’s  strengths  include  its  relative 
simplicity  and  suitability  to  use  in  calibration 
facilities,  and  qviick  response,  typically  100 
microseconds  or  less.  However,  the  major  as- 
sumption, that  the  positive  response  of  the 
gage  is  equal  and  opposite  to  the  negative 
response  of  the  gage,  is  not  completely  ac- 
curate; pressure  preloading  of  the  gage-to- 
mount  interface  and  hysteresis  can  cause  sig- 
nificant differences  between  the  magnitudes  of 
the  pressurization  and  depressurization  pul- 
ses. 

POSITIVE-PRESSURE-STEP  METHOD 

Figure  1 illustrates  the  basic  configuration 
of  the  positive-pressure-step  device.  A large 
pressure  reservoir  is  connected  to  a much 
smaller  test  chamber  by  a fast-acting  ball 
valve.  The  test  chamber  is  equipped  with  seve- 
ral gage  ports  and  a vacuum  port  which  aids 
in  filling  the  reservoir  and  setting  the  baseline 
pressure  in  the  test  chamber.  Reservoir  pres- 
sure is  provided  by  a (Harwood)  hydraulic 
high-pressure-generating  panel.  The  pressure 
is  monitored  by  a static  reference  gage  whose 
calibration  is  traceable  to  NIST. 

The  large  ratio  of  reservoir  to  test-chamber 
volume  minimizes  the  overall  system  pressure 
drop  while  generating  the  pressure  step,  there- 
by reducing  stress  loading  on  the  system  com- 
ponents (particularly  the  static  reference 
gage).  The  maximum  operating  design  pres- 
sure is  150,000  psi,  and  the  action  time  (10- 
90%  of  peak  pressure)  is  well  under  one  mil- 
lisecond. Outputs  of  both  the  static  reference 
gage  and  the  ballistic  test  gages  are  monitored 
during  the  test.  The  final  steady-state  output 
of  the  reference  gage  is  taken  as  the  true  value 
of  the  pressure-step  maximum.  Short-  and 
long-term  monitoring  of  the  test  gage  outputs 
establish  a relationship  between  dynamic  and 
steady-state  responsive  behavior. 

PERFORMANCE  ANALYSIS 

Examples  of  both  short-  and  long-term  res- 
ponses of  a commercial  piezoelectric  pressure 
gage  to  a positive-going  pressure  step  are  pre- 
sented in  Figures  2 through  5.  Figure  2 shows 
a typical  pressure  versus  time  history  for  a 
75,000  psi  pulse  acquired  over  10  seconds. 
Monitoring  the  test  gage  response  for  10  sec- 
onds allows  correlation  to  the  steady-state  res- 
ponse of  the  static  reference  gage.  Figure  3 
shows  a 20-millisecond  window  of  the  same 
75,000  psi  pulse;  Figures  4 and  5 present  pro- 
gressively shorter  windows. 
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Overlaying  the  traces  indicates  that  the 
steady-state  vsQue  of  the  pressure  step  is  rap- 
idly achieved  and  held,  after  initial  oscillations 
have  died  out.  Similar  oscillations  have  been 
reported  in  lower-pressure  devices  by  both 
Smith  and  Dykstra^’®;  one  may  conclude  that 
these  oscillations  are  caused  by  actions  within 
the  pressure-generating  system  and  the 
mount,  not  by  the  gage  itself. 

One  important  application  of  the  step  cali- 
brator is  comparing  the  responses  of  different 
types  of  pressure  gages  to  the  same  input. 
Figures  6 and  7 show  pressure  versics  time 
histories  of  both  a piezoelectric  gage  and  a 
strain-gage  measxiring  a positive-going  100,000 
psi  pressure  step.  The  traces  are  virtually 
identical,  exhibiting  initial  system  oscillations 
which  decay  quickly. 

The  calibrator  may  also  be  used  to  analyze 
the  behavior  of  experimental  gages  by  compar- 
ing their  output  with  a known  standard.  Fig- 
ure 8 shows  pressure  versus  time  history  of  a 
developmental  pressure  gage  subjected  to  a 
125,000  psi  pressure  pulse.  This  curve  ex- 
hibits, as  does  the  known  standard,  a smooth 
pressure  rise  and  stable  output  after  the  peak 
pressure  has  been  attained. 

The  positive-step  pressure  generator,  in 
copjvmction  with  conventional  deadweight  cal- 
ibration methods,  can  be  useful  in  tracing  dy- 
namic gage-response  problems.  Figure  9 ex- 
hibits the  pressure  history  of  a developmental 
pressure  gage  exposed  to  a 130,000  psi  pres- 
svire  step.  The  trace  indicates  a clear  upward 
drift  after  the  step  is  complete.  The  same  gage 
had  exhibited  good  response  behavior  on  static 
calibration  using  the  deadweight  system.  In 
gim  tests,  however,  the  device  read  5%  low. 
From  the  dynamic  calibration  tests,  it  became 
clear  that  it  took  tens  of  milliseconds  for  the 
gage  output  to  reach  a maximum.  This  could 
not  be  detected  by  a static  calibration  techni- 
que. On  dissecting  the  gage,  it  was  found  that 
the  bond  between  the  strain  patch  and  the 
gage  body  was  faulty.  A defect  of  this  type 
could  account  for  the  differences  in  static  ver- 
sus dynamic  behavior  noted  in  testing. 

The  equipment  and  method  of  dynamic 
high-pressure  pulse  generation  for  calibration 
and  test  purposes  has  the  following  advan- 
tages: 

1.  Direct  on-site  calibration  traceable  to  NIST. 
No  involved  calibration  calculation  models  are 
required. 

2.  Calibration  is  not  affected  by  changes  in  test- 
chamber  volume. 

3.  The  equipment  does  not  require  special  foun- 
dation provisions. 

4.  The  equipment  does  not  require  special  align- 
ment during  instellation. 
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FIGURE  3.  Pressure-Time  History  of  a 75,000  psi 
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FIGURE  4.  Pressure-Time  History  of  a 75,000  psi 
dynamic  pressure  step;  5-millisecond  window. 


FIGURE  5.  Pressure-Time  History  of  a 75,000  psi 
dynamic  pressure  step:  2-millisecond  window. 
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FIGURE  9.  Response  of  a Faulty 
Developmental  Pressure  Gage  to  a 
130,000  psi  positive  pressure  step. 


4.  The  equipment  does  not  require  special  align- 
ment during  installation. 

5.  Mechanical  parts,  critical  to  accurate  calibra- 
tion, are  not  subjected  to  severe  impact  shock 
or  stress. 

6.  Operation  does  not  involve  expendable  parts  or 

supplies. 

7.  Peak  pressure  can  be  maintained  to  monitor  the 
final  steady-state  output  of  the  test  gages. 

8.  A negative-going  pressure  pulse  is  available  as 

an  option. 

The  design  parameters  of  the  equipment 
can  be  varied  easily  so  that  greater  accuracies 
and  still  higher  pressures  can  be  attained.  The 
device  can  be  turned  to  knotty  problems  in 
internal  ballistics  as  they  emerge.  Similar 
equipment  could  be  useful  in  fundamental  stu- 
dies in  physics,  such  as  thermodynamics,  reac- 
tion and  diffusion  rates,  and  physical  proper- 
ties of  materials  under  dynamic  pressure  load- 
ing. The  speed  of  compression  loading  is  fast 
enough  to  approach  isentropiccompressions.Q 
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FIGURE  10.  Dynamic 
Positive  Step  Pressura  Gen- 
erator. Trunion  rrKHjnted 
pressure  vessel  assembly 
with  pressure  generatbn 
control  panel. 
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DYNAMIC  PRESSURE  CALIBRATOR 


(drop  hammer  method) 


J.  Riegebauer 


AVL  List  GmbH  Graz,  Austria 


Abstract: 

The  dynamic  pressure  calibrator  functioning  according  to  the  principle 
of  the  drop  hammer  method,  allows  the  generation  of  pressure  curves, 
which  are  similar  in  amplitude  and  time  (10  - 800  MPa,  2.5  msec)  to 
phenomena  e.g.  occurring  in  ballistic  measurement  techniques.  The 
paper  explains  the  set-up  of  the  calibrator  and  the  possible  error 
influences  on  the  peak  pressure  are  discussed  in  detail.  A calibration 
method  for  obtaining  traceability  is  introduced  and  prospects  as  to 
further  development  work  for  reaching  primary  standard  are  outlined. 
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Introduction 


The  course  of  the  gas  pressure  is  an  important  characteristic  magnitude  for 
engine  indication  and  internal  ballistic  measurements  techniques.  In  the  past  few 
years  piezoelectric  pressure  transducers  have  proved  to  be  especially  reliable 
for  measuring  the  course  of  this  characteristic  magnitude. 

Based  on  its  high  resonance  frequency  (100  - 280  kHz),  the  piezoelectric 
pressure  transducer  is  well' suited  for  measuring  these  highly  dynamic  events 
(ms-range).  It  is  of  utmost  importance  to  exactly  know  the  transducer 
characteristics,  in  particular,  its  sensitivity  behavior  for  these  measurement 
requirements.  For  obtaining  the  sensitivity  characteristics  a calibration  method  is 
to  be  applied  which  largely  corresponds  to  the  practical  measurement  procedure 
(see  [3]  fig.  4). 


From  the  above  the  following  requirements  can  be  derived  for  a dynamic 
pressure  calibrator: 

calibration  curve  similar  in  amplitude  and  time  e.g.  up  to  800  MPa, 

2-20  msec  to  the  phenomenon  to  be  measured 

short  load  duration  positively  effects  life  time  of  transducers 

high  reproducibility  and  accuracy 

quick  and  easy  operation 

to  make  the  calibration  method  traceable 

to  become  primary  standard  for  dynamic  pressure 

The  present  classical  calibration  method  with  a static  dead  weight  tester  can  not 
fulfill  the  requirement  on  the  calibration  procedure  to  correspond  to  the  highly 
dynamic  measurement  events. 

However,  the  dead  weight  tester  is  a primary  standard  and  the  required  pressure 
can  be  generated  very  precisely,  which  is  an  essential  advantage. 

During  the  past  few  years,  great  efforts  have  been  made  by  various  institutes  to 
also  obtain  a primary  standard  for  dynamic  pressure.  So  far,  this  aim  has  not 
been  reached. 

The  following  different  methods  are  known  for  generating  dynamic  pressure: 

ballistic  pressure  generators  (e.g.  air  gun,  closed  vessel) 
the  principle  of  the  shock  wave  tube 
negative  or  positive  going  pressure  step  (^isec-range) 
drop  hammer  methods  etc. 

From  these,  AVL  chose  the  drop  hammer  method  for  developing  the  dynamic 
pressure  calibrator. 
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Principle  of  pressure  generation  and  calibrator  set-up 


The  impact  mass  M drops  from  the  falling  height  (see  fig.  1)  on  the  piston  m, 
which  seals  the  hydraulic  pressure  chamber.  The  piston  and  pressure  chamber 
assembly  represent  the  measuring  head  MH  equipped  with  the  pressure 
transducers  to  be  calibrated.  At  the  impact  of  the  impact  mass  the  kinetic  energy 
is  converted  into  compression  energy  by  ingressing  of  the  piston  into  the 
pressure  chamber.  During  transition  from  the  compression  phase  into  the 
expansion  phase  the  maximum  pressure  p is  reached. 

max 


feUijg 
Epol.  = 


H 


1 p 

— M * V ^ g*h*M 

2 


M — falling  ttbss 
g — earth  gravity 
V — irrpact  velocity 
A — piston  area 
Epol.  - potential  erergy 
E^jj^  - kiretic  energy 


^kom.- 
X 

H - 


ocnpression  energy 
piston  travel 
felling  leigfe 


pressure  chamber 
measuring  head 


foondaticn 


Fig.  1 Principle  of  pressure  generation 
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Set-up  of  the  calibrator 


B621-FU  B621-CS  COMPUTER  PLOTTER 


Fig.  2 Complete  set-up  of  the  B621  DPC  Dynamic  Pressure  Calibrator 


Fig  3 shows  all  essential 
B621DPC  components 
and  lists  the  most 
important  magnitudes 

influencing  the  pressure 
generation 


Fig.  3 Schematic  view  of 
magnitudes  influencing 
the  pressure  generation 
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Betow  some  of  these  magnitudes  will  be  discussed  in  detail. 


B621  FU  Falling  Unit 


Fig.  4 Schematic  arrangement  of  the  B621  FU  Falling  Unit 
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Requirements  on  the  place  of  installation 

When  choosing  the  place  of  installation  the  following  conditions  have  to  be  met: 

the  concrete  foundation  has  to  be  made  on  solid  ground 

the  place  of  installation  must  not  have  cellaring  or  rooms  lying  underneath 

the  room  must  be  dust-free  and  free  of  combustion  gases 

room  temperature  20'  C +/-  T C 

total  space  requirement  approx.  3.3  m x 3.0  m 

power  supply  220  Volt,  50  Hz  or  1 1 0/120  Volt,  60  Hz 

compressed  air  5 - 12  bar  (70  - 170  psi) 


Influence  of  gravitational  acceleration  at  the  place  of  installation  on  the 
peak  pressure: 

The  impact  mass  M drops  from  the  falling  height  H due  to  gravitational 
acceleration  in  first  approximation  in  free  fall  and  reaches  the  impact  velocity  v^ 
directly  before  hitting  the  piston.  This  velocity  and  the  weight  of  the  impact  mass 
depend  on  the  gravitational  acceleration.  The  influence  of  the  gravitational 
acceleration  at  the  place  of  installation  on  the  falling  height  H for  reaching  the 
same  impact  velocity  as  at  the  place  of  standard  (45°  latitude,  0 m above  sea 
level,  comparable  with  Graz,  Austria)  is  shown  in  fig.  5. 


Fig.  5 Influence  of  gravitational  acceleration 
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Requirements  on  the  concrete  foundation 

The  B621  DPC  is  positioned  on  a concrete  foundation  with  cast-in  steel  plate. 
The  foundation  acts  as  counter  mass  during  the  impact  of  the  impact  mass  and 
makes  the  impact  procedure  and  thus  the  pressure  generation  largely 
independent  from  installation  conditions  such  as  nature  of  ground  composition, 
stiffness  and  stability  (see  fig.  6). 


Fig.  6 Constructional  drawing  of  concrete  foundation 


At  the  preparation  of  the  foundation  the  geometric  dimensions  can  only  be  kept 
within  certain  tolerances.  These  deviations  result  in  a dispersion  of  the  total 
mass  by  the  nominal  value  and  thus  in  a dispersion  of  the  peak  pressure  which 
can  be  considered  as  systematic  error  from  calibrator  to  calibrator. 
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The  influence  of  the  total  mass  of  the  foundation  on  the  peak  pressure  and  on 
the  ingression  depth  of  the  piston  is  shown  on  fig.  7. 


X [mm] 


M [kg] 


A p/p(0)  [%] 


Fig.  7 influence  of  foundation  mass  on  the  peak  pressure 

The  stiffness  of  the  ground  on  which  the  concrete  foundation  is  installed  also 
influences  the  peak  pressure.  Therefore  make  sure  that  the  place  of  installation 
has  no  cellarings  or  rooms  lying  underneath  (change  in  peak  pressure  of  approx. 
5-10%). 
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Mass  of  the  impact  mass 


Manufacturing  tolerances;  AM  / Mq  < +/-  0.140% 

Weighing  at  the  Bureau  of  Standards:  AM  / Mg  < +/-  0.001  % 


Mass  of  the  built-in  impact  mass 
including  trigger  cable 


AM  / Mg  < +/-  0.080  % 
Ap  / Pg  < +/-  0.040  % 


Height  adjustment  and  zeroing 


When  a pressure  value  is  selected  it  is  converted  to  the  relevant  falling  height. 
The  impact  mass  is  positioned  by  means  of  a stepping  motor  (400  steps  per 
revolution)  and  the  electronic  control  unit  at  this  falling  height.  The  falling  height 
is  calculated  beginning  from  the  piston  surface  i.e.  this  position  has  to  be 
determined  after  switching  on  the  system  (zeroing).  The  positioning  repeatability 
of  the  impact  mass  to  the  falling  height  selected  is  +/-  0.02  mm,  this  value  was 
verified  by  a laser  interferometer.  A falling  height  of  14  mm  corresponds  to  a 
relative  uncertainty  of  the  height  adjustment  or  the  pressure  generation  of 


AH/Hg  = +/-  0.15% 
Ap  / Pg  = +/-  0.08  % 


This  uncertainty  is  less  than  +/-  0.1  % for  falling  heights  over  300  mm. 
(Ap  / Pg  = +/-  0.03  % ) 
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Measuring  the  impact  velocity 


A light  screen  system  measures  the  characteristic  time  during  falling  of  the 
impact  mass.  By  means  of  this  time  the  falling  behavior  of  the  impact  mass  can 
be  checked  and  extrapolated  to  the  impact  velocity,  (see  fig.  8) 

The  uncertainty  at  measuring  the  impact  velocity  is  approx.  +/-  0.4  % 


0 100  200  300  400  500  600  700 

Falling  Height  [mm] 


Fig.  8 Measurement  uncertainty  of  impact  velocity  v^ 


Influence  of  friction  between  guiding  rods  and  guide  sleeve  of  the  impact 
mass  and  of  air  friction  on  the  impact  velocity  v^ 


When  the  measured  velocities  are  compared  with  the  theoretically  calculated 
velocities  it  deviates  by  approx.  +1-2 -b  %. 

Out  of  that  the  air  friction  has  a quota  of  approx.  +/-  0.03  %.  As  the  velocities  are 
measured  only  the  measurement  uncertainty  is  relevant  for  error  considerations 
(see  fig.  8). 
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B&21  MH/HP  Measuring  Head  / High  Pressure 


piston  gage 
I 


locking  coupling 


bleeding  coupling 


w 


bleeding  channel 


sealing  surface 

filling  chamber 
filling  channel 


sealing  surface 

high  pressure  chamber 
transducer  vent 

0-rings 


■ sealing  screw 


Fig.  9 Set-up  of  the  High  Pressure  Measuring  Head 
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Pressure  curves  generated  with  B621 DPC 


Time  [msl 

Fig.  1 0 Pressure  curves  generated  with  B621  DPC 

Influence  of  the  measuring  head  volume  on  the  peak  pressure 

The  peak  pressure  generated  depends  on  the  measuring  head  volume,  the 
dependence  is  shown  on  fig.  1 1 (experimentally  determined). 

KH  [mm] 


p/p(0) 

Fig.  1 1 Dependence  of  peak  pressure  on  measuring  head  volume 
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When  manufacturing  the  measuring  head  it  is  therefore  necessary  to  specify  the 
exact  pressure  transducer  type  and  to  know  the  dimensions  of  the  front  part.  Any 
deviation  or  change  of  the  original  volume,  e.g.  use  of  different  sealing  rings  (not 
in  line  with  the  specification),  thermal  shock  diaphragms,  mechanical  filters  etc. 
leads  to  a systematic  error. 

The  uncertainty  of  the  measuring  head  volume  due  to  manufacturing  tolerances 
at  known  pressure  transducer  types 

AV  = +/-  20  mm3 
AV/Vq  = +/-2.4% 

Ap  / Po  = +/-  0.7  % 

corresponds  to  an  uncertainty  piston  height: 

AKH  = +/-0.26mm 

When  measuring  the  manufacturing  dimensions  with  a micrometer  the 
uncertainty  can  be  reduced  to 


AV  = +/-  2.6  mm3 
AV/Vq  = +/-0.30% 

Ap  / Po  = +/-  0.09  % 

When  the  measuring  head  had  already  been  calibrated  the  piston  gage  always 
ensures  the  same  piston  height.  The  uncertainty  of  this  method  for  adjusting  the 
piston  height  with  this  piston  gage  is 


AKH  =+/-0.02mm 
A V = +/- 1.6  mm3 
A V / Vo  = +/_  0.20  % 
Ap/po  = +/-0.05% 


Dead  volume  of  pressure  transducers 

Pressure  transducers  with  the  same  installation  dimensions  can  have  a different 
dead  volume  (see  fig.  12)  in  the  front  part,  e.g.  the  transducers  of  the  type 
AVL  5QP6000T,  KISTLER  621 1 , KISTLER  6203  have  a different  dead  volume  of 
approx.  4 mm3- 
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When,  for  example,  these  transducers  are  exchanged  during  measurement 
series  without  considering  this  fact,  it  will  cause  a relative  change  of  the 
measuring  head  volume  or  the  peak  pressure  by  approx. 


AV  / Vo  = +/-  0.47  % 
Ap  / Po  = +/-  0.13  % 


transducer 


dummy 


— r ^ 

nj 


Fig.  12  Dead  volume  of  pressure  transducers  or  dummies 


Use  of  dummies 

Should  the  dummy  not  have  exactly  the  shape  of  the  pressure  transducer  it  will 
result  in  a relative  change  of  the  measuring  head  volume  or  the  peak  pressure 
by  approx. 


AV/Vo=  +/-1.0-1.4% 
Ap  / Po  = +/-  0.3  - 0.4  % 


A further  change  in  peak  pressure  derives  from  the  different  stress  behavior 
between  pressure  transducers  and  dummies  (not  important  at  the  static 
calibration)  which  can  be  proven  experimentally  (see  fig.  13).  Out  of  four 
pressure  transducers  of  the  same  type  one  was  defined  as  reference  pressure 
transducer  and  a measurement  series  was  taken.  In  the  2nd  measurement 
series,  the  three  pressure  transducers  were  replaced  by  three  dummies  with 
their  geometric  dimensions  exactly  corresponding  to  the  pressure  transducer 
shapes.  The  reference  pressure  transducer  was  not  dismounted  from  the  bore. 
The  difference  of  the  results  from  the  two  series  was  higher  than  +/-  0.8  %.  The 
simplified  model  of  the  pressure  transducer  and  the  dummy  shown  on  fig.  13 
was  used  for  estimating  the  influence  of  stress  property  (Young's  modulus  E). 
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The  results  can  be  seen  on  table  1 . 


Fig.  13  Simplified  model  of  pressure  transducer  and  dummy  and  schematic 
view  of  experimental  set-up 

Table  1 E... Young's  modulus,  x displacement 

E/GPa  x/pjn  5x/Xo/% 


85.000  0.992  0.0 

84.575  0.996  0.4 

84.150  1.001  0.9 

80.750  1.037  4.4 

Uncertainty  when  determining  the  piston  area 


AA/Ao=  2-AD/Do=  2-0.005/10  = +/-  0.001  % 


This  value  is  true  for  the  piston  not  exposed  to  pressure. 
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Contaminated  contact  surface  of  the  measuring  head  bottom 

When  the  contact  surface  between  measuring  head  bottom  and  base  bottom  of 
the  B621FU  is  not  cleaned,  e.g.  from  drops  of  liquid,  chip  grains,  etc.  it  will  lead 
to  a peak  pressure  error.  Fig.  14  shows  this  error  (experimentally  determined). 

della  p/p(0)  [%) 


Fig.  14  Inclined  contact  surface  of  the  measuring  head 


Air  bubbies  in  the  high  pressure  chamber  of  the  measuring  head 

Air  bubbies  influence  the  shape  of  the  generated  pressure  curve  and  its  peak 
pressure.  Furthermore,  the  reproducibility  of  zeroing  is  deteriorated  which  is  an 
additional  uncertainty  for  the  height  adjustment.  Constructional  measures  on  the 
filling  station  of  the  measuring  head  largely  ensure  bubble-free  filling. 


Influence  due  to  leakage  of  the  measuring  head  on  its  volume 

100  repetitions  at  100  MPa  and  at  600  MPa  were  made.  The  piston  height  was 
measured  with  a micrometer  before  starting  the  measurement.  After  completion 
of  the  measurement  series  a difference  of  50  )jjn  in  the  piston  height  was 
measured,  this  corresponds  to  a relative  change  in  volume  or  pressure  of 


AV/Vo=  +/-0.4% 
Ap/po  = +/-0.1  % 
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As"  approx.  15  repetitions  are  made  at  a typical  calibration  cycle,  this  error 
influence  is  very  low  i.e.  less  than  +/-  0.015  % . 

Influence  of  the  ambient  temperature  on  the  measuring  head 

Change  of  piston  area 

A temperature  change  of  2°C  results  in  a relative  change  of  the  piston  area  of 
approx.  0.005  % . 

Change  of  measuring  head  volume 

A temperature  change  of  2°C  results  in  a relative  change  of  the  measuring  head 
volume  of  approx.  0.007  % . 

Change  of  compressibility 

A temperature  change  of  2°C  results  in  a relative  change  of  the  compressibility  of 
the  high  pressure  liquid  of  approx.  0.6%  . 

The  total  influence  of  the  temperature  change  on  the  peak  pressure  can  be  seen 
on  fig.  15  (experimentally  determined) 

Ap  / Po  = +/-  0.1 5 % per  -/+  1 °C 


delta  p/p(0)  [%] 


Fig.  15  Temperature  influence  on  peak  pressure 
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Electronic  measurement  chain  for  the  dynamic  calibration  (especially  for 
piezoelectric  pressure  transducers): 

Piezoelectric  pressure  transducers  deliver  electric  charges  at  compressive  load 
(see  fig.  16).  Charge  amplifiers  convert  this  charge  into  voltage,  which  is 
proportional  to  the  charge.  As  the  dynamic  calibration  takes  place  within  approx. 
2-5  msec,  the  voltage  signal  has  to  be  digitized  by  a high-speed  A/D  converter 
(12  bit  resolution,  1 psec  sampling  rate)  for  this  fast  event  and  the  measured 
values  have  to  be  stored  in  the  memory  e.g.  16  kByte  (transient  recorders,  digital 
storage  oscilloscopes). 


Measurement  Chain  for  Dynamic  Calibration 


Pressure 


'tv 


Piezoeleclric 
Pressure  Transducer 


Charge 

Amplilier 


Analog/Digital 

Converter 


MEMORY 


Kv 


Analog  Display 
Pressure  Curve 


Transient  Recorder 


Fig.  16  Principle  of  the  electronic  measurement  chain  for  the  dynamic  pressure 
measurement 


A highly  accurate  voltage  source  U and  a highly  accurate  capacitor  C generate 
the  reference  charge  Q = C*U  for  calibrating  the  measurement  chain. 


Accuracy  of  the  charge  amplifier  and  AD-Converter 
(range  600  - 20000  pC) 

+/-  0.35  % 

Uncertainty  +/-  0.03  % 

Accuracy  of  the  charge  calibrator 

Voltage  source  (B692.A01);  +/-  0.050  % 

or 

Voltage  source  (B201):  +/-  0.100  % 

calibration  capacitor  (1 000  pC)  +/-  0.1 00  % 
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Uncertainty  of  the  electronic  measurement  chain 


+/-  0.23  % 


Total  repeatability  within  a measurement  series 

Four  pressure  transducers  were  installed  into  the  measuring  head  and  50 
repetitions  were  made  for  each  pressure  step.  The  pressure  transducers  were 
not  dismounted  during  the  complete  measurement  procedure. 


Table  2:  Repeatability  within  a series 


Step  a Range  Range  / 


[MPa] 

[MPa] 

[MPa] 

[MPa] 

[%] 

100 

101.2 

0.04 

0.2 

0.19 

300 

301.5 

0.07 

0.3 

0.10 

500 

501.4 

0.20 

0.7 

0.15 

600 

604.5 

0.70 

2.7 

0.43 

800 

792.9 

0.90 

2.6 

0.33 

Total  reproducibility  with  mounting/dismounting  of  pressure  transducers 
(30  repetitions  per  pressure  step,  9 combinations  with  three  different  types  of 
pressure  transducers) 


A reference  transducer  (AVL  5QP6000T)  was  mounted  in  one  bore  of  the 
measuring  head.  It  was  not  dismounted  during  the  measurement.  Similar  types 
of  pressure  transducers  were  randomly  mounted  or  dismounted  in  the  remaining 
three  bores  (see  table  3). 
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Tabte  3: 


Reproducibility  within  a series  by  exchanging  transducers 
(reference  transducer) 


step 

[MPa] 

P 

' mean 

[MPa] 

o 

[MPa] 

Range 

[MPa] 

Range/P, 

[%] 

200 

201.2 

0.7 

2.8 

1.3 

400 

398.5 

0.6  , 

2.3 

0.6 

600 

598.9 

1.1 

3.5 

0.6 

Methods  for  calibrating  the  calibrator 

Mathematical  model  [2] 

A simplified  model  (spring-dashpot  model)  allows  to  derive  the  relation  between 
the  falling  height  and  the  maximum  pressure  (see  figure  17).  The  mathematical 
description  of  this  model  results  in  a nonlinear  differential  equation  of  second 
order  (nonlinear  spring  and  damping  characteristics). 


Fig.  17  Mathematical  model 


M • i'(t)  + D[x(t)]  • x(t)  + C[x(t)]  • x(t)  = 0 

x(t)-C[x(t)] 

A 
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The_experimental!y  derived  spring  and  damping  characteristics  allow  a numerical 
solution  of  the  equation  (see  fig.  18). 


Fig. 


Characteristic  of  Spring  Stiffness 

Force  of  Spring  [kN] 


Characteristic  of  Damping 

Force  of  Damping  [kN] 
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The-  result  is  the  mathematical  dependency  of  falling  height  to  generated 
pressure,  see  table  4 and  fig.  19/20. 


Table  4;  Results 


Step 

H 

P P 

meas  eval 

e ' 

^gale 

T 

A 

T 

B 

[bar] 

[mm] 

[bar] 

[bar] 

[%] 

[ms] 

[m/s] 

[ms] 

[ms] 

1000 

13.58 

949 

991 

0.9 

13.87 

0.505 

1.303 

2.635 

2000 

48.90 

1912 

2004 

-0.2 

6.47 

0.973 

1.250 

2.535 

3000 

102.16 

2903 

2989 

0.4 

4.41 

1.401 

1.214 

2.470 

4000 

171.28 

3911 

4001 

0.0 

3.39 

1.810 

1.180 

2.408 

5000 

254.96 

4927 

4998 

0.1 

2.78 

2.200 

1.152 

2.355 

6000 

352.61 

5956 

6000 

0.0 

2.37 

2.577 

1.128 

2.306 

7000 

464.42 

6994 

7004 

-0.1 

2.08 

2.942 

1.106 

2.261 

8000 

591.30 

8049 

8031 

-0.4 

1.85 

3.302 

1.085 

2.221 

p[MPa] 


0 100  200  300  400  500  600  700 

Falling  Height  [mm] 


Fig.  19  Dependence  pressure  - falling  height 
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p [MPa] 


V ^[m/s] 

Fig.  20  Dependence  pressure  - impact  velocity 


When  using  the  mathematical  model  for  the  calibration  of  the  calibrator  there  are 
uncertainties  at  the  determination  of  the  and  times  of  approx. 

+/-  5 - 10  jisec.  Considering  the  uncertainties  of  the  other  influencing 
magnitudes,  such  as  impact  mass  weight,  piston  area  and  impact  velocity,  the 
uncertainty  of  the  mathematical  model  for  the  pressure  calculation  is 
approx.  +/-  3.5%  . 

Mathematical  representation  of  relation  between  pressure  and  falling 
heights 

When  the  falling  heights  for  the  measured  pressure  values  are  known,  these 
pairs  are  used  for  the  polynomial  calculation  of  fourth  order  according  to  the 
least  squares  method  (see  table  5). 

H{p)  = ag  + ai  -p  + a2‘p2  + ag-p^  + a4*p^ 


Based  on  this  mathematical  calculation  also  other  pressure  values  can  be 
selected.  The  relevant  falling  height  are  automatically  calculated  by  the  B621CS 
Control  System  from  the  stored  polynomial  coefficients  aQ,  ...,  a^  . 
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Table  5:  Coefficients  for  polynomial  of  fourth  order 


Step  H Pe^ai  (H  - 

H 

[bar]  [mm]  [bar]  [mm]  [%] 


1000 

13.58 

995 

13.54 

0.24 

2000 

48.90 

2004 

49.07 

-0.33 

3000 

102.16 

2992 

101.81 

0.32 

4000 

171.28 

4001 

171.60 

-0.18 

5000 

254.96 

4998 

254.91 

0.02 

6000 

352.61 

6000 

352.47 

0.04 

7000 

464.42 

7004 

464.53 

-0.02 

8000 

591.30 

8011 

591.27 

0.00 

ao  = -6913.0721-10^5 
ai  = 1 667.9072-1 0-<56 

a2  = 1281.3057-10-08 
as  = -7798.4908-10-13 
04  = 3835.7080-10-17 


As  the  pair  of  values  (pressure  versus  falling  heights)  does  not  lie  exactly  on  the 
polynomial  curve,  a systematic  error  of  maximally  +/-  0.5  % in  the  height 
adjustment  (for  falling  heights  < 100  mm)  is  caused.  With  regard  to  the  peak 
pressure  this  represents  a relative  error  of  +/-  0.25  % . 
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Transfer  calibration  method  as  a traceable  method 


One  selected  reference  transducer  and  three  test  transducers  are  mounted  on 
the  transfer  measuring  head  into  well  defined  bores  with  accurate  torque  (see 
fig.  21) 


Fig.  21  Set-up  of  transfer  measuring  head 

Classic  static  calibration  cycles  are  applied  and  the  data  retrieved.  The  same 
transducers  on  the  same  measuring  head  are  exposed  to  dynamic  calibration 
pulses,  eliminating  all  minor  error  influences  caused  by  e.g.  bore,  adapter, 
torque  etc.  When  the  results  of  the  pressure  transducers  to  be  tested  and  the 
ones  of  the  reference  transducer  are  compared,  they  should  be  in  accordance 
within  the  uncertainty  of  the  electronic  measuring  chain  (lapprox.+/-  0.2  %)  . 
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However,  in  practice  deviations  of  individual  pressure  transducers  are  observed, 
which  lie  far  beyond  this  value  fig.  22  , see  [3], 


■o  Transducer  1 


ro 

C 

E 

o 

c 

E 

o 

L— 

c 

_o 

fO 

"> 

lU 

Q 


Nominal 

pressure 

kbar 


Fig.  22  a)  Same  pressure  - different  result  b)  Same  pressure  - same  result 
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Irvcase  of  further  calibrations  of  pressure  transducers  e.g.  in  another  measuring 
head  at  least  one  bore  for  the  selected  reference  pressure  transducer  has  to  be 
available. 

This  method  can  also  be  used  to  calibrate  the  calibrator.  In  that  case  static  and 
dynamic  results  have  to  show  a difference  of  less  than  +/-  1 % . The  measured 
pressure  values  of  the  reference  transducer(s)  are  used  for  computing  the 
polynomial  of  4th  order  which  describes  the  dependency  between  pressure  and 
falling  height. 

This  method  has  already  been  applied  at  a customer  in  cooperation  with  the 
national  Bureau  of  Standards.  The  customer  defined  four  specific  pressure 
transducers  as  reference  transducers.  With  these  transducers  and  the  relevant 
transfer  measuring  head  the  static  and  dynamic  calibration  was  carried  out 
without  mounting  or  dismounting  the  transducers.  All  four  reference  pressure 
transducers  show  within  0.8  % the  same  pressure  value  at  each  individual 
pressure  step. 

Prospects 

Improvement  of  the  total  reproducibility  and  accuracy 
larger  foundation  mass 

measure  temperature  of  measuring  head  with  relevant  pressure 
correction 

better  mathematical  representation  of  the  relation  between  pressure  and 
falling  height  (representation  with  polynomial  of  4th  order) 

Introduction  of  an  adapter  measuring  head 

The  measuring  head  is  provided  with  four  bores  for  the  adapters,  which  can  be 
manufactured  for  any  transducer  type  in  such  a way  that  the  measuring  head 
volume  always  remains  constant  when  combining  different  pressure  transducer 
types.  Further  advantage  is  the  equal  adapter  material  for  calibration  (static  and 
dynamic)  and  in  practical  use. 


Introduction  of  a fifth  pressure  transducer  bore 


for  firmly  mounting  the  reference  pressure  transducer  in  the  measuring  head, 
this  transducer  will  not  be  dismounted  and  will  be  calibrated  statically  in  the 
measuring  head,  if  necessary  (traceable  method). 
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Developing  measurement  methods  for  obtaining  the  primary  standard 

Starting  from  the  definition  of  the  adiabatic  compressibility  k (the  pressure  is 
generated  so  quickly  that  no  heat  transfer  takes  place  during  this  time),  the 
pressure  generated  can  be  calculated: 

k(p)  =-1A/-dV/dp  (V  = a*x  + Vo,  dV  = A'dx) 

integration  ; i 

p(V)  = - I 1/k(p)  • A / (a-x  + Vo)  dx 
0 

From  this  equation  it  can  be  seen  that  when  knowing  the  compressibility  of  the 
liquid,  the  volume  of  the  measuring  head  and  the  piston  travel,  the  pressure 
generated  can  be  calculated.  The  difficulty  for  obtaining  such  primary  standard 
lies  in  the  insufficient  accuracy  of  the  presently  known  compressibility  values  (at 
1 00  - 800  MPa)  and  in  the  piston  travel  which  is  not  easy  to  measure. 

Measuring  the  acceleration  of  the  impact  mass  during  the  impact  with 

Bruehl  & Kjaer  piezoelectric  acceleration  transducers 
Polytec  laser  vibration  meter  / interferometer 


Fig.  23  shows  the  experimental  set-up  as  applied  for  measuring  the  course  of 
acceleration,  the  velocity  and  the  travel  of  the  impact  mass  during  the  impact. 
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The  courses  of  velocity  and  travel  measured  by  the  laser  vibration  meter  were 
used  for  calculating  the  acceleration  course  (see  fig.  24).  The  maximum 
acceleration  from  the  acceleration  measurements  was  used  for  further 
calculations  (see  fig.  25). 


V [m  sec""^] 


^ Laser  Measurements 

X [mm] 


Fig.  24  Results  of  laser  measurements 
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Accelerometer  Measurements 


Fig.  25  Results  of  acceleration  measurements 

The  peak  pressure  p^^  was  calculated  by  using  the  formula 

p^^  = M * a^3^  / A M impact  mass 

a acceleration 

A piston  area 

It  corresponds  to  the  peak  pressure  measured  with  the  pressure  transducer  in 
the  range  of  < 4 % . 

Conclusion 

The  error  influences  of  the  dynamic  pressure  calibrator  were  explained  and 
discussed  in  detail.  Furthermore,  the  present  level  of  calibration  methods  for 
calibrating  the  pressure  calibrator  was  outlined  and  prospects  for  further 
development  work  in  order  to  reach  primary  standard  for  dynamic  pressure  were 
given. 

1 would  like  to  thank  Mr.  J.  Winkler  for  his  valuable  ideas,  Mr.  R.  Kristjan  for  his 
technical  assistance  and  Mrs.  U.  Wutti  for  the  translation  of  this  paper  into 
English  and  Mr.  R.  Boschanig  for  proof-reading,  as  well  as  all  members  of  the 
advertising  department  for  their  support. 
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ABSTRACT 

After  a brief  review  of  the  different  methods  permitting  the  determination  of  the  transfer 
function  of  a transducer,  this  paper  presents  the  means  of  dynamic  calibration  of  pressure 
transducers  in  a gaseous  condition,  working  in  the  laboratory  of  dynamic  metrology  at  the  Ecole 
Nationale  Superieure  d’Arts  et  Metiers  of  Paris,  and  more  especially  the  shock  tubes  and  the  fast 
opening  devices. 

The  operation  of  the  aperiodic  pressure  generators  is  described,  and  the  limits  of  use  in 
amplitude  and  frequency  are  given  and  also  some  examples  of  calibration. 

RESUME 

Apres  un  bref  r6sum6  des  diff^rentes  m6thodes  permettant  la  determination  de  la  fonction 
de  transfert  d’un  capteur,  cet  article  pr^sente  les  moyens  d’6talonnage  dynamique  des  capteurs 
de  pression  en  milieu  gazeux,  en  service  au  laboratoire  de  M6trologie  dynarrique  a I’Ecole 
Nationale  Superieure  d’Arts  et  Metiers  de  Paris,  et  plus  specialement  les  tubes  k choc  et  les 
disposititfs  a ouverture  rapide. 

Le  fonctionnement  des  capteurs  de  pression  aperiodique  est  decrit,  et  les  limites 
d’utilisation  en  amplitude  et  en  frequence  sont  donnees  ainsi  que  quelques  examples 
d’etalonnage. 
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1  INTRODUCTION 


In  a static  state,  a pressure  transducer  is  characterized  by  its  sensitivity  which  is  the  ratio  of 
the  variation  of  the  output  over  the  variation  of  the  input.  The  sensitivity  is  practically  constant 
within  the  range  of  the  transducer. 

In  a dynamic  state,  when  the  input  varies  with  time,  keeping  the  above  mentioned  ratio, 
and  comparing  it  with  the  sensitivity,  some  differences  appear. 

A transducer  used  in  a dynamic  state  cannot  be  characterized  only  by  its  sensitivity,  the 
objective  of  dynamic  calibration  is  the  determination  of  the  transfer  function  (gain  and  phase  in 
function  of  the  frequency)  able  to  describe  the  working  of  a pressure  transducer  in  its  static  use 
or  in  its  dynamic  use. 

To  have  access  to  this  transfer  function,  two  important  methods  are  currently  used  : 

- the  harmonic  method 

- the  transient  method. 

To  use  these  two  methods,  two  important  typ>es  of  means  are  used  : 

- periodic  pressure  generators  which  give  a sinusoidal  signal 

- aperiodic  pressure  generators  which  give  a step  or  a pressure  impulse. 


2  METHODS  OF  DYNAMIC  CAUBRATION 

In  the  harmonic  tests,  the  transducer  is  submitted  to  a sinusoidal  input  pressure  of  a given 
amplitude  and  a given  frequency.  The  ratio  of  the  amplitude  of  the  transducer  output  over  the 
transducer  input  gives  the  gain  curve  and  the  phase-angle  between  the  output  and  the  input,  the 
phase  curve  for  a given  frequency.  Through  varying  the  input  frequency  point  by  point,  the 
transfer  function  can  be  determined. 

In  the  transient  tests,  the  transducer  is  submitted  to  a pressure  signal  test  (impulse  - step) 
and  the  response  is  recorded.  Then  the  transfer  function  H(v)  is  calculated  point  by  point  and  is 
defined  as  the  ratio  of  the  Fourier  transform  of  the  transducer  output  S(v)  over  a Fourier  transform 
of  the  input  E(v) 


S(v)  _ s(t)  e dt 

^ e(t)  e dt 


3  MEANS  OF  DYNAMIC  CAUBRATION 

3.1  Periodic  pressure  generators 

There  are  many  types  of  periodic  generators  which,  in  most  cases,  consist  of : 

- a cavity  in  which  the  pressure  changes  sinusoidally  and  where  the  transducer  to  be 
calibrated  is  placed 

- a pressure  modulator  (variation  of  volume  in  the  cavity  or  variation  of  the  mass  inside  the 
cavity) 

- a modulation  motor  with  direct  control  in  frequency  (electrodynamic  or  piezoelectric 
actuator)  or  indirect  control  in  frequency  with  a cinematic  electric  motor  which  converts  the 
circular  motion  into  sinusoidal  linear  motion). 
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Most  periodic  generators  work  in  gaseous  conditions  and  their  range  of  operation  is  in  the 
low  pressures  and  low  frequencies.  In  actual  fact,  the  generated  amplitudes  are  low  because  the 
modulation  factor  must  be  low  in  order  not  to  entail  too  high  a distorsion  rate.  The  frequencies  are 
also  relatively  low  for  the  pressure  In  the  cavity  to  remain  constant  (important  ratio  of  the 
wave  length  over  the  main  size  of  the  cavity). 

Finally,  one  must  bear  in  mind  that  there  is  practically  no  absolute  periodic  generator;  it  is 
always  necessary  to  use  a reference  transducer  or  considered  as  such. 

3.2  Aperiodic  pressure  generators 

There  are  essentially  two  main  types  of  pressure  step  generators  which  work  in  gaseous 
conditions ; 

- the  shock  tube 

- the  fast  opening  device. 

The  shock  tube  is  the  most  frequently  used  means  and  its  application  concerns  the  high 
frequencies. 

The  limitation  in  low  frequency  is  essentially  a function  of  the  duration  of  the  step  pressure. 
The  range  in  amplitude  is  in  low  pressure  as  well  as  high  pressure. 

The  fast  opening  device  is  the  extension  of  the  shock  tube  in  the  low  frequency. 

It  allows  a step  function  of  long  duration  to  be  generated,  but  the  rise  time  is  greater  than 
the  one  met  in  the  shock  tube.  Its  application  range  goes  from  quasi-continuous  to  medium 
frequency  limit  which  depends  a lot  on  the  means. 

The  amplitude  range  is  the  same  as  with  the  shock  tube. 


3.3.  A few  remarks  on  generators 

Periodic  pressure  generators  have  a limited  use,  as  the  frequency  and  pressure  range  are 
very  narrow.  Besides,  to  cover  that  narrow  range,  the  number  of  generators  to  be  designed  is 
excessively  high.  However,  user  requirements  far  exceed  that  application  range. 

Aperiodic  pressure  generators  with  only  two  main  types  of  equipment  cover  the  high  and 
low  pressures  anyway  and  also  the  high  and  low  frequencies.  In  spite  of  the  obvious  difficulties  in 
carrying  out,  exploiting  and  interpreting  the  results  obtained,  the  Laboratoire  of  M6trologie 
Dynamique  of  E.N.S.A.M.  PARIS  is  orientated  towards  the  development  of  aperiodic  means  which 
seem  to  meet  user  requirements  better. 


3.  FAST  OPENING  DEVICE 

3.1  Operation 

A fast  opening  device  consists  of  three  components  ; 

- a small  chamber  at  pressure  PI  where  the  transducer  to  be  calibrated  is  placed. 

- a large  chamber  at  P2 

- a system  which  separates  the  passage  of  two  chambers  and  creates  a fast  opening 
between  the  two. 

The  chamber  at  pressure  P1  being  less  important  than  the  chamber  at  pressure  P2,  it 
results  that  when  the  opening  device  is  actuated,  the  pressure  in  the  small  chamber  changes 
from  P1  to  P2. 

The  transducer  in  the  latter  is  submitted  to  a step  pressure  of  amplitude  P2  - P1 . 
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Because  P2  can  be  greater  or  smaller  than  Pi , the  excitation  can  be  a positive  step  or  a 
negative  step. 


Transducer 
Small  chamber  (Pi) 

Opening  system 


Large  chamber  (P2) 


3.2  Limits  of  use 
- Amplitude  limits 

Concerning  the  high  pressures,  the  fast  opening  devices  are  only  limited  in  step  generation 
by  the  pressure  resistance  of  two  chambers  and  by  the  available  power  for  controlling  of  the 
opening  system. 

Concerning  the  low  pressures,  the  fast  opening  devices  can  generate  lower  steps  than  in 
the  case  of  the  shock  tube  without  too  great  much  difficulty. 

- Frequency  limits 

Concerning  the  low  frequencies,  the  fast  opening  devices  are  only  limited  by  the  duration 
of  the  recording  of  the  transducer  response  to  be  calibrated  at  a pressure  step,  as  the  level  time 
can  be  as  long  as  we  want. 

The  limits  concerning  the  high  frequencies  depend  especially  on  the  opening  system  and 
on  the  small  chamber  volum. 

In  order  to  determine  the  frequency  limit,  it  is  necessary  to  measure  the  “shape"  of  the 
transition  between  pressures  Pi  and  P2  with  a reference  transducer  dynamically  tested  at  the 
shock  tube  and  to  determine  the  transfer  function  of  the  means. 

In  these  conditions,  we  can  effect  a calibration  using  an  absolute  method  (without  a 
reference  transducer)  or  through  a comparison  method  (with  a reference  transducer).  In  the 
latter  case,  the  uncertainties  are  greater  but  the  limits  in  frequency  are  higher. 

Figure  2 gives  a step  response  example  obtained  on  a fast  opening  device  and  figure  3 the 
transfer  function  of  this  one. 


4.  SHOCK  TUBE 
4.1  Operation 

A shock  tube  has  two  chambers  with  a constant  section  separated  by  a membrane.  The 
high  pressure  chamber  (HP)  contains  the  drive  gas,  and.  the  low  pressure  chamber  (LP)  the 
working  gas. 
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Just  before  the  breaking  of  the  membrane,  the  drive  gas  is  in  state  4 and  the  working  gas 
in  state  1 . When  the  membrane  breaks  at  time  t = 0,  waves  represented  in  a distance  (x)  - time  (t) 
diagram  propagate  in  chambers  HP  et  LP  (fig.  4). 


1 1 

Chamber  HP 

1 J 

Chamber  LP 

f ^ 

1 

Fig.  4 ; Shock  tube  Distance  (x)  / Time  (t) 


In  chamber  HP,  a rarefaction  wave  sends  the  gas  from  the  state  4 to  the  state  3 and 
propagates  in  chamber  HP.  In  chamber  LP,  a shock  wave  sends  the  working  gas  from  state  1 to 
state  2 and  propagates  at  the  bottom  of  chamber  LP.  The  shock  wave  compresses  the  gas  it  is 
in  contact  with  and  imparts  it  with  a given  speed.  A contact  wave  is  aiso  impelled  to  the  bottom 
of  chamber  LP  and  separates  the  drive  gas  from  the  working  gas.  The  pressures  and  the 
speeds  are  equal  on  either  side  of  this  wave,  however,  the  temperatures  are  different.  When  the 
first  rarefaction  wave  goes  to  the  bottom  of  chamber  HP,  it  reflects  on  the  solid  wall  in  a 
rarefaction  wave  and  the  pressure  falls  again  behind  it.  When  the  shock  wave  reaches  the  end  of 
chamber  LP,  it  reflects  on  the  solid  wall  as  a shock  wave  and  the  gas  goes  from  state  2 to  state 
5.  This  shock  wave  recompresses  the  gas  which  is  in  front  of  It  and  stays  at  rest. 

A pressure  transducer  placed  on  the  wall  of  chamber  BP  is  initially  at  pressure  PI ; at  the 
passage  of  the  incident  shock,  its  pressure  climbs  to  P2.  For  a transducer  placed  at  the  bottom 
of  chamber  LP,  the  initial  pressure  is  always  PI , but  the  pressure  goes  directly  to  P5  during  the 
reflection  of  the  incident  shock. 

During  a limited  time,  a transducer  placed  on  the  wall  (or  at  the  bottom)  is  submitted,  to  a 
step  of  amplitude  P2  - P1  (or  P5  - P1);  then,  the  transducer  ’sees"  the  continuation  of  the 
phenomenon  (the  reflected  shock  wave  and/or  the  reflected  rarefaction). 

The  figure  5 shows  an  example  of  the  signal  obtained  at  the  bottom  of  the  chamber  LP  for 
a pressure  ratio  about  1,3.  Figure  6 gives,  for  this  same  pressure  ratio,  the  signal  issued  from  a 
transducer  placed  on  the  wall. 
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Rg.  5 : Step  pressure  in  the  end  of  shock  tube  TCR 


Fig.  6 : Step  pressure  in  the  side  wall  of  shock  tube  TCR 
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4.2  Limits  of  use 


- Ampirtude  limits 

The  amplitude  limit  (high  or  low)  of  the  pressure  step  generated  in  a shock  tube  has  at 
inception  : 

- a purely  theoretical  limit 

- a limit  due  to  disregard  of  the  hypothesis  corxjitions  (real  operation) 

- a technological  limit. 

The  theoretical  limit  is  not  a problem  because  in  all  cases  the  real  operation  of  the  shock 
tube  sets  in  before  it  is  reached.  As  a matter  of  fact,  in  order  that  the  flow  in  the  shock  tube  be 
kept  sufficiently  close  to  the  theoretical  flow  (equation  of  state  in  real  gas  close  to  the  equation  of 
state  in  perfect  gas),  the  pressure  ratio  P4/P1  must  be  as  close  to  the  unit  as  possible. 
Nevertheless,  in  order  to  obtain  high  amplitudes,  it  is  necessary  to  increase  the  pressure  in 
chamber  LP. 

With  regard  to  technology,  obtaining  high  amplitudes  for  creating  high  pressures  is  difficult 
and  the  safety  measures  to  be  taken  make  its  use  tricky. 

Concerning  the  low  pressures,  the  difficulties  come  from  the  membrane.  Its  breakage  is 
often  unsatisfactory  and  ctenges  the  quality  of  the  pressure  step  a great  deal. 

- Frequency  limits 

The  frequency  limits  in  the  dynamic  calibration  in  shock  tubes  are  two-fold  : 

- the  limit  in  high  frequency,  which  is  function  of  the  rise  time  of  the  pressure  step  effectively 
generated 

- the  limit  in  low  frequency  which  depends  essentially  on  the  duration  of  the  level  of  the 
pressure  step. 

The  position  of  the  calibration  section  (wall  or  bottom  of  the  tube  BP)  is  to  be  taken  into 
account  as  well  as  the  shape  and  the  physical  dimensions  of  the  sensitive  surface  of  the  tested 
transducer. 

- High  frequency  limits 

Theoretically,  the  "thickness*  of  a shock  wave  (the  necessary  distance  for  the  pressure  to 
go  from  P1  to  P2  or  from  Pi  to  P5)  is  very  small  (some  microns)  and  moreover,  the  speed  of  the 
shock  wave  is  great,  the  rise  time  is  generally  assumed  to  be  below  the  microsecond  for  an 
imaginary  observer  placed  at  a definite  point  on  the  wall  or  at  the  bottom  of  the  shock  tube. 

In  the  case  of  a calibration  at  the  wall  of  the  shock  tube,  the  actual  rise  time  depends  on  the 
speed  of  the  incident  shock  wave  and  on  the  diameter  of  the  transducer. 

This  limit  is  the  higher  as  the  speed  of  the  incident  shock  is  important  and  the  diameter  of 
the  transducer  is  small. 

At  the  bottom  of  the  shock  tube,  the  rise  time  is  affected  nevertheless  by  the  curve  of  the 
incident  shock  wave  and  by  the  possible  positioning  defects  of  the  transducer. 

These  two  phenomena  contribute  to  increasing  the  effective  rise  time  of  the  pressure  step 
but  in  proportions  less  important  than  at  the  wall  of  the  shock  tube. 


Low  frequency  limit 

The  low  frequency  limit  in  a shock  tube  is  due  to  the  finite  duration  of  the  pressure  level,  for 
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a calibration  at  the  wall  or  at  the  bottom  of  the  tube. 


The  duration  of  the  level,  and  consequently  the  low  frequency  limit,  is  essentially  a function 
of  the  length  of  a shock  tube. 

Figure  7 gives  a step  response  example  of  a transducer  obtained  on  shock  tube  and  figure 
8 the  transfer  function  of  this  one. 

5 MEANS  OF  ENSAM  PARIS 


During  the  last  few  years,  ENSAM  PARIS  has  developed  with  the  help  of  the  Bureau 
National  de  M6trologie  (BNM)  several  means  in  order  to  satisfy  the  pressure  frequency  domain 
needed  by  the  users. 

Shock  tubes 


The  shock  tubes  developed  by  the  laboratory  are  ; 

- shock  tube  TC10  which  chambers  HP  and  LP  are  made  to  work  at  10  bars  and  which 
enables  to  generate  pressure  steps  of  amplitude  from  0,1  to  6 bars. 

- shock  tube  TCI  00  is  able  to  work  in  the  range  of  pressure  from  5 to  80  bars. 

- shock  tube  TC200  is  able  to  work  with  a static  pressure  up  to  200  bars. 

- shock  tube  TCR  works  in  low  pressure  (1  to  4 bars)  but  its  step  duration  is  about  50  ms. 


Fast  opening  device 

The  fast  opening  devices  developed  at  the  laboratory  are  ; 

- fast  opening  device  DOR  10  works  with  a static  pressure  up  to  10  bars  and  is  able  to 
generate  step  pressure  with  a rise  time  of  about  4 ms. 

- fast  opening  device  DOR  100  works  at  100  bars  and  the  rise  time  of  the  step  pressure  is 
about  0.4  ms. 

- fast  opening  device  DOR200  have  the  same  performance  as  DOR100  but  the  static 
pressure  is  up  to  200  bars. 


6 CONCLUSION 


The  laboratoire  de  M6trologie  Dynamique  of  E.N.S.A.M.  PARIS  has  developed  a coherent 
equipment  to  calibrate  dynamic  pressure  transducers  in  gaseous  condition.  The  whole  means 
are  regrouped  in  the  pressure  frequency  diagram  given  figure  9. 

With  shock  tubes  and  fast  opening  devices,  the  laboratory  can  calibrate  sensors  in  pressure 
range  from  a few  mbars  to  a few  hundred  bars  and  in  frequency  range  from  a few  Hz  to  a few 
tens  kHz.  In  addition,  all  the  means  can  be  intercompared  with  the  reference  shock  tube  TCR. 
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Fig.  9 : Frequency  / Pressure  Domain  of  means  for  dynamic  calibration 
of  pressure  transducer  at  ENSAM  Paris 
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Diatomic  gas  molecules  have  a funda- 
mental vibrational  motion  whose  fre- 
quency is  affected  by  pressure  in  a 
simple  way.  In  addition,  these 
molecules  have  well  defined  rotational 
energy  levels  whose  populations  provide 
a reliable  measure  of  the  thermodynamic 
temperature.  Since  information  con- 
cerning the  frequency  of  vibration  and 
the  relative  populations  can  be  deter- 
mined by  laser  spectroscopy,  the  gas 
molecules  themselves  can  serve  as  sen- 
sors of  pressure  and  temperature. 
Through  measurements  under  static 
conditions,  the  pressure  and  temperature 
dependence  of  the  spectra  of  selected 
molecules  is  now  understood.  As  the 
time  required  for  the  spectroscopic  mea- 
surement can  be  reduced  to  nanosec- 
onds, the  diatomic  gas  molecule  is  an 


excellent  candidate  for  a dynamic  pres- 
sure/temperature primary  standard. 

The  temporal  response  in  this  case  will 
be  limited  by  the  equilibration  time  for 
the  molecules  to  respond  to  changes  in 
local  thermodynamic  variables.  Prelimi- 
nary feasibility  studies  suggest  that  by 
using  coherent  anti-Stokes  Raman  spec- 
troscopy we  will  be  able  to  measure  dy- 
namic pressure  up  to  10*  Pa  and 
dynamic  temperature  up  to  1500  K with 
an  uncertainty  of  5%. 
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1.  Introduction 

With  modern  laser  diagnostic  techniques,  it  is 
possible  to  characterize  the  pressure  {P)  and  tem- 
perature (T)  of  a gas  at  the  molecular  level.  The 
measurement  times  for  these  techniques  are  such 
that  the  response  to  changes  in  T and  P is  limited 
only  by  the  fundamental  relaxation  and  transport 
processes  of  the  molecular  system.  This  provides 
the  basis  for  a new  approach  to  the  calibration  of 
transducers  used  in  the  measurement  of  dynamical 
P and  T.  The  essence  of  dynamic  calibrations  is  the 
determination  of  the  time  dependent  response  of 
the  transducer,  which  requires,  at  a minimum,  the 
application  of  a stimulus  with  known  time  depen- 
dence, i.e.,  a “standard”  dynamic  source. 


If  one  were  to  rely  on  conventional  sensors 
(whose  response  functions  are  not  a priori  known) 
to  characterize  the  dynamic  source,  an  inescapable 
circularity  emerges  from  the  preceding  paragraph. 
Approaches  to  solution  of  this  problem  have  tradi- 
tionally [1]  relied  on  some  form  of  calculable 
source.  In  essence,  this  is  a source  some  properties 
of  which  can  be  determined  from  accurate  mea- 
surements, for  example  of  quasi-static  values  of  P 
and  T and  time  rate  of  change  of  position,  and 
whose  time  dependent  P and  T is  then  derived  from 
an  appropriate  theoretical  prescription,  e.g.,  from 
hydrodynamics  for  sound  propagation  or  fluid  me- 
chanics for  shock  waves  (with  appropriate  equa- 
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tions  of  state  for  isentropic  or  adiabatic  expan- 
sions). It  must  be  recognized  that  every  theory  re- 
lies to  some  degree  on  idealizations  and  that  any 
laboratory  realization  of  a dynamic  source  is  non- 
ideal. Thus,  sources  of  dynamical  P and  T cannot 
be  accurately  known  from  theory  alone;  measure- 
ment of  the  “standard”  source  always  is  required. 

Ideally,  in  the  maintenance  of  national  standards 
one  seeks  to  relate  the  measured  quantity  to  a con- 
stant of  nature,  maintained,  for  example,  in  the  en- 
ergy levels  of  an  isolated  atom  or  molecule.  We  are 
proposing  this  type  of  approach  for  the  develop- 
ment of  a “standard”  source  for  dynamic  P and  T. 
The  essence  of  our  approach  is  to  combine  the  very 
best  in  calculable  generators,  fast  transducers,  and 
high-speed  digital  data  acquisition  systems  with  a 
new,  fundamental  measurement  approach.  The  lat- 
ter relies  on  the  use  of  laser-based  diagnostic  tech- 
niques, developed  over  the  past  10  years,  to 
determine  the  P and  T of  the  dynamic  system.  The 
unique  characteristics  of  the  optical  techniques  are: 

•T  and  P are  derived  from  measurement  of  the 
optical  transitions  between  the  atomic  or 
molecular  energy  levels  of  the  constituents  of 
the  dynamic  source,  i.e.,  the  atoms  or  molecules 
are  the  fundamental  transducers  of  the  local  P 
and  T environment 

•optical  measurements  can  be  accomplished 
with  a single  laser  pulse  of  nanosecond  dura- 
tion, with  the  consequence  that  the  “response 
time  of  the  transducer”  reduces  to  the  equilibra- 
tion time  for  the  atoms  or  molecules  (in  the  in- 
teraction region)  to  respond  to  changes  in  the 
local  thermodynamic  variables 

•optical  measurements  can  be  accomplished 
within  harsh  environments  by  means  of  trans- 
mitted or  reflected  laser  beams  and,  for  multiple 
beam  techniques,  spatial  resolution  within  the 
source  volume  can  be  defined  by  the  regions  of 
overlap  of  these  beams,  e.g.,  mm^  dimensions. 

Our  approach  relies  mainly  on  the  use  of  nonlin- 
ear Raman  spectroscopies,  since  these  have  consis- 
tently been  shown  to  provide  useful  diagnostic 
spectra  in  very  short  times  with  high  spatial  resolu- 
tion[2].  The  spectrum  determined  with  these  non- 
linear Raman  approaches  is  the  simplest,  best 
understood,  and  most  highly  characterized  of  any 
optical  diagnostic  technique.  Comparisons  of  spec- 
tra observed  for  systems  in  known  (static)  states  of 


P and  T with  the  predictions  of  theory  provide  a 
high  degree  of  certainty  in  the  use  of  these  data  for 
P and  T measurement. 

The  purpose  of  this  paper  is  to  describe  the  non- 
linear Raman  optical  measurements  that  can 
provide  the  new  primary  standard  for  dynamical  P 
and  T.  This  description  will  include  information  on 
the  P and  T dependence  of  the  spectrum  and  a brief 
consideration  of  the  important  elements  of  a mea- 
surement system  which  can  be  applied  to  a dy- 
namic source.  For  the  purposes  of  this  discussion 
we  do  not  consider  the  dynamic  source  in  any  de- 
tail; however,  the  information  we  present  is  consid- 
ered applicable  to  a suitably  designed  shock  tube 
source.  The  state-of-the-art  in  optical  diagnostics  is 
now  at  a point  where  accurate  measurement  of 
such  a “standard”  dynamic  source  is  possible.  Ac- 
curacy limits  of  the  order  of  5%  for  the  metrologi- 
cally  significant  range  of  P up  to  10*  Pa  and  T up  to 
1500  K appear  achievable. 

In  the  following,  we  begin  with  an  operational 
description  of  the  use  of  nonlinear  Raman  spec- 
troscopy for  P and  T measurement,  drawing  from 
the  already  established  data  base  on  the  T and  P 
dependence  of  observed  spectra.  We  will  then  out- 
line the  elements  of  a measurement  system  for  a 
dynamic  source.  Areas  needing  significant  instru- 
mental development  are  included  in  this  discussion. 
Some  questions  with  regard  to  the  P and  T depen- 
dence of  the  Raman  spectrum  which  need  further 
fundamental  research  also  are  highlighted.  The 
presentation  style  is  intended  to  be  descriptive 
rather  than  rigorous;  for  completeness,  more  de- 
tailed information  on  the  T and  P dependence  of 
nonlinear  Raman  spectra  is  included  in  the  Ap- 
pendix. 

2.  Nonlinear  Raman  Optical  Diagnostics 

The  proposed  approach  to  measurement  of  the  T 
and  P of  the  dynamic  source  is  coherent  anti-Stokes 
Raman  spectroscopy  (CARS)[3].  In  its  most  simple 
realizations  this  technique  uses  two  lasers,  termed 
the  pump  and  the  Stokes  beams,  whose  frequency 
difference  is  selected  to  be  in  resonance  with  a pure 
vibrational  transition  of  a diatomic  gas  molecule 
contained  in  the  source  medium.  The  nonlinear  in- 
teraction of  the  electric  fields  of  these  lasers  with 
the  molecules  of  the  source  medium  generates  a 
third,  laser-like  beam,  termed  the  anti-Stokes  beam, 
which  carries  the  information  about  the  molecular 
system,  in  particular  about  its  T and  P. 
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A useful  arrangement  of  these  beams,  which  pro- 
vides a high  degree  of  spatial  resolution,  is  shown 
in  figure  1.  In  this  configuration,  termed  BOX- 
CARS, the  interaction  volume  is  defined  by  the 
region  of  overlap  of  the  two  pump  beams, 
(derived  from  one  laser  source),  and  the  Stokes 
beam,  k^.  Sample  volumes  of  millimeter  and  submil- 
limeter dimension  are  remotely  accessible  in  this 
arrangement. 

The  information  about  the  local  T and  P envi- 
ronment of  the  molecules  in  the  interaction  volume 
is  determined  from  the  spectrum  of  the  generated 
anti-Stokes  beam  (designated  by  k^^  in  fig.  1).  The 
spectrum  is  obtained  by  measuring  the  power  of 
the  anti-Stokes  beam  as  a function  of  the  frequency 
difference  between  the  pump  and  Stokes  lasers. 
Considering  for  the  moment  a static  system,  a spec- 
trum can  be  obtained  by  measuring  this  power  as 
we  change  the  frequency  difference  between  a tun- 
able narrowband  Stokes-laser  and  a fixed  fre- 
quency narrowband  pump-laser.  Since  this  is 
simply  a power  measurement,  we  can  essentially 
eliminate  the  use  of  traditional  spectroscopic  in- 


Figure  1.  The  approximate  geometrical  arrangement  of  the 
pump,  (subscript  0),  Stokes  (subscript  s),  and  generated  anti- 
Stokes  (subscript  as)  beams  in  a CARS  experiment.  The  sample 
region  is  at  the  intersection  of  the  crossing  beams.  The  phase 
matching  condition  for  (folded)  BOXCARS  [11]  is  indicated. 


struments  (e.g.,  prism  or  grating  spectrometers) 
and  retrieve  an  undistorted  measure  of  the  informa- 
tion imparted  by  the  molecular  system.  The  nar- 
rowband Stokes  and  pump  laser  sources  can 
readily  be  made  of  essentially  delta-function-like 
bandwidth  for  this  application. 

The  special  conditions  for  static  systems  and  nar- 
rowband lasers,  described  in  the  last  paragraph, 
have  been  achieved  in  the  laboratory  in  order  to 
determine  the  fundamental  molecular  response,  i.e., 
its  spectrum,  under  known  conditions  of  T and  P. 
This  has  been  accomplished  for  certain  ranges  of 
these  variables  and  for  a few  selected  molecular 
systems  [4].  We  will  illustrate  the  basics  of  spectro- 
scopic T and  P determinations  by  describing  spec- 
tra derived  from  these  studies. 


3.  Temperature  Dependence  of  CARS 
Spectra 

In  figure  2 we  show  CARS  spectra  of  pure  N2  as 
a function  of  T with  the  pressure  held  fixed  at  1.0 
atm.  The  horizontal  axis  is  the  frequency  difference 
between  the  pump  and  Stokes  lasers  and  the  verti- 
cal axis  is  the  (calculated)  power  in  the  anti-Stokes 
beam  (in  an  unspecified  arbitrary  unit).  These  spec- 
tra are  referred  to  as  vibrational  C 'branch  spectra, 
because  the  optical  transition  involves  a change 
only  in  the  vibrational  quantum  state  (quantum 
number  v)  and  no  change  in  rotational  state  (i.e.,  no 
change  in  the  rotational  quantum  number,  J)  [5]. 
The  relevant  states  and  modes  of  motion  are  sche- 
matically illustrated  in  figure  3. 

Returning  to  figure  2,  we  see  that  there  are  many 
maxima  in  the  power  as  a function  of  frequency 
difference.  Each  of  these  arises  from  a pure  vibra- 
tional transition  which  originates  in  a different  ro- 
tational state  y.  The  vibrational  frequency  depends, 
to  a small  degree,  on  the  rotational  state  because 
the  rotation  of  the  molecule  results  m a slight 
stretching  of  the  bond  length  producing  a small 
change  in  the  forces  binding  the  molecule  and  a 
concomitant  change  in  the  vibrational  frequency. 
The  vibrational  frequency  (we  use  the  traditional 
spectroscopic  unit  cm“‘,  1 cm“'^30  GHz,  in  the 
figure)  is  approximately  2329.91  cm~‘  for  the  J — Q 
rotationless  state;  the  value  decreases  approxi- 
mately according  to  0.017387(7+1),  which  places 
the  7=10  transition  at  :^2328.00  cm“',  less  than  a 
0.1%  change.  The  strength  (integrated  area)  of  an 
individual  transition  is  a function  of  the  population 
difference  between  the  initial  (v  =0,  7)  and  the  final 
(v  = l,  7)  states  of  the  transition.  This  population 
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Figure  2.  Calculated  CARS  spectra  for  the  Ni  vibrational  Q 
branch  as  functions  of  T for  fixed  P {=\  atm).  The  horizontal 
axis  is  the  frequency  difference  between  the  pump  and  Stokes 
lasers.  The  vertical  axis  is  a measure  of  the  CARS  power. 
Although  the  absolute  units  of  this  power  are  arbitrary,  the  rela- 
tive magnitudes  as  a function  of  T are  accurately  represented. 
Selected  transitions  and  bands  of  the  complete  spectrum  are  in- 
dicated. In  the  bottom  panel,  the  spectral  region  from  0(16) 
thru  0(0)  is  shown  on  an  expanded  frequency  scale. 

dependence  in  the  relative  strengths  of  the  transi- 
tions as  a function  of  the  rotational  level,  J,  is  the 
basis  for  temperature  determination,  since  for  sys- 
tems in  thermodynamic  equilibrium  the  state  popu- 
lations are  functions  only  of  the  temperature.  The 
vibrational  ^-branch  spectrum  is  very  useful  for 
measuring  T because  there  are  essentially  no  cor- 
rections to  apply  in  order  to  relate  the  strength  of  a 
fully  resolved  transition  to  the  population  differene 
and  therefore  to  the  T [6].  As  is  seen  in  figure  2,  the 


^ vibrate 
rotate 


R(N-N) 


Figure  3.  The  solid  curve  schematically  represents  the  potential 
energy  of  the  ground  electronic  state  versus  internuclear  separa- 
tion. Vibrational  and  rotational  energy  levels  (quantum  numbers 
V and  J,  respectively)  also  are  indicated  (not  to  scale).  A molec- 
ular transition,  0(2),  associated  with  the  vibrational  0 branch  is 
indicated  by  the  arrow. 

transitions  for  higher-/  states  increase  in  strength 
with  increasing  T;  this  simply  mimics  the  popula- 
tion shifts  to  higher  energy  states  associated  with 
increasing  T. 

Strictly  speaking,  the  temperature  determined 
from  the  relative  populations  of  the  rotational  lev- 
els should  be  called  a “rotational  temperature.”  In 
like  manner,  a “vibrational  temperature”  can  be  de- 
termined from  a measurement  of  the  relative  popu- 
lations of  the  vibrational  levels.  We  observe  in  the 
higher- T spectra  in  figure  2 that  there  are  spectral 
maxima  for  transitions  labeled  v = l— >v  = 2.  This  is 
a vibrational  Q branch  which  originates  in  the  first 
excited  vibrational  level,  v = l,  and  terminates  in 
the  second  excited  level,  v = 2.  At  sufficiently  high 
T there  is  a significant  population  in  the  v = 1 state 
and  this  transition  becomes  observable.  This  transi- 
tion is  totally  analogous  to  that  discussed  above 
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which  initiated  in  the  vibrational  ground  state, 
v==0.  Comparison  of  the  integrated  areas  between 
the  1-^2  and  0-^1  transitions  gives  a measure  of  the 
“vibrational  temperature.”  The  rotational  tempera- 
ture of  the  vibrationally  excited  state  also  can  be 
determined.  The  assumption  of  local  thermody- 
namic equilibrium  can  thus  be  tested,  since  all 
levels  should  yield  the  same  thermodynamic  tem- 
perature in  the  equilibrium  situation. 

4.  Pressure  Dependence  of  CARS  Spectra 

First  we  note  that  the  integrated  area  of  the  en- 
tire spectrum  is  a function  of  the  number  of 
molecules  with  which  the  intersecting  laser  beams 
interact.  This  feature  often  is  used  as  a means  of 
measuring  species  concentration  in  a diagnostic  en- 
vironment [2,3].  Each  molecular  species,  e.g.,  N2, 
O2,  CO,  H2,  etc.,  has  a separate  vibrational  reso- 
nance because  the  resonant  frequency  is  a sensitive 
function  of  the  binding  forces  and  the  masses  of  the 
atoms  comprising  the  molecule.  Thus,  intercom- 
parison of  the  areas  of  these  different  Q branches 
can  be  used  as  a measure  of  the  relative  numbers  of 
each  molecule  in  the  sample  volume.  Because  it  is 
very  difficult  to  make  accurate  measurements  of 
absolute  intensity,  it  has  been  found  that  the  abso- 
lute intensity  of  a transition  is  not  a good  measure 
of  the  density  or  pressure  of  a sample. 

Fortunately,  there  are  good  measures  of  the 
pressure  of  a sample  which  can  be  recovered  from 
the  0-branch  spectrum  of  some  diatomics.  We  il- 
lustrate these  by  spectra  of  the  Q branch  of  D2, 
which  has  been  extensively  studied  in  our  labora- 
tory. Figure  4 presents  CARS  spectra  calculated 
from  results  of  these  experimental  studies.  The  fea- 
tures of  primary  interest  in  this  figure  are  the  reso- 
nance frequencies  and  widths  of  the  transitions. 
Figure  4 illustrates  the  important  fact  that  these 
transitions  change  their  width  (broaden)  and  their 
resonance  frequency  (shift  to  lower  values)  with 
increasing  pressure.  In  the  pressure  range  shown, 
this  broadening  and  shifting  are  linear  with  pres- 
sure. For  the  0(1)  transition  of  D2,  the  broadening 
rate  is  0.0012  cm“Vatm  (we  use  the  conventional 
half  width  at  half  peak  height  as  our  measure  of 
width)  whereas  the  shift  rate  is  —0.0019  cm  “'/atm, 
i.e.,  the  line  shifts  more  rapidly  than  it  broadens. 
We  see  also  that  the  transitions  in  the  D2  0-branch 
spectrum  remain  isolated,  non-overlapped,  up  to 
100  atm.  We  thus  identify  a very  appealing  ap- 
proach to  pressure  measurement  in  that  it  is  tied  to 


the  measurement  of  the  frequency  positions  and 
widths  of  molecular  transitions. 

The  situation  described  above  for  the  D2  0 
branch  should  be  contrasted  with  the  observations 
of  pressure  broadening  and  shifting  in  the  0 branch 
of  N2.  The  broadening  rates  of  the  N2  transitions 
are  typically  30-40  times  larger  than  that  of  the  D2 
0(1)  transition.  A shifting  rate  for  a N2  transition  is 
typically  10  times  smaller  than  its  broadening  rate, 
thereby  making  accurate  frequency  determination 
less  certain  with  increasing  pressure.  Additionally, 
we  note  that  the  interline  spacing  in  the  N2  0 
branch  is  approximately  l/60th  that  of  D2.  As  a 
consequence,  a transition  such  as  0(10)  would 
broaden  to  overlap  most  of  the  other  transitions 
shown  in  figure  2 at  a pressure  of  100  atm  at  room 
temperature.  This  overlap  of  transitions  leads  to 
important  changes  in  the  appearance  of  the  spec- 
trum which  are  discussed  briefly  in  the  Appendix. 
At  this  point  it  suffices  to  say  that  the  N2  0 branch 
involves  a complicated  spectral  distribution  func- 
tion which  has  the  consequence  that  pressure  de- 
terminations would  generally  be  less  reliable  than 
those  derived  from  the  approach  based  on  the  0 
branch  of  D2. 


Figure  4.  Calculated  CARS  spectra  for  the  7=0  and  7=  1 tran- 
sitions of  the  vibrational  Q branch  of  pure  Di  at  T =295  K.  The 
zero  of  intensity  for  each  P is  shifted  by  an  amount  proportional 
to  P.  The  dotted  vertical  lines  indicate  the  resonance  frequency 
of  the  0(1)  line  at  each  P.  The  solid  sloping  line  drawn  through 
the  vertical  lines  is  thus  an  indication  of  the  linear  with  R shift  of 
the  resonance  frequency.  At  P = 0 the  0(0)  line  is  at  2993.57 
cm“‘  with  the  higher  7 transitions  at  lower  frequencies  approxi- 
mately by  the  amount  1.056J(7-t- 1). 
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To  this  point  we  have  restricted  consideration  to 
static  systems  and  very  high  resolution  measure- 
ments with  narrowband  lasers.  We  have  displayed 
the  results  of  careful  measurements  under  known, 
static  conditions  of  T and  P which  demonstrate  that 
temperature  can  be  determined  by  the  measure- 
ment of  the  relative  intensities  of  molecular  transi- 
tions and  that  pressure  can  be  determined  from  the 
positions  and  widths  of  molecular  transitions.  We 
turn  now  to  discuss  briefly  the  measurement  ap- 
proach required  to  characterize  a dynamic  system. 

5.  Single  Shot  Diagnostic  Measurements 

First  we  point  out  that  we  are  interested  in  very 
rapid  measurements,  in  the  range  of  microseconds 
or  less,  in  order  to  characterize  the  full  bandwidth 
of  transducers.  Next  we  note,  from  the  previous 
discussion,  that  a relatively  large  amount  of  spec- 
tral information  must  be  acquired  in  order  to  char- 
acterize T and  P.  These  requirements  of  time  and 
spectral  range  (at  relatively  high  resolution)  result 
in  the  selection  of  a technique  in  which  all  the  spec- 
tral elements  are  measured  at  once,  i.e.,  a single 
shot  measurement.  Single  shot  measurements  are 
the  strong  suit  of  nonlinear  Raman  diagnostics  be- 
cause we  can  employ  pulsed  laser  sources  and  ob- 
tain high  optical  power  levels  which,  by  virtue  of 
the  inherent  nonlinear  response  of  these  techniques, 
results  in  the  generation  of  strong,  readily  detected, 
CARS  signals. 

The  most  straightforward  approach  to  single 
shot  measurement  is  to  replace  the  scanned,  nar- 
rowband Stokes  laser,  considered  above,  with  a 
broadband  Stokes  laser  whose  frequency  width  is 
sufficient  to  provide  the  entire  Raman  spectrum  of 
interest  in  a single  shot.  For  example,  the  entire  60 
cm*'  region  in  figure  2 would  be  measured.  We 
note  that  a similar  spectral  bandwidth  would  be 
required  for  temperature  measurement  via  the  D2 
Q branch  because  of  the  larger  interline  spacing  in 
the  D2  Q branch,  cf  figure  4.  The  measurement  of 
the  spectrum  then  requires  multiple  detectors  cou- 
pled to  an  instrument  which  disperses  different  fre- 
quency components  in  the  anti-Stokes  beam  to  the 
separate  spatial  locations  of  the  detectors.  This 
combination  of  a spectrometer  and  detector  we 
refer  to  as  a multichannel  system. 

6.  Instrument  Development 

The  basic  tools  for  these  single  shot  measure- 
ments do  exist.  The  necessary  laser  systems  have 


been  realized  in  the  research  laboratory.  Adapta- 
tion of  these  to  specific  diagnostic  situations  will  be 
required.  Similarly,  a large  variety  of  spectrome- 
ter/detector systems  are  available  and  this  should 
allow  for  reasonably  straightforward  realization  of 
the  specific  systems  required  for  our  application. 
Because  of  the  requirements  for  large  spectral 
ranges  and  high  resolution,  noted  above,  some  mul- 
tiplexing of  multichannel  systems  will  be  necessary. 
More  details  on  the  instrumentation  are  discussed 
in  the  following. 

6.1  Lasers 

The  basic  requirements  for  the  laser  sources  are 
well  within  the  state-of-the-art  of  laser  technology. 
The  combination  of  a single-frequency,  pulsed 
Nd:YAG  oscillator  with  a series  of  amplifiers  and  a 
frequency  doubling  crystal  will  provide  an  ade- 
quate Raman  pump  laser.  This  laser  also  will  serve 
as  a pumping  source  for  amplifying  the  required 
broadband  dye  laser  for  the  Stokes  beam.  The  os- 
cillator for  this  Stokes  laser  will  require  special  de- 
sign, tailored  specifically  for  the  transitions 
selected  for  optical  characterization  of  the  dynamic 
source. 

6.2  Spectrometer/Detector  System 

This  system  also  will  be  tailored  to  the  transi- 
tions chosen.  It  will  be  unique  in  that  it  will  have  a 
high  resolution  (0.001-0.05  cm~’)  system  for  fre- 
quency and  linewidth  determination  along  with  the 
more  conventional  low  resolution  (0. 1-0.5  cm“’) 
channel  for  temperature  determinations.  A high 
resolution  multichannel  system  has  not  yet  been  at- 
tempted in  diagnostic  applications.  Our  ability  to 
selectively  add  a component  such  as  D2  to  serve  as 
a pressure  transducer  in  the  working  medium  al- 
lows us  to  consider  this  approach.  Incorporating  a 
significant  fraction  of  N2  in  this  medium  allows  for 
a number  of  independent  checks  on  the  tempera- 
ture and  on  the  state  of  thermodynamic  equi- 
librium. Depending  on  the  pressure  range,  we  may 
be  able  to  determine  the  temperature  from  the  D2 
spectrum.  An  approach  under  consideration  is  to 
simultaneously  measure  the  Q branch  of  D2  for 
pressure  determination  and  the  pure  rotation  S 
branch  (transitions  within  a vibrational  state  in- 
volving a change  of  +2  in  the  rotational  quantum 
number  only)  of  N2  to  determine  the  temperature 

[7]. 

The  discussion  in  the  last  few  paragraphs  indi- 
cates one  area  of  the  research  necessary  to  develop 
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a standard  source  of  dynamic  P and  T.  There  are 
additional  fundamental  studies  required  to  assure 
accurate  knowledge  of  the  T and  P dependence  of 
CARS  spectra  and  ultimately  to  enable  the  selec- 
tion of  a measurement  approach  optimized  for  ap- 
plication to  the  dynamic  source.  We  turn  now  to  a 
brief  discussion  of  these. 

7.  Fundamental  Studies 

The  most  important  fundamental  question  to  be 
answered  in  this  work  is  the  dependence  of  the 
spectra  of  our  specific  dynamic  system  on  T and  P. 
What  is  unique  and  new  about  the  system  we  pro- 
pose is  that  it  is  by  necessity  a mixed  gas  system, 
for  example  D2  contained  at  low  concentrations  in 
Nj.  Except  for  our  recent  studies  (see  ref.  [4]  and 
below),  there  is  relatively  little  work  on  mixed  gas 
systems.  Additionally,  we  are  interested  in  these 
systems  over  quite  large  ranges  of  P and  T. 

There  are  a number  of  motivations  for  the  selec- 
tion of  a mixed  gas  system.  First,  no  one  molecule 
presents  all  the  necessary  spectral  characteristics 
which  allow  an  adequate  characterization  of  both 
T and  P.  Of  the  many  single  component  systems 
considered,  D2  comes  closest  to  fulfilling  this  re- 
quirement. D2  provides  a good  measure  of  P be- 
cause it  has  very  narrow,  isolated  lines.  This 
characteristic  makes  it  difficult  to  simultaneously 
measure  the  relative  strengths  of  all  the  0 -branch 
transitions  to  obtain  an  accurate  measure  of  T.  Fur- 
thermore, it  is  much  harder  (i.e.,  it  takes  much 
longer)  to  establish  thermodynamic  equilibrium  for 
D2  than  it  is  for  N2  because  of  the  much  larger  gaps 
between  the  energy  levels  in  D2.  Large  departures 
from  thermodynamic  equilibrium  can  degrade  the 
performance  of  a dynamic  source  and  the  accuracy 
of  the  spectroscopic  measurements  themselves. 
The  use  of  pure  D2  has  a number  of  other  limita- 
tions; for  example,  a shock  tube  dynamic  source 
would  be  very  limited  in  dynamic  range  by  the  use 
of  such  a light  gas.  We  note  also  that  the  cost  of 
running  a pure  D2  system  would  be  prohibitive,  as 
also  might  the  potential  safety  hazards  associated 
with  pure  D2  at  high  pressures  or  high  tempera- 
tures. 

The  use  of  a mixed  gas  system  is  potentially  an 
advantage  for  pressure  measurements.  This  state- 
ment is  derived  from  our  observation  that  mixed 
gas  systems  such  as  D2:X,  with  X— He,  Ar,  and  N2, 
have  a larger  shift  to  width  ratio  than  pure  D2. 
There  remain  some  fundamental  questions  in  this 
regard,  however,  because  we  have  observed  line 


shape  asymmetries  in  the  spectrum  of  D2:Ar  and 
H2:Ar  [8].  We  also  know  that  these  asymmetries 
are  functions  of  T.  These  observations  need  to  be 
fully  understood  in  order  to  reliably  calibrate  P\ 
furthermore,  these  studies  must  be  extended  to  the 
D2:N2  system.  An  adequate  spectral  model  and  the 
pertinent  molecular  parameters  must  be  determined 
and  tested  over  the  complete  range  of  P and  T of 
interest  for  the  dynamic  source. 


8.  Characterization  of  the  Dynamic  Source 

8.1  Accuracy  Levels  in  T and  P 

The  accuracy  level  for  dynamic  pressure  mea- 
surements using  the  proposed  optical/molecular 
approach  is  estimated  to  be  approximately  5%. 
This  estimate  is  based  on  our  ability  to  measure  the 
widths  and  shifts  of  the  D2  C -branch  transitions. 
For  the  mixed  gas  systems  proposed  as  standards, 
the  lines  shift  at  least  as  rapidly  than  they  broaden. 
Thus  our  accuracy  goals,  imply  that  we  can  locate 
the  line  center  to  approximately  l/10-l/20th  of  the 
line  width  and  that  we  can  make  line-width  mea- 
surements at  approximately  the  5%  level.  Achiev- 
ing these  measurement  accuracies  in  a single  shot 
experiment  is  one  of  the  development  goals  of  this 
work. 

Accuracy  in  temperature  measurements  depends 
on  our  ability  to  measure  the  relative  intensities  of 
spectral  features  associated  with  different  molecu- 
lar quantum  states.  In  practice,  this  comes  down  to 
comparing  an  observed  spectrum  to  predicted 
spectra  such  as  those  shown  in  figure  2.  In  a single 
shot  experiment  this  comparison  must  include  the 
effects  of  a finite-resolution,  spectroscopic  instru- 
ment function.  The  accuracy  of  these  comparisons 
for  single  shot  measurements  can  be  of  the  order  of 
5%  [3,4].  This  is  the  target  level  of  accuracy  for 
temperature  measurements  on  the  dynamic  source. 
The  simultaneous  and  independent  measurement  of 
pressure  will  help  eliminate  systematic  errors  in  ei- 
ther of  these  measurements. 

8.2  Temporal  Evolution 

It  is  important  to  recognize  that  the  nonlinear 
laser  diagnostic  techniques  which  allow  us  to  ob- 
tain a 10  ns  “snap-shot”  of  the  local  T and  P envi- 
ronment of  a molecule,  can  not  easily  be  extended 
to  provide  a continuous-in-time  measurement 
(“movie”)  of  the  evolution  of  T and  F in  a dynamic 
source.  These  techniques  do  provide  a primary 
standard  which  can  be  used  to  calibrate  measure- 
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ment  approaches  which  yield  good  relative  mea- 
sures of  T and  P with  rapid  and  continuous  tempo- 
ral response.  Until  this  point,  we  have  had  no 
primary  standard  which  could  assure  the  accuracy 
of  these  measurements.  The  primary  standard 
would  be  used  to  accurately  pin  the  P and  T values 
at  representative  points  in  the  temporal  evolution 
of  the  source  during  each  calibration  run.  Further, 
the  nonlinear  optical  techniques  would  be  used  to 
hold  the  absolute  calibration  of  the  source  over 
long  periods  of  use. 

9.  Conclusions 

We  have  described  a new  approach  to  the  mea- 
surement of  T and  P of  a dynamic  source.  The  mea- 
surement is  based  on  the  fundamental  properties  of 
molecules,  specifically  on  the  energies  and  popula- 
tions of  the  vibrational  and  rotational  levels  of  di- 
atomic gas  molecules.  The  nonlinear  optical 
technique,  coherent  anti-Stokes  Raman  spec- 
troscopy (CARS),  has  been  proposed  as  the 
method  for  acquiring  the  information  from  the 
molecular  system.  This  technique  has  the  advan- 
tages of  very  rapid  (nanosecond)  measurement 
times,  small  sampling  volumes,  and  a well  under- 
stood and  verified  T-  and  P-dependent  spectrum. 
The  development  of  this  measurement  system  to 
provide  a reference  standard  for  dynamic  calibra- 
tions of  T up  to  1500  K and  P up  to  10*  Pa  with  a 
5%  accuracy  appears  highly  feasible. 


10.  Appendix 

10.1  Raman  Scattering 

The  Raman  spectrum  [5]  arises  from  inelastic 
scattering  in  which  an  incident  photon  of  fre- 
quency Vo  is  “absorbed,”  a scattered  photon  of  fre- 
quency Vo±Av  is  “emitted,”  and  the  molecule 
undergoes  a transition  among  its  internal  states. 
The  energy  difference  between  the  scattered  and 
incident  photons,  hAv,  is  equal  to  that  between  the 
internal  states  of  the  scattering  molecule.  More 
concisely,  Raman  scattering,  as  we  will  consider  it 
in  this  discussion,  will  refer  to  a two-photon  transi- 
tion originating  in  the  ground  electronic  state,  with 
intermediate  states  being  far  off-resonance,  one- 
photon-allowed  electronic  transitions.  The  initial 
and  final  states  of  this  transition  are  molecular  vi- 
bration-rotation levels  in  the  ground  electronic 
state. 


10,2  Nonlinear  Optics 

The  fact  that  Raman  scattering  is  a two-photon 
transition  makes  it  essentially  a nonlinear  optical 
phenomenon.  Many  nonlinear  Raman  spec- 
troscopy techniques  simply  involve  impressing  (on 
the  molecules)  optical  fields  corresponding  to  the 
two  photons  of  the  Raman  transition.  In  the  fol- 
lowing we  will  briefly  discuss  Raman  scattering  in 
the  context  of  nonlinear  optics  [3,9].  In  this  theory 
the  response  of  a macroscopic  system  to  an  applied 
electromagnetic  field,  E,  is  written  in  a form  such 
as 

n CC  x^')£:  + x^'>£:£-hx®F  (1) 

where  11  is  the  polarization  density  induced  by  the 
action  of  the  fields  and  x*”^  is  the  n th  order  suscep- 
tibility (a  tensor  of  rank  « -|- 1).  The  relationship  in 
eq  (1)  is  intended  to  indicate  that  we  are  concerned 
with  vector  quantities  and  that  the  response  of  the 
system  involves,  in  general,  all  orders  of  the  elec- 
tric fields.  The  electromagnetic  fields  and  polariza- 
tion densities  appearing  in  this  expression  typically 
are  assumed  to  be  monochromatic  (angular  fre- 
quency o))  plane  waves  (wave  vector  k)  with  the 
representation 

E = ^ 1/2  F((rf,)exp(A:,-r  — oj,  r)-f  c.c.  (2) 

i 

(c.c.  denotes  complex  conjugate)  with 


where  n,  is  the  index  of  refraction  and  c the  speed 
of  light.  With  this  representation  of  the  fields  and 
polarization,  the  susceptibilities  are  written  as  func- 
tions of  the  frequencies  of  the  associated 
monochromatic  components  of  the  field.  Nonlinear 
Raman  effects  occur  when  the  frequency  differ- 
ences and  electric  field  directions  of  the  applied 
fields  can  drive  a Raman  transition. 

10.3  Coherent  Anti-Stokes  Raman  Spectroscopy 

In  particular,  the  nonlinear  Raman  effect  known 
as  coherent  anti-Stokes  Raman  spectroscopy 
(CARS)  arises  from  the  3rd -order  terms  with  elec- 
tric field  components  £)((yo)>  jE*(ct)o),  and  [£'/(o)as)]* 
(the  j,k,l  are  spatial  indices  and  the  “*”  denotes 
complex  conjugate)  such  that  a polarization  density 
at  frequency  oJas  with  wave  vector  is  generated 
with  the  conditions 
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oi^  — loio—o),  (4a) 

and 

K^  = 2ko-K  (4b) 

These  equations  are  essentially  statements  of  con- 
servation of  energy  and  momentum.  The  latter 
condition  indicates  that  a phase  matching  geometry 
must  be  established  between  the  incident  “pump” 
field  (&)o)  and  the  incident  “Stokes”  field  (ws)  in 
order  to  generate  a signal  beam,  the  anti-Stokes 
component  (coas).  The  corresponding  susceptibility 
has  2-photon,  Raman  resonances  when 

U)o  — ^^molecular  (5) 

with  Awmoiecuiar  being  a molecular  transition  angular 
frequency  (27rAv  = Aco^oiecuiar)- 


and  anti-Stokes  electric  field  vectors  in  parallel  ori- 
entations is  defined  by  the  following  relations  [10] 

fll(u^as)  ^ X^ni(  U)as5^0?U)o,  UJ5) 


[E,(a>o)]^[E.(caJ]* 

(8) 

(9) 

S(a))  = /'(G)-' Apo/ 

(10) 

G =G{oi)  = i{oi  I -oij}+N  W 

(11) 

I 

0 

3 

II 

3 

(12) 

{Oij)jj-  = 2TTC  V{J)  bjj- 

(13) 

{^po)jf  = [Peq(v  =1,7)-  Peq(  V = 0,7)]6y^ 

(14) 

10.4  Raman  ^-Branch  Resonance 

For  example,  the  fundamental  (i.e.,  v=0— »v  = l) 
vibrational  Q branch  is  composed  of  all  transitions 
from  the  ground  to  first  excited  vibrational  level 
without  change  in  rotational  quantum  number;  the 
molecular  transition  frequency  in  this  case  is 

AcUn^olecular  = ^ " ‘ [£  ( V = 1 , J)  - £ ( V = 0,/)] . (6) 

This  relation  yields  a multiple  line  spectrum  for  the 
Q branch  because  the  rotational  energy  level  spac- 
ing changes  with  the  vibrational  state.  The  transi- 
tion frequencies  for  the  vibrational  0 -branch  of  a 
simple  diatomic  molecule  such  as  N2,  CO,  or  D2 
can  be  expressed  in  terms  of  familiar  molecular 
constants  [5]  as 

v(J)^Vo-aJ{J+V)  + fiJ\J  + \f  (7) 

with  Vo  the  frequency  of  the  pure  vibrational  transi- 
tion between  the  rotationless  (J  = 0)  states  and  with 
the  tte  and  expressing,  respectively,  the  change 
in  the  moment  of  inertia  and  (lowest  order)  cen- 
trifugal distortion  with  vibrational  state.  In  keeping 
with  convention,  we  will  express  frequency  in  the 
spectroscopic  unit  the  wavenumber,  cm”',  which  is 
the  frequency  in  Hz  divided  by  the  speed  of  light  in 
cm/s  [v/c  =a)/(2'7rc)]. 

10.5  CARS  Susceptibility 

The  CARS  susceptibility  resonance  for  the 
isotropic  vibrational  Q branch  with  laser,  Stokes, 


IXrrP-r,  ] ' .... 

afi  Z(^(a,„-a)o)^^(a,y,  + a>o)i 

with  A' the  number  density,  / a column  vector  of  I’s 
(superscript  “t”  signifying  transpose),  I the  unit  ma- 
trix, o}j  and  Apo  diagonal  matrices  of  transition  fre- 
quencies and  (equilibrium)  population  differences 
respectively.  A nonresonant  contribution,  to 
the  third  order  susceptibility  (see  below)  is  in- 
cluded in  eq  (9).  Equation  (10)  is  a short-hand  ma- 
trix equation  for  the  frequency  dependence  of  the 
(Raman)  resonant  part  of  the  susceptibility  which 
involves  the  relaxation  matrix  W.  This  matrix  car- 
ries all  the  information  about  the  effects  of  colli- 
sions on  the  spectrum;  we  will  discuss  it  later. 

10.6  Raman  Cross-Section 

First  we  note  that  eq  (15)  is  an  approximate  ex- 
pression for  the  Raman  scattering  cross-section.  It 
explicitly  shows  the  two,  one-photon  dipole  matrix 
elements,  jUmn.  connecting  the  initial  (/={v=0,/}) 
and  final  (A={v  = l,y})  molecular  levels  via  inter- 
mediate electronic  states,  r [the  subscripts  {mn } 
stand  for  sets  such  as  {ri}  or  {fr}  in  eq  (15)].  This 
form  of  the  cross-section  is  used  because  all  pho- 
tons (laser,  Stokes,  and,  anti-Stokes)  satisfy  the  re- 
lation 

mn  (16) 

with 

oi„„^fi-'[E{m)-E{n)].  (17) 
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As  a consequence  of  these  conditions,  no  denomi- 
nator in  eq  (15)  is  near  zero,  i.e.,  we  are  far  from 
any  one  photon  resonance. 

The  quantity  inside  the  squared-magnitude  sym- 
bols in  eq  (15)  is  an  expression  for  the  optical  polar- 
izability matrix  element  for  the  transition  {//}.  The 
optical  polarizability  is  a 2nd-rank  tensor  quantity 
which,  for  the  systems  of  interest  here,  has  two 
irreducible  tensor  components.  These  are  a scalar, 
rotationally  invariant  part  of  rank  zero,  and  a trace- 
less, symmetric  2nd-rank  tensor  part,  known  as  the 
anisotropy.  Reference  to  the  “isotropic  Q branch” 
indicates  consideration  of  that  part  of  the  interac- 
tion associated  with  the  scaler  part  of  the  polariz- 
ability matrix  elements. 

The  squared-magnitude  of  the  matrix  element  of 
the  scaler  portion  of  the  polarizability  is  linearly 
dependent  on  the  vibrational  quantum  number  but 
essentially  independent  of  the  rotational  quantum 
number  [6].  As  a result,  the  contribution  of  an  indi- 
vidual Q{J)  transition  to  the  isotropic  Q branch  is 
almost  exactly  proportional  to  the  population  dif- 
ference in  eq  (14). 


10.7  Nonresonant  Contribution  to  x 

Reference  to  eq  (1 1)  will  show  that  the  contribu- 
tion of  the  Raman  resonance  to  the  susceptibility 
has  both  real  and  imaginary  components.  Far  off 
the  Raman  resonance,  when  |<u  — oiy  | > \ NWfj  \ for 
all  {J',J},  this  contribution  becomes  pure  real,  cf. 
eq  (9).  When  one  includes  all  the  possible  field 
combinations  implicit  in  eq  (1),  there  are  a large 
variety  of  far  off-resonance  two-photon  transitions 
which  add  together  to  provide  a constant-in-fre- 
quency real  susceptibility  called  the  nonresonant 
background.  This  is  the  source  of  in  eq  (9). 

10.8  CARS  Intensity 

The  polarization  density  in  eq  (8)  is  a source 
term  in  Maxwell’s  equations  which  are  solved  to 
yield  the  anti-Stokes  field.  For  copropagatmg 
(along  z),  plane  wave,  linearly  polarized  (along  x) 
pump  and  Stokes  excitation,  the  approximate  solu- 
tion is  a plane,  anti-Stokes  wave,  with  electric  field 
parallel  to  x(x  = l) 


with 


X‘=^^^=nxun  (20) 

and  with 

8k=K,,-k^  (21) 


expressing  the  phase  matching  condition.  The 
wave  vectors  on  the  RHS  of  eq  (21)  are  those  of 
the  CARS  polarization,  eq  (4b),  and  the  generated, 
radiative  CARS  field,  eq  (18).  The  intensity  (power 
density)  associated  with  this  field  is 


/(<^J 


167c*G)L/(o)P[/(tOo)]^lx^Pfsiii(5fa/2)y 


(22) 


Equations  (18)-(22)  involve  some  idealizations 
and  approximations,  such  as  plane  wave  excitation 
and  weak  field  limits,  which  can  easily  be  removed 
in  the  theory  or  to  some  extent  achieved  in  the 
laboratory,  or  both,  see  for  example  reference  [3]. 
These  equations  do  provide  a basis  for  understand- 
ing the  essential  features  of  the  signal  generation 
process,  namely,  an  intensity  proportional  to:  the 
squared  magnitude  of  a third  order  nonlinear  sus- 
ceptibility, x(  — (^0,  o>o,  —dis),  the  square  of  the 
pump  intensity,  the  probe  intensity,  and  a phase 
matching  factor. 


10.9  Spatial  Resolution  and  Phase  Matching 

For  systems  with  8k=^0,  e.g.,  low  density  gases, 
spatial  resolution  and  signal  enhancement  is 
achieved  simultaneously  in  the  copropagating  ge- 
ometry, discussed  above,  by  using  focussed  pump 
and  Stokes  lasers.  Gaussian  transverse-mode  lasers 
yield  nearly  Gaussian  anti-Stokes  beams  generated 
from  spatial  volumes  of  diameter  less  than  the 
pump  focal  diameter  and  with  axial  distances  of  the 
order  of  6zo,  with  Zq  the  familiar  confocal  parameter 
of  the  focussed  pump  beam  [3].  For  a 5-mm  diame- 
ter {do)  pump  focussed  by  a 20-cm  (f)  focussing 
lens,  the  approximate  dimensions  of  the  sample  vol- 
ume are  found  from  the  equations  for  the  focussed 
beam  diameter,  Iwq, 


ET  = ) exp(^as2  - ) 


(18) 


4^/. 

TT  do 


2wo=— "r  = 28x  lO-" 


cm 


(23) 


£'(a)asrZ)=  yCARs  and  for  the  confocal  parameter 

2 

sin(8^/2)  zo=^^=0.108  cm.  (24) 

bk/1  ^ ^ 


42 

204 


Volume  95,  Number  1,  January-February  1990 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


Spatial  resolution  can  be  enhanced  by  crossing  the 
pump  and  Stokes  beams,  for  example  as  shown  in 
figure  1.  The  requirements  of  phase  matching  for 
signal  enhancement,  c.f  eqs  (21)  and  (22),  result  in 
the  folded  BOXCARS  [11]  arrangement  shown  in 
the  figure.  This  arrangement  will  yield  a spatial  res- 
olution nearer  to  the  dimensions  2wo  and  Izq.  Since 
the  folded  BOXCARS  geometry  yields  good  spa- 
tial separation  of  the  anti-Stokes  beam,  this  usually 
is  the  spectroscopically  preferred  arrangement. 
However,  many  other  configurations  are  possible 
[2]. 


10,10  CARS  Spectrum 

The  spectrum  results  from  the  frequency  depen- 
dence of  the  squared  magnitude  of  X-  From  eqs  (9) 
and  (10),  we  can  write  the  general  form 

I X I ' = (;^nr)2  + 2 xk  X^’^  + Xr^  + Xr^  (25) 

with  the  subscript  R referring  to  the  (Raman)  reso- 
nant part  of  the  susceptibility  and  the  prime  denot- 
ing real  and  the  double  prime  denoting  imaginary 
part,  respectively.  If  we  limit  consideration  to  a 
single,  presumed  isolated  (our  specific  meaning  for 
this  term  will  be  clarified  below)  transition,  at  cuy, 
the  resonant  part  of  the  susceptibility  is 


(fa)  —COj)-i-l  Ty 

^<^^(co-cojy  + r/ 


(26) 


ment  of  the  relaxation  matrix  introduced  in  eq  (11). 
It  is  a measure  of  the  rate  of  collisionally-induced 
exchange  between  transitions  (see  below).  Equa- 
tion (28)  simply  states  that  if  the  rate  of  exchange 
between  transitions  is  much  smaller  than  the  fre- 
quency separation  of  the  transitions,  then  they  can 
be  considered  isolated.  Conversely,  if  collisions  are 
mixing  lines  at  a rate  commensurate  with  their  sep- 
arations, we  will  see  line  interferences.  These  are 
described  by  the  full  relaxation  matrix  equation. 
The  corrections  for  multiple  line  spectra  which 
arise  when  we  cannot  ignore  the  nondiagonality  of 
the  relaxation  matrix  in  eq  (11)  are  discussed  in  the 
following. 


10.12  The  Relaxation  Matrix 

10,12.1  Pressure  Shift  and  Broadening  In  gen- 
eral, the  diagonal  elements  of  the  relaxation  matrix 
describe  the  familiar  pressure  broadening  and  shift- 
ing of  spectral  lines  [12].  Thus,  the  resonant  fre- 
quency and  linewidth  in  eq  (26)  are 


fa)y=fiy  + 5yP  (29) 

[with  the  convention  that  the  quantity  fiy  is  the 
zero  pressure  value  derived  from  accurately 
known  molecular  constants,  see  for  example  eq  (7)] 
and 

ry=7yR  (30) 


with  r y the  linewidth  of  the  pressure  broadened 
line.  The  sum  of  the  last  two  terms  in  eq  (25)  yields 
a Lorentzian  lineshape,  i.e.. 


1 

-ojjy+r/- 


(27) 


The  presence  of  a nonzero  nonresonant  back- 
ground leads  to  asymmetric  lineshapes  via  the  sec- 
ond term  in  eq  (25).  The  first  term  yields  a uniform 
nonzero  signal  level  which  can  be  dominant  if  the 
concentration  of  the  resonant  species  is  low. 


10,11  Collisional  Line  Interference 

A transition  can  be  considered  isolated  if  the 
condition  expressed  by  the  relation 


I 


HA 


(fajy  -fajy) 


<^1 


(28) 


is  satisfied.  The  quantity  Hy-yis  an  off-diagonal  ele- 


Equations (29)  and  (30)  implicitly  define  pressure 
shift  and  pressure  broadening  coefficients,  8 and  y 
respectively. 

The  change  and  spread  of  the  observed  reso- 
nance frequency  are  associated  with  intermolecular 
collisions.  Encounters  with  other  molecules  alter 
the  binding  forces  and  therefore  the  resonance  fre- 
quencies of  the  optically  active  molecule.  In  the 
pressure  regimes  of  our  interest,  such  encounters 
are  limited  to  binary  collision  events  such  that  the 
duration  of  the  encounter  is  much  shorter  than  the 
time  between  collisions.  In  simplest  terms,  the  shift- 
ing coefficient  expresses  the  average  rate  of  colli- 
sion-induced phase  shift  of  the  oscillator.  Similarly, 
the  broadening  coefficient  expresses  the  rate  at 
which  the  oscillator  is  dephased  by  the  collisions. 

These  effects  are  expressed  as  proportional  to 
pressure  because  the  probability  per  unit  time  that  a 
collision  produces  a phase  shift  in  the  observed  os- 
cillation frequency  can  be  written  as  the  product 
A-Vav-cr  of  the  number  density,  N (molecule/cm^), 
the  average  velocity,  Vav  (cm/s),  and  a cross- 
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section,  cr  (cm^),  appropriate  to  shifting  or  broad- 
ening for  the  collision.  Such  an  expression  is 
analogous  to  that  for  the  collisional  momentum 
transfer  which  defines  the  pressure  in  the  gas. 

We  now  turn  to  consideration  of  the  off-diagonal 
elements  of  the  relaxation  matrix. 

10.12.2  Line  Mixing  and  State-to-State  Rates 
In  the  familiar  “impact  limit”  of  collision  physics 
(aspects  of  which  were  enumerated  in  the  preced- 
ing section,  see  also  ref.  [12]),  the  elements  of  W are 
proportional  to  appropriately  averaged  products  of 
two  scattering  (“5 -matrix”)  operators.  For 
molecules  such  as  Nj  and  CO,  it  has  been  shown 
that  there  is  very  weak  dependence  of  the  scatter- 
ing matrix  elements  on  the  vibrational  states.  This 
reveals  itself  in  very  weak  shifting  and  a line  broad- 
ening which  is  independent  of  vibrational  (e.g., 
1^2  vs  0— »1)  or  overtone  (e.g.,  0-^2  vs  0-^1) 
branch.  The  absence  of  strong  vibrational  state  de- 
pendencies in  the  S' -matrix  elements  yields  a relax- 
ation matrix  for  the  isotropic  Q branch  which  can 
be  expressed  entirely  in  terms  of  state-to-state  rates 
for  rotational  energy  transfer.  The  off-diagonal  ele- 
ment, Wfj  for  these  conditions,  is  simply  the  nega- 
tive of  the  inelastic  rate  for  a transition  from  J to 
In  terms  of  the  discussion  in  the  last  section,  we  can 
see  that  the  occurrence  of  an  inelastic  collision  ter- 
minates the  oscillation,  i.e.,  equivalently  totally 
dephases  the  oscillator.  Very  little  shifting  would 
be  expected  under  these  conditions.  Thus,  the  diag- 
onal term,  the  line  broadening  coefficient,  reduces 
to  the  total  inelastic  rate  out  of  a state  J . 

Because  of  the  nonzero  off-diagonal  compo- 
nents, the  spectrum  resulting  from  the  relaxation 
matrix  equation  will  show  the  phenomenon  of  col- 
lisional narrowing  of  the  Q branch  with  increasing 
pressure.  The  important  criterion  for  the  narrow- 
ing is  a rapid  rate  of  rotational  energy  transfer  rela- 
tive to  the  frequency  difference  between  the  QiJ) 
transitions,  i.e.,  the  relation  in  eq  (28)  is  not  satis- 
fied. Systems  such  as  N2  and  CO  show  significant 
narrowing  at  1 atm  and  room  temperature,  because 
they  have  relatively  small  rotational  energy  gaps 
with  concomitantly  large  rates  of  rotational  energy 
transfer,  and  small  interline  spacings  within  the  Q 
branch.  In  contrast,  the  spectra  of  D2  displayed  in 
figure  4 can  be  considered  isolated,  independent 
lines.  We  illustrate  the  effects  of  interference  in  the 
following. 

10.13  ^-Branch  Spectra  of  N2 


line  model  (fVj-j  = 0 for  for  the  spectrum  at 

10  atm  is  also  included  in  the  figure.  Clearly,  ne- 
glect of  the  collisional  interference  terms,  Wj-j, 
would  lead  to  large  errors  in  pressure  determina- 
tions. The  accuracy  of  theoretically  predicted  spec- 
tra, such  as  those  shown  in  figure  5,  has  been  tested 
against  spectra  measured  on  samples  in  the  pressure 
range  up  to  100  atm  and  for  temperatures  from  295 
up  to  2400  K (the  upper  limit  in  T depends  on  P 
because  of  the  pressure  limits  of  high-T  furnaces 
and  cells),  see  reference  [4]. 

10,14  Energy  Gap  Laws  for  Rotational  Energy 
Transfer 

The  fundamental  basis  for  calculation  of  the 
spectra  shown  in  figure  5 is  an  energy-gap  rate  law 
[13]  for  the  individual  state-to-state  rates  for  rota- 
tional energy  transfer.  As  discussed  above,  knowl- 
edge of  the  state-to-state  rates  allows  specification 
of  the  entire  relaxation  matrix  and  therefore 


Figure  5.  Calculated  CARS  spectra  for  the  pure  Nj  Q branch  at 
295  K.  The  intensity  at  each  pressure  has  been  multiplied  by  the 
indicated  factor  and  plotted  on  a unit-normalized  scale.  An  iso- 
lated line  model  which  considers  only  the  diagonal  elements  of 
the  relaxation  matrix  is  included  for  comparison  at  10  atm. 


Pressure  dependent  spectra  of  pure  N2  at  room  T 
are  shown  in  figure  5.  For  comparison,  an  isolated 
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computation  of  the  spectrum  [from  eqs  (9)-(l  1)  and 
(22)].  The  most  widely  used  form  of  the  rate  law 
for  specification  of  the  relaxation  matrix  is  the 
modified  exponential  gap  (MEG)  law  [14]  which 
parameterizes  the  upward  with  J <J')  state- 

to-state  rate  according  to 


^fj  = aof(Ty\ 


\+A  EjjykTb' 
1 +A  EiJykT 


(UrAEyymv  /-/3{£(/')-£(/)}\ 
\\+AE{J)/kT)  y kT  ) 

(31) 

with 


= l-exp(-w) 

> \Tj  l-exp(-mr/ro)‘ 


(32) 


The  downward  rate  is  found  from  the  principle  of 
detailed  balance,  i.e., 


Pe,(/)  = Wfj-  (33) 

The  quantity  A is  calculated  from  the  Lennard- 
Jones  intermolecular  potential  parameters  for  each 
molecular  collision  system.  The  parameters  ao,  6', 
8,  /S,  and  m are  determined  from  comparisons  to 
accurately  measured  line  broadening  coefficients, 
jj,  (functions  of  J and  T)  and  to  experimentally 
determined  pressure  dependent  spectra. 

There  remain  some  significant  research  issues  in 
the  prediction  of  N2  isotropic  0-branch  spectra, 
such  as:  some  refinement  of  the  rate  law  model 
seems  important  at  high  temperature  [15]  and  the 
appropriate  spectral  and  rate  law  model  for  situa- 
tions where  N2  is  found  in  a foreign  gas  host  must 
be  experimentally  examined. 


10.15  Q-Branch  Spectra  of  the  Hydrogens 

The  state-of-the-art  for  prediction  of  the  spec- 
trum of  molecular  systems  such  as  D2(or  H2):X 
[X  = D2(or  H,),  He,  Ar,  NJ  is  not  as  fully  devel- 
oped as  for  N2:N2.  This  situation  arises  because 
there  is  a large  vibrational  state  dependence  in  the 
5 -matrix  elements.  Such  a vibrational  state  depen- 
dence leads  to  the  addition  of  a significant  “elastic 
vibrational  dephasing,”  discussed  earlier,  contribu- 
tion to  the  linewidth  and  line  shift  in  these  systems. 
The  elastic  vibrational  dephasing  component,  al- 
though understood  formally  from  theory,  is  not 
easily  calculated  and  does  not  have  a simple  predic- 
tive basis  similar  to  that  just  presented  for  the  rota- 
tionally  inelastic  component. 


The  significance  of  elastic  dephasing  for  the  hy- 
drogens is  Judged  relative  to  their  rotationally  in- 
elastic contributions.  For  the  hydrogenic  systems, 
the  presence  of  elastic  vibrational  dephasing  is  evi- 
denced by  pressure  shifting  coefficients,  8y  [cf  eq 
(29)],  which  are  of  the  order  of  and  often  a few 
times  larger  than  the  line  broadening  coefficients, 
7j  [cf  eq  (30)].  Because  of  the  large  rotational  en- 
ergy gaps  in  these  systems  the  inelastic  portion  of 
the  total  line  broadening  is  relatively  weak  (espe- 
cially for  the  foreign  gas  collision  partners).  Small 
inelastic  rates  and  large  vibrational  dephasing  both 
lead  to  a decrease  in  the  line  interference  resulting 
from  the  solution  of  eqs  (10)  and  (11).  Finally,  there 
is  large  rotation-vibration  interaction  in  the  hydro- 
gens, e.g.,  tte  in  eq  (7)  is  1.08  for  D2,  more  than  60 
times  larger  than  the  of  N2.  As  a result  of  all 
these  factors,  the  hydrogens  will  display  essentially 
isolated-line  type  spectra  with  shift  and  width  mea- 
surements yielding  good  measures  of  the  pressure. 
The  spectra  displayed  in  figure  4 illustrate  the  pres- 
sure sensitivity  of  the  pure  gas  system,  those  in  fig- 
ure 6 indicate  the  potential  sensitivity  for  the  D2:Ar 
system.  The  spectra  given  in  this  figure  are  based 
on  actual  measurements  which  to  this  point  have 
been  extended  only  to  50  atm.  Realizing  the  poten- 
tial of  this  system  requires  more  study  as  we  will 
discuss  in  the  following. 


Figure  6.  Calculated  CARS  spectra  for  the  7=0  and  7=  1 tran- 
sitions of  D2  contained  at  10%  concentration  in  Ar  at  295  K. 
The  spectra  are  based  on  experimental  measurements  which  ex- 
tend only  up  to  50  atm.  The  extrapolation  to  100  atm  is  based  on 
the  lower  pressure  data.  See  caption  and  compare  to  figure  4. 
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The  complete  J,  P,  and  T dependence  of  the 
spectra  of  D2(H2):X  is  yet  to  be  determined.  Some 
preliminary  data  on  pure  H2  up  to  1000  K and  50 
atm  has  been  obtained.  An  empirical  fitting  law 
based  on  the  MEG  law  for  the  rotationally  inelastic 
contribution  and  a simple  power  law  in  T for  scal- 
ing the  (assumed  J independent)  vibrational 
dephasing  contribution  has  been  applied  with  some 
success  to  these  data.  More  work  is  needed  for  all 
the  systems  of  interest,  in  order  to  have  a predic- 
tive basis. 

In  particular  there  remain  some  important  scien- 
tific questions  with  respect  to  the  foreign  gas  sys- 
tems. Recent  work  has  revealed  anomalous 
concentration  dependent  lineshapes  for  the  D2(and 
H2):Ar  system  [8].  These  anomalies  are  associated 
with  a large  (collision)  speed  dependence  in  the 
cross-sections  for  line  shifting  and  a small  propen- 
sity of  D2:Ar  collisions  to  result  in  speed  changes. 
As  a result,  the  lines  are  inhomogeneously  broad- 
ened and  rather  asymmetric  at  very  low  concentra- 
tions of  the  active  molecule  in  the  Ar  host.  More 
study  of  the  concentration,  pressure,  and  tempera- 
ture dependence  of  the  D2:Ar  lineshape  is  required. 

10.16  Doppler  Broadening 

The  line  broadening,  shifting,  and  interferences 
discussed  above  are  appropriate  to  pressure  and 
temperature  regimes  in  which  the  rate  of  velocity 
changing  collisions,  Nv^^cti^,  with  ctd  the  “diffusion 
cross-section,”  is  larger  than  the  Doppler  width, 
~ I /cq  — /:s|  Vav,  of  the  Raman  transition.  Doppler 
broadening  arises  from  the  free  streaming  of  the 
molecules,  and,  at  low  pressures,  it  can  be  used  for 
measurement  of  translational  temperature.  For  the 
hydrogens,  in  particular,  this  source  of  line  broad- 
ening is  significant  at  pressures  up  to  10  atm  at  high 
temperatures,  e.g.,  1000  K.  Research  on  the  D2:X 
systems  (X=:D2,  H2,  He,  CH4,  Ar),  under  condi- 
tions where  both  velocity  changing  and  normal 
pressure  broadening  collisions,  as  discussed  above, 
contribute  to  the  line  shape  is  in  progress.  A reli- 
able theoretical  description  appears  possible,  espe- 
cially for  relatively  light  collision  partners,  using  a 
simple  diffusion-like  model  solution  to  kinetic  equa- 
tions for  line  formation.  However,  more  work,  in 
particular  for  the  heavy  collision  partners  and  high 
temperatures,  remains  to  be  done  for  these  systems. 

10.17  Summary 

In  the  above,  we  have  briefly  summarized  essen- 
tial aspects  of  the  theory  of  nonlinear  Raman  spec- 


troscopy as  seen  from  the  viewpoint  of  one  who 
wishes  to  apply  the  technique  to  temperature  and 
pressure  measurements.  This  has  involved  some 
discussion  of  elements  of  the  theory  of  spectral  line 
formation  for  the  relevant  Raman  transitions.  Al- 
though the  level  of  detail  and  rigor  of  this  discus- 
sion has  been  limited,  the  intention  is  simply  to 
indicate  that  we  have  a good  understanding  of  the 
basic  physics  which  underlies  our  measurement  ap- 
proach. A great  deal  of  work  has  gone  into  build- 
ing this  understanding  and  developing  a reliable 
predictive  basis  for  temperature  and  pressure  mea- 
surements. It  is  now  possible  to  transfer  the  static 
pressure  and  temperature  scale  to  the  molecular 
level  to  serve  as  a primary  standard  for  dynamic  P 
and  T. 
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